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ABSTRACT 



This thesis develops a Lanchester-type force-on- force combat model 
simulating small-unit amphibious operations. The model commences with 
a ship-to-shore assault of aggressor forces mounted onboard Landing 
Vehicle Assault craft moving against a defensive force ashore. Once 
the ship-to-shore phase of combat is completed, the model continues 
to simulate land combat further inland between the assaulting aggressor 
forces and other defensive forces occupying key terrain. 

The main thrust of the thesis is to alleviate some of the problems 
associated with the inherent abstractness of Lanchester-type combat 
models; specifically, to develop "user-friendly 1 ' input-data and output 
structure, and more thorough documentation of the model's algorithms 
to provide a model which would be more easily understood and utilized 
by students of combat modeling. 
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I. INTRODUCTION 



A. OVERVIEW 

This thesis develops a Lanchester-type force-on-force combat model 
simulating small-unit amphibious operations. The model commences with 
a ship-to-shore assault of aggressor forces (e.g., a U.S. Marine Infantry 
Battalion), mounted onboard Landing Vehicle Assault craft (LVA) moving 
against a defensive force ashore located in fixed positions along the 
coast the aggressor force is attempting to occupy. Once the ship-to- 
shore phase of combat is completed, the model continues to simulate 
land combat further inland between the assaulting aggressor forces and 
other defensive forces occupying key terrain. 

The main thrust of the thesis is to alleviate some of the problems 
associated with the inherent abstractness of Lanchester-type combat 
models (see [Ref. 1]), and specifically to integrate and enhance work 
done in previous models, to develop "user-friendly" enhancements, and 
more thorough documentation of algorithms to provide a model which would 
be more easily understood and utilized by students of combat modeling. 

B. BACKGROUND AND GENERAL MODEL 

1 . Overview 

The small -unit 'amphibious operation combat model presented in 
this thesis is the result of the integration and enhancement of two 
independent combat models. The first model is a ship-to-shore combat 
model which models a ship-to-shore assault conducted by landing vehicle 
assault craft against fixed enemy positions ashore. The second model 



9 



is a land combat model which models a land assault conducted by LVA 
forces on a beach, against fixed enemy positions further inland. 

2. Original Ship-to-Shore Combat Model 

The ship-to-shore combat model used as a basis for this thesis 
was presented in a thesis by David L. Chadwick [Ref. 2]. It modeled 
the amphibious assault of five waves of LVA against a defensive force 
composed of tanks and antitank guided missiles (ATGM) in fixed positions 
ashore. Attrition was modeled using Lanchester area-fire and aimed-fire 
equations. The purpose of developing such a model was to determine tne 
optimal design characteristics of LVA in an amphibious assault for a 
given combat scenario. The optimal design of an LVA was considered 
to be that design which produced the lowest level of LVA attrition for 
the given combat scenario. 

3. Original Land Combat Model 

The land combat model used as a basis for this thesis was 
developed in Joseph Smoler's thesis [Ref. 3]. It modeled land combat 
conducted by three aggressor force units utilizing tanks assaulting 
three defensive force units armed with Tube-Launched, Optical -Guiaed, 
Wire-Controlled missiles (TOW's) in fixed positions. The location of 
the land combat was the Fulda Gap region in West Germany. Attrition 
was modeled using Lanchester aimed-fire equations. The purpose of 
Smoler's thesis was to develop a basic small-unit land combat model for 
determining optimal defensive unit locations for a given combat scenario. 
The optimal locations of the defensive units were considered to be those 
locations which provided the lowest level of attrition of the defensive 
units for the given combat scenario. 
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4. The Enhanced Land Combat Model 



An enhanced version of Smoler's land combat model was developed 
by Glenn Mills in his thesis [Ref. 4]. The enhancements developed by 
Mills added flexibility to Smoler's land combat model by providing 
user selected options which could be employed depending upon the abili- 
ties and desires of the model's user. The enhancements included the 
option of altering the aggressor force's attack routes enabling the user 
to study not only the optimal defensive unit locations, but the optimal 
aggressor force attack routes for the given defensive unit locations as 
wel 1 . 

A second enhancement was the option of selecting a stochastic 
attrition-rate coefficient. This introduced the element of randomness 
into the model's attrition algorithm providing a more realistic approach 
to modeling a unit's fighting effectiveness. 

The third enhancement is the option of providing alternate 
defensive positions so that the defensive units could move to more 
defensible terrain once their original positions had become untenable. 

5. The Original Small -Unit Amphibious Warfare Model 

The original small -unit amphibious warfare model used as a 
basis for this thesis was developed by Soon Dae Park in his thesis 
[Ref. 5]. Park's model attempted to conceptualize the flow of events 
of an amphibious assault by first running the ship-to-shore model, 
followed immediately by running the land combat model. The analysis 
of this model as a class project served as the catalyst for the develop- 
ment of the small -unit amphibious operation combat model presented in 
this thesis. 
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6. The Analysis of Park's Model 



The class project conducted by Clay Grubb, Robert Larson, and 
this author had as its purpose the analysis of Soon Dae Park's small - 
unit amphibious operation combat model. The results of the analysis 
revealed the value of Park's thesis in providing a general scheme of 
events for the modeling of small -unit amphibious operations. The results 
also identified some enhancements that could be applied to his conceptua- 
lized model that would integrate the ship-to-shore and land combat models 
into a singular small-unit amphibious operation combat model. The develop- 
ment and application of these enhancements to Park's model served as the 
foundation for this thesis, and the development of the model presented. 

C. MAJOR GOAL AND OBJECTIVES OF THE THESIS 

1 . Major Goal of the Thesis 

The overall goal of the thesis is the development of a small-unit 
amphibious operation combat model. It will be based on the integration 
and enhancement of the two combat models discussed in the previous section 
of this chapter. There are three underlying objectives of the thesis 
which will guide the development of the model toward the accomplishment 
of this goal . 

2. Objectives of the Thesis 

a. Integration of Independent Combat Models 

The first objective in the development of the model was to 
integrate two initially independent combat models into a singular con- 
tinuous flow combat model. This was accomplished by first allowing 
force levels at the completion of the ship-to-shore phase of combat to 
be used as the initial force levels in the land phase of combat. 
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Secondly, it was recognized that four combat modelers con- 
tributed to the resulting model presented by Park in his thesis. As 
such, four individualized FORTRAN coding techniques were reformulated 
into one style to provide a more tractible small-unit amphibious oper- 
ation combat model . 

b. User-Friendly Input-Data and Output Structure 

The second objective of the thesis was to provide a user- 
friendly combat model. It is a major contention of this thesis that 
combat modelers have not adhered closely to this principle when pro- 
viding combat models for the United States military. Furthermore, it 
is believed that the lack of concern given to this approach of combat 
modeling is a major reason for the less than unanimous reception that 
combat models have received by the United States military as tools for 
training its commanders aid staffs. Therefore, the model presented in 
this thesis was designed and documented with the user's capabilities 
and needs in mind as opposed to those of the programmer. 

c. Student-Oriented Combat Model 

The third objective of the thesis was to provide the student 
of combat modeling with a combat model which was easily understood and 
studied. As a result, the model presented in this thesis was designed 
with a low level of complexity to allow the student with little or no 
experience in combat modeling to understand more easily the combat 
modeling theory and its application. 
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II. MODEL ENHANCEMENTS 



A. OVERVIEW 

This thesis had as its goal the development of a small-unit amphibious 
operation combat model. Guided by the three objectives discussed in 
Chapter One, five modeling enhancements were aoplied to the two original 
combat models serving as the foundation for the resulting small-unit 
amphibious operation combat model presented in this thesis. The enhance- 
ments provide for the proper integration of the ship-to-shore and land 
combat models. In addition, they have contribjted to the development 
of a more user-friendly combat model which can be used to assist combat 
modeling students in their understanding the theory of combat modeling 
and its appl ication. 

B. INTEGRATION OF SHIP-TO-SHORE AND LAND COMBAT MODELS 

The intent of the model presented here is to view the amphibious 
assault as a continuous process made up of two phases of combat (ship- 
to-shore, and land combat) where the land combat phase is dependent 
upon the outcome of the ship-to-shore combat pnase of the model. 

Implementation of this enhancement called for the creation of a 
new variable, Total Landing Force Ashore (TLF), which would accumulate 
the surviving landing force of each assault wave as it reached the beach. 
This total landing force ashore would than be redistributed into three 
main assault units for the land combat phase of the model. The 
rationale for the redistribution of forces is based on realistic 
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military doctrine which is to maintain a well-balanced force when the 
strength and location of the enemy is unknown to the assaulting forces 
(as is assumed in the model). 

Since the manner in which defender force levels are determined by 
the ship-to-shore and land combat models appears to be quite realistic, 
the defending force level as modeled by Soon Dae Park was used as input 
to the land combat phase. In particular, if the aggressor force had been 
successful in routing the defending forces situated on the beach, 
defending forces situated further inland naturally would be impelled 
to defend the remaining terrain still in their possession, it should 
be noted that the size of these defending forces further inland is in 
option of the user which in itself can be varied for analysis of variant 
battle scenarios. 

C. AGGRESSOR FORCE ATTRITION— SHIP-TO-SHORE PHASE 

Attrition in Lanchester-type combat modeling is based upon the 
expected percentage of the original force remaining at a given point 
in time. The expected percentage of forces remaining then car be 
restated in terms of a real number to represent the expected number of 
forces remaining. This method of computing reduced force levels is 
considered to be quite appropriate when modeling land combat, and was 
implemented by Chadwick in his ship-to-shore combat model to simulate 
LVA attrition. However, use of Lanchester equations to model such 
vehicular attrition of a vehicle at sea was determined to be inappro- 
priate. Where it is a reasonable assumption that a disabled vehicle 
on land still can contribute something toward the final outcome of 
the battle if any of its weapons systems or onboard troops survive. 
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an LVA that is disabled at sea is of no use to the amphibious assault 
and subsequent land combat phase. The LVA will be recovered, and on- 
board troops brought to the landing site after the assault has taken 
place. 

Chadwick was not concerned with this distinction due to his model's 
purpose of modeling LVA attrition in terms of ship-to-shore movement 
only. Therefore, he simply utilized Lanchester equations in modeling 
LVA attrition resulting in fractionalized losses of assaulting LVA's. 
However, if a ship-to-shore combat model is to be properly integrated 
with a land combat model, only whole numbers of LVA's ashore should be 
used as input. Hence, an enhancement was made to the model. 

The approach was to find the integer value of the number of surviv- 
ing LVA's in each assault wave, and then sum these values resulting 
in the total landing force ashore (TLF) . The fractional portion remain- 
ing was considered to be those LVA's disabled at sea and unable to 
participate in the land combat phase of the operation. 

D. STOCHASTIC ATTRITION-RATE COEFFICIENT MODIFICATION 

Mill's land combat model allowed the user the option of selecting 
either deterministic or stochastic attrition-rate coefficients to be 
used in assessing the attrition of opposing forces. The justification 
for utilizing stochastic attrition-rate coefficients to model force-on- 
force attrition rates was based upon the assumption that the attrition- 
rate coefficient is a random quantity measuring a unit's fighting 
ability, and can be estimated before any given battle. 

This can be illustrated by considering the expected value of a 
random variable. For example, assume a probability distribution is 
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selected for the random variable such that the expected value of the 
random variable is equal to the deterministic attrition-rate coefficient 
set for all units. When a random sample is taken from this distribution, 
the individual values assigned the random variable will serve as 
individual unit attrition-rate coefficients, where the sample mean 
will serve as the overall force attrition-rate coefficient. The result 
is that the overall force attrition-rate is equal to the sample mean, 
which is approximately equal to the population mean of the random 
variable. Recalling that this population was selected with a mean 
that equalled the deterministic attrition-rate coefficient, units now 
have their own individual attrition-rate coefficients, while the force 
attrition-rate coefficient has remained close to the intended value of 
the deterministic attrition-rate coefficient. This is more realistic 
than the deterministic option since each unit would be expected to 
have a different level of effectiveness, which necessarily would imply 
different attrition rates while maintaining one overall force level 
attrition rate. 

The attrition-rate coefficient, A.jj, is used as the measure of the 
rate a firer in Unit i attrits a target in Unit j. This has been 
likened to the fighting effectiveness of a particular Unit i. Obviously, 
this is a variable quantity influenced by a myriad of factors to include 
esprit de corps, past history of success or failure, prior exposure to 
combat, weather, quality of leadership, etc. The intent of such a 
basic model as this is to attempt to capture the overall effect of 
these factors by developing a distribution of a unit's initial fighting 
capabilities (specifically, to develop a distribution of A — ' s for the 
unit) . 
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Mills proposed a distribution based upon a quadradic function which 
would produce a symmetric distribution with a mean value of approximately 
0.55. This distribution restricted a unit's maximum effectiveness to 
only 80 percent of its maximum capable effectiveness level. It also 
implied that the average unit in combat will only perform at 55 percent 
of its maximum effectiveness level at any given time. 

A more plausible way of assigning a distribution to the A--'s might 

* J 

be a truncated Normal Distribution limited to values between 0.00 and 
1.00. However, this approach would leave little flexibility in terms 
of modeling variant scenarios since the opposing forces always would 
have attrition-rate coefficients associated with that particular dis- 
tribution whenever the stochastic option was selected. This restriction 
is due to the programming constraints encountered in attempting to 
implement variant truncated Normal Distributions in the model. There- 
fore, a Beta Distribution was selected for use in the model. 

The natural range of the Beta Distribution is from 0.00 to 1.00 
thereby alleviating the burden of constructing a truncated distribution. 
Furthermore, its two scaling parameters, P and Q, can be selected readily 
and input by the user to construct virtually any variant of the Beta 
Distribution so desired. The specific values selected for P and Q would 
parameterize the distribution of the A^.'s according to the user's par- 
ticular combat scenario without the burden of reprogramming the distri- 
bution on each successive run of the model. 

The density function for the Beta Distribution is as follows: 
f(x) = X P ' 1 (1-X) Q_1 for 0.0 i x i 1.0 



Therefore, a P=21 and Q=7 would yield a distribution of A-.'s with a 

1 J 

mean of 0.75. This says that a unit with an A., of 0.75 is operating 
at 75 percent of its potential effectiveness. Whereas, a P=7 and Q=21 
would yield a distribution of A^’s with a mean of 0.25, indicating that 
a unit is operating at 25 percent of its potential effectiveness. 

To illustrate the flexibility of this approach in determining sto- 
chastic attrition-rate coefficients. Figure 2-1 is provided displaying 
the distribution of A—'s that would be obtained when the user alters 
the parameters of the Beta Distribution. The user now can model a 
strong elite force using, for example, parameter values P=21 , Q=7 , or 
model a weak and poorly lead force using parameter values P=7 and Q=21 , 
depending upon the particular battle scenario the user is analyzing. 

While the Beta Distribution used in this thesis is different than 
the Quadradic Distribution used by Glenn Mills, the implementation of 
this distribution for the attrition-rate coefficients is exactly the 
same as originally modeled. Since it was assumed earlier that the fight- 
ing effectiveness of each unit is a random quantity prior to a given 
battle, it is only necessary to obtain a realization of the random 
variable for each unit prior to the initialization of the battle. This 
realization, A°,. , is determined by the user-supplied inputs P and Q, 

^ j 

and subsequent calls to a Beta Distribution Random Deviate Generator 
[Ref. 6]. Therefore, an attrition-rate coefficient is computed for 
each unit using the following equation: 

A?, x (1 - r/r ) 2 for 0 i r ± r 
ij ' e e 

A ij 

, 0 for r g ± r 
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BET8 DENSITY 



f(x) = X P_1 (1 - X) 2 * 1 







X X 



Figure 2-1. Beta Density 
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Where: A°j = Realization of unit’s fighting effectiveness 

r = Current range between firer and target 
r e = Maximum effective range of a firer' s weapon 
This function was utilized because it is a function of both range and 
thus creating a different attrition-rate curve for each unit, 
depending on that unit's effectiveness level prior to the battle. A 
graphic illustration of an attrition-rate coefficient curve for an 
A°j equal to a mean of 0.75 from the Beta Distribution where P=21 and 
Q=7, and the maximum effective range r , of 3000 meters would look 
like Figure 2-2. 




Figure 2-2. Attrition-Rate Coefficient Curve for A?. = 0.75 

and r = 3000m 1J 

e 

To illustrate the effect that this stochastic attrition option has 
on the outcome of the model, two runs of the model were made using this 
option, while varying the Beta Distribution parameter values for both 
forces on each run. All other characteristics of both forces were left 
unaltered. In the first run, the aggressor forces were modeled to 
operate at 75 percent of their potential effectiveness, and the defending 
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forces were modeled to operate at 25 percent of their potential effec- 
tiveness. The battle outcome, as listed in Table 2-1, indicates that 
the aggressor forces won the battle. In the second run, the potential 
effectiveness of the opposing forces was reversed. The aggressor forces 
were now modeled to operate at 25 percent of their potential effective- 
ness, and the defending forces were modeled to operate at 75 percent of 
their potential effectiveness. The battle outcome, as listed in Table 
2-1, indicates that the battle was terminated due to the opposing forces 
being too close. The aggressor forces were unable to overrun the defend- 
ing forces, as was the case in the first run, which was due solely to the 
change in the potential fighting effectiveness of the opposing forces. 

A change in the battle outcome was expected; however, to what degree 
that change would be was unknown. The fact that the defender forces 
were unable to win the battle on the second run, while having a much 
higher level of effectiveness, indicates that other characteristics of 
the opposing forces were also playing an important role in the battle 
(e.g., types of weapons employed, original force levels, speed of attack, 
etc . ) . 

Through the use of the stochastic attrition option, the user now 
has the capability of studying one more facet of combat (i.e., potential 
fighting effectiveness), and can analyze to what degree different fight- 
ing effectiveness levels will have on final battle outcome. 

E. USER-FRIENDLY I/O STRUCTURE 

A significant and important part of writing a computer program for 
a combat model is to provide for the input and output of data to and 
from the program. It is my belief that one of the major factors 
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contributing to the lukewarm reception, in general, that combat models 
have received by the United States military is due, in part, to the 
poorly designed input-data and output structure of the combat models. 

The primary user of those models, the military commander, normally 
finds it difficult, to decipher the myriad of input-data requirements , 
or the voluminous output from combat models that supposedly were de- 
signed for the commander's use. It is a contention of this thesis that 
if more attention was given to the development of user-friendly input- 
data/output requirements, that more interest would be generated toward 
the use of such models in training military commanders. Therefore, an 
enhancement was made to the input-data and output requirements of the 
model to demonstrate a method of alleviating this problem. 

1 . User-Friendly Input Structure 

A readymade input data file was provided with the model to serve 
as a guide for entering all of the required data in the correct format 
required by the model (see Appendix C). Each variable requiring input 
for the model has been listed in the sample input file with sufficient 
space provided for ensuring that data is entered in the correct format. 
This file, therefore, provides the unfamiliar user of the model with the 
opportunity to utilize the model with only a limited knowledge of the 
model's algorithm and input requirements. This type of user-oriented 
input requirement will alleviate some of the apprehension that an unfami - 
iar user of the program might have, and might actually act as a catalyst 
in increasing the amount of use the model receives. 
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2. User-Frieifdly Output Structure 

Indecipherable output, or too much output from a model, can be 
just as much of a deterrent to a model's use as complex input require- 
ments can be. This point was brought out by Ye S. Venttsal in her dis- 
cussion of good combat models: 

It is advisable in such "training" modeling of combat actions 
that the commander receive information from the computer not in 
the form of mean characteristics averaged over a set of realiza- 
tions, but rather in the form of only one specific realization, 
on the basis of which a decision is in fact made. [Ref. 7] 

To paraphrase Venttsal, the combat model output must be clear, 
concise, and identifiable to the military commander. Furthermore, it 
must answer the questions that were originally asked by the user--speci- 
fically, who won and why? 

The model output was therefore restructured to provide a concise 
listing of what input parameters were entered into the program for pro- 
cessing, and a concise and understandable output summary of what occurred 
throughout the battle (see Appendix F). Additional" y, a new feature was 
introduced into the model which gives the user the option of viewing 
either a detailed time-step battle summary, or just a final battle summary 
of what occurred in the running of the model. 

F. DOCUMENTATION AND PROGRAM FORMAT 

Two of the objectives of the model presented in this thesis were 
first to serve as an example of the way in which combat models should 
be designed to be user-friendly to ensure their acceptance and use in 
training military commanders; and secondly, to serve as a model for 
combat-modeling students so that they might acquire a better understand- 
ing of how combat models ought to be programmed into a computer. 
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It already has been discussed how the user-friendly I/O structure assists 
the user of the model . However, proper structuring of programs for read- 
ability and good documentation is equally necessary to ensure readability 
and understanding by students and analysts. 

The FORTRAN program presented in this thesis which integrated and 
enhanced the ship-to-shore and land combat models is an amalgamation of 
subroutines originally written by different people, with their own unique 
style of programming. The interweaving of these four styles of program- 
ming throughout the program seriously detracted from the smooth flew of 
program structure and readability desired when analyzing the computer 
program. Therefore, an enhancement was made to the model: the program 
was restructured so that it would follow one basic style of programming 
(see Appendix B). New labeling and structuring of formatted statements 
and nested FORTRAN functions were provided to make the computer program 
more readable. This restructuring should assist the student in under- 
standing the program flow, and provide an incentive to those interested 
students to develop future enhancements to the model . 

In addition to developing one style of programming, more detailed 
documentation of variable definitions and descriptions of program flow 
were added to the program. The purpose of this documentation was to 
have the program serve as a reference to itself in order that the reader 
would not be forced to refer to various manuals outside of the program 
each time an explanation of the functioning of a particular aspect of 
the program is desired. 
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III. CURRENT MODEL DESIGN 



A. OVERVIEW 

The small -unit amphibious operation combat model presented in this 
thesis consists of the integration and enhancement of a ship-to-shore 
combat submodel, and a land combat submodel. Both of the original sub- 
models were similar in design, basing force attrition on Lanchester- 
type expected-value equations. As presented earlier, enhancements to 
both submodels reduced the differences in design of these submodels, 
molding them into what may be called a small -unit amphibious operation 
combat model. Figure 3-1 provides the scheme for the sequence and 
general flow of events in the overall model. 

It should be noted that although the ship-to-shore and land combat 
models are quite similar, they still have their own unique characteristics 
in modeling certain events that take place throughout the battle. There- 
fore, in discussing the model as a whole, the two phases of the battle 
will be addressed separately, and those events which are of particular 
interest in each phase of combat will be elaborated on in order that the 
reader might acquire an overall appreciation of the contributions each 
submodel makes to the overall model. 

B. SHIP-TO-SHORE PHASE 

1 . Overview 

Since the objectives of the thesis are to provide a user-friendly 
tractible combat model, a number of broad assumptions have been made 
regarding the exact method of employment of the LVA in the ship-to-shore 
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phase of the amphibious operation. First of all, it is envisioned that 
for command and contro purposes, as well as for mine clearing operations, 
there will exist LVA approach lanes as depicted in Figure 3-2, along which 
columns of LVA will transit a 25-mile distance to shore from the Amphibious 
Task Force (ATF). The 25-mile distance is based upon recent requirements 
studies indicating that in future amphibious operations, due to the 
increased lethality of anti -ship missiles and long range artillery, it 
will be necessary to increase the Amphibious Task Force standoff distance 
to approximately 25 miles from shore to reduce the vulnerability of the 
amphibious shipping against this anticipated threat [Ref. 8]. Secondly, 
it is assumed that a maneuver area will exist within which the columns of 
waves of LVA will form into a conventional landing formation composed of 
waves of landing craft as prescribed by current doctrine. 

The two previous assumptions set the stage for the primary assump- 
tion used in computing LVA force level attrition: that is, direct fire 

weapons will be assumed to be the primary anti-LVA threat .-- specifically, 
modified versions of current tank and antitank guided missiles (ATGM) 
assets. Although in reality some attrition of LVA can be expected in 
the maneuver area, it will be assumed that the critical exposure period 
will be that portion of time in the ship-to-shore movement that the first 
assault wave comes within 5,000 meters of the shore defenses until, up to, 
and including the arrival of the last assault wave ashore. Figure 3-3 is 
a flowchart depicting the general sequence of events of the ship-to-shore 
phase model . 

2. LVA Movement Conceptualization 

Two tactical decision variables were utilized for modeling LVA 
ship-to-shore movement: 
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Figure 3-3. Generalized Flowchart for Ship-To-Shore Phase 
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— TBW is the decision variable for the time between successive 



waves. As TBW is shortened, coordination problems will arise resulting in 
confusion on the beach due to insufficient time provided for an assault 
wave to move inland prior to the next wave's arrival. The level of con- 
fusion generated by a short TBW must be balanced against the cost of not 
having sufficiently rapid initial buildup of offensive forces ashore. 

— RD is the distance from the shoreline that each wave will 
commence the transition from planning model to displacement mode. This 
process will be termed a sequential wave transition since each of the 
assault waves sequentially performs the mode transition. This is illus- 
trated in Figure 3-4. The reason for this transition is due to engineer- 
ing stability requirements that this displacement configuration be 
achieved prior to crossing the surf line. The obvious effect of this 
transition is that exposure time to close-in direct-aimed fire will be 
created. 

3. Overall Force Structure 

The model aggregates the combat organizations involved in the 
ship-to-shore phase of the amphibious operation into several homogeneous 
combat units. Each unit is characterized by certain offensive and 
defensive capabilities in comparison to each of the other units. 

Table 3-1 illustrates the combat organizations which have been 
explicitly modeled. The force level of each unit was represented by 
state variables as indicated. The initial force level for each unit is 
input-data to the model. This, therefore, permits the user to investi- 
gate alternative wave composition options as well as various defensive 
scenarios without having to make modifications to the model algorithm. 
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Table 3-1. State Variables Representing Combat Organizations 



Combat Organization 


State Variable 


Shore Defenses - Tank Assets 


DT 


Shore Defenses - ATGM Assets 


DS 


Incoming assault Waves of LVA 
representing waves 1 thru 5 


WV(I) 1=1 thru 5 


A cumulative combat force 
comprised of those Marine 
grounc units which have 
arrived on the beach and 
debarked their LVA 


TLF 


r ire support assets of the 
amphibious task force 

! 


ATFFS 



The tactical combat interactions that exist between these nine combat 
units 'within the overall force structure are illustrated in Figure 3-5. 

4. Shore-Defense Scenario 

The defensive scenario utilized in the model includes a force 
comprised of both tank (DT) and antitank guided missiles (DS). Tank 
ana AT GM units are emplaced 75 meters inland of the waterline at an 
elevation of approximately 5-10 meters. The model does not explicitly 
maneuver or emplace individual tanks or ATOM systems within each unit, 
but aggregates the cumulative effects of the individual vehicles and 
weapons within each category. 

a. Defensive Unit Force Levels 

The state variables DT and DS represent the total unit 
"strengths" in each of the defensive unit categories. Specifically, a 
DT=3 indicates that within the shore defenses there exists a unit of 
tanks having a total combat effectiveness equivalent to three continuously 
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Figure 3-5. Force Interrelationships 
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firing independent weapon systems. A similar interpretation is applicable 
to the state variable DS. 

b. Defensive Fire Allocation 

The two categories of direct-fire weapons are assumed to 
engage targets (incoming LVA) according to a predetermined tactical scheme. 
The defensive "plan" was parameterized as follows: 

(1) Window of Engagement . Each weapon category was assigned 
an engagement window as illustrated in Figure 3-6. Onl> those LVA located 
within the range window could be fired upon by the shore defensive forces. 
The windows are designated by the following input parameters: 





TANK 


ATGM 


Maximum Engagement Range 


TENGMX 


SENGMX 


Minimum Engagement Range 


TENGMN 


SENGMN 



(2) Engagement Rules. Additional defensive tactical criteria 
are implemented into the model logic according to the following rules of 
engagement: 

--A defensive weapon may only engage the two closest incom- 
ing waves if more than two waves of LVA are at any time located within the 
weapon's engagement window. 

--If only one wave of LVA is present in a weapon's engage- 
ment window, defensive fires of that particular weapon type will be dis- 
tributed uniformly against the surviving LVA in that wave. 

--If two waves of LVA are both contained within the engage- 
ment window, defensive fires of that particular weapon type will be 
distributed according to a tactical allocation submodel. A weighting 
factor, (DEFWT) , input by the user is utilized in establishing the 
proportion of the total weapon strength to be allocated against the 
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Figure 3-6. Engagement Window Parameters 



surviving LVA's in each of the two waves. Specifically, if DEFWT(l) = 2, 
and DEFWT (2) = 1, then each surviving LVA in the closer of the two incom- 
ing waves would be allocated twice as much fire as surviving LVA in the 
seaward wave. For example, if waves three and four were both located 
within the ATGM engagement window, then the proportion of DS's fire 
allocated to surviving LVA in wave three would be: 



DEFWT(l) x WV(3) 

x DS 

DEFWT(l) x WV(3) + DEFWT (2) x WV(4) 

where: WV(3) is the state variable for the 

current number of survivors in wave 3 

c. Attrition-Rate Coefficient Computation 

The classical Lanchester hypothesis for aimed-fire attrition 
is that the casualty rate of a unit is proportional to the size of the 
opposing force. If a Unit "A" is being engaged by a Unit "D", this 
action may be expressed by the differential equation: 

dA 

- = x D 

where: Beta^ is called the Lanchester attrition-rate 

coefficient 

It is assumed that this functional relationship holds for 
each pairing (firing unit, target unit) over the small time interval 
dt. The credibility of the model relates the performance characteristic 
data together with the tactical and physical configurations for each of 
the combat units to derive the attrition-rate coefficients. 

It was decided to express the attrition-rate coefficients as 
the product of the rate of fire (ROF) and the single shot kill probability 
(P(k)) as follows: 
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Beta DA = P < k >DA x R0F DA 

where: DA represents a Unit "D" firing on a Unit "A" 

More complicated models exist [Refs. 9 and 10], however, for the purposes 
of the modeling of the ship-to-shore LVA and defender attrition, this 
method was deemed sufficient. 

Attrition-rate coefficients as described above were derived 
for each pairing (defensive weapon, target) yielding the ten variables: 



Beta DT-WV(I) R0F DT-WV(I) X P ^ k ^DT-WV(I) 



I = 1-5 



and 



Beta DS-WV(I) “ R0F DS-WV(I) X P ^DS-WV(I) 1 ~ 1-5 
A switch mechanism is incorporated into the rate of fire (ROF) factor by 
implementing the functional relationship: 



ROF 



D -WV(I) =< 



if WV(I) is located outside 
the engagement window 



if WV(I) is located within 



-j.gP the engagement window 



where: TBF (Time Between Firings) is evaluated by 

TBF = Aim-Reload Time + Pr . 0je clilf Velocity - 

Target Speed 



It should be noted that the relatively slow projectile velocities repre- 
sentative of anticipated ATGM assets in the future does cause such 
velocities to become significant in this computation. 

The second factor used in determining each attrition-rate 
coefficient is the single-shot kill probability (P(k)). It is assumed 
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that a hit by a large caliber projectile would constitute a "kill" in 
that it most likely would inflict damage serious enough either to sink 
the LVA, or render it immobile, thus eliminating it from contributing 
to the buildup ashore. A second assumption is that both defensive 
weapon systems addressed would exhibit normally distributed, unccrrelated 
horizontal and vertical errors. Typical dispersion data, both mean and 
standard deviation, for the Tank and ATGM is required as input-data for 
the hit probability computations. 

The suppressive effects of incoming fire upon each of the 
defensive units was considered a significant factor with respect to its 
effect upon the survivability of the incoming assault waves of LVA. It 
was assumed that the suppressive effect would significantly reduce a 
unit's rate of fire, and also increase the error standard deviation. The 
modeling of these suppression effects was accomplished by assigning a 
relative suppression factor (SUPFAC) in the interval 1, 2, to both the 
Tank and ATGM units. This factor was determined according to the 
following guidelines: 

SUPFAC = 1 No incoming fires (i.e., the defensive 
unit casualty rate is zero) 

SUPFAC = 2 Maximum incoming fires (i.e., the defen- 
sive unit casualty rate is comparable to 
that realized upon full allocation of the 
ATF fire support assets) 

It was assumed that the aim-reload time ^\RTM) would be 
increased by approximately 50 percent under the conditions represented 
by a SUPFAC = 2.0. Within the ROF submodel, this is expressed by the 
1 inear relationship: 

ARTM SUp = ARTM N0Nsup x (0.5 + SUPFAC / 2.0) 
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Additionally, it was assumed that up to 100 percent increase 
in the error standard deviation could be expected under a maximum suppres- 
sion environment, hence: 

ERROR STD SUp = ERROR STD N0Nsup x SUPFAC 
d. Defensive Breakpoint 

It was assumed that if during the course of the ship-to-shore 
movement phase the defensive forces suffered a cumulative loss in excess 
of 70 percent of their initial force strength, the remaining shore defenses 
would withdraw and commencement of the land combat phase of the battle 
would take place. 

5 . LVA Assault Wave Conceptualizatio n 

The model is programmed to handle up to r ive incoming waves of LVA. 
The initial composition of these waves is input by the user by means of 
the variable WVINT. There are no limitations on the number of LVA's 
capable of being in each wave. However, the use 7 ’ is advised that the 
model was intended to model small -unit amphibious operations only. 

0 

a. Wave Posture 

Model functions RNG, HT, and SPD are called upon within the 
model logic to generate the range, height, and speed, respectively, for 
each assault wave as time is incremented throughout the course of the 
ship-to-shore movement phase. The input of tactical employment para- 
meters TBW and RD in conjunction with the physical design parameters 
SPDMAX , SPDMIN, HTMAX , and HTMIN for the LVA being evaluated uniquely 
determines the exact range offshore and vehicle configuration (planning/ 
displacement) for each of the five waves. This information then is imple- 
mented in the rate of fire and hit probability calculations. 
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b. Ground Forces Ashore 

As each assault wave arrives at the beach, the surviving 
strength of that wave is transferred to the variable TLF (Total Landing 
Force Ashore). TLF represents a ground combat force equal to that trans- 
ported by the number of LVA survivors having arrived ashore. Once 
established, the TLF engages the two defensive units allocating its 
fires between the two defensive weapon categories in the same proportion 
as the number of surviving Tanks and ATGM's--that is: 



tlf dt 




DT 


x TLF 


DT 


+ DS 


tlf ds 




DS 


x TLF 


DT 


+ DS 



The casualty rates applied against the DT and DS state that 

survivor variables are determined by means of the Lanchester aimed-fire 

attrition-rate coefficients WBETA ti c nT and WBETA X1 c ,, c by the 

TLF - DT TLF - DS ^ 

equations: 



dDT 

dt 

dDS 

dt 



wbeta tlf _ DT x tlf dt 



-WBETAtlf . DS X TLF 



DS 



The computation of these WBETA coefficients is not performed 
within the model utilizing the detailed rate of fire and P(hit) arguments 
described previously, since in the original LVA assault model developed 
by Chadwick, these parameters were considered to be insignificant in 
relation to the overall model. Chadwick assumed his assault model would 
be used as an auxiliary model to a higher-level model, and would receive 
values for these coefficients from that model. 
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6. ATF Fire Support Conceptualization 

The impact of the amphibious task force's fire support assets 
contribute significantly to the combat effectiveness of the shore 
defense units. Characterizing each of the two defensive force units 
by a simple "located" or "not located" attribute, the attrition rates 
realized by these force units can be simplified substantially by the 
following approach. 

a. "Not Located" Shore Defenses 

At the commencement of the model it is assumed that the 
defensive units DT and DS are emplaced on shore at locations unknown to 
the ATF. The units initially are engaged as "not located" targets by 
area fire for which the following Lanchester area-fire equations are 
appl i cable: 



dDT 



dt 


= (ALPHA dt x 


ATFFS) 


dDS 


A4- 


= (ALPHA ds x 


ATFFS) 



The terms in parentheses on the right hand side of the equa- 
tions are to be considered a generalized input parameter. The combat 
effectiveness of the ATF fire support assets is also to be considered 
relatively constant during this segment of combat time, and thus it is 
possible to synthesize these input factors by examining the attrition 
losses due to area realized in a previous full-scale model calibration run. 
b. "Located" Shore Defenses 

Once a particular defensive unit has initiated its engagement 
of incoming waves of LVA it is considered located. At this point it is 
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assumed that the ATF fire support organization will engage that defensive 
unit through the use of aimed fire. Again it is assumed that the loss 
rate will be in accordance with the Lanchester hypothesis for aimed fire, 
that is: 

dDT 

= Beta nT x ATFFS 

dt Ul 

dDS 

= Beta nc x ATFFS 

dt U5 

It should be noted that the right hand side of both of the 
equations is to be regarded as synthesized factors to be calibrated 
from a previous high-resolution application. 

C. LAND COMBAT PHASE 
1 . Overview 

The land combat phase, like the ship-to-shore phase, has been 
modeled after broad assumptions have been made concerning the type of 
forces modeled and force attrition. These assumptions are quite similar 
in nature to those assumptions made in the ship-to-shore phase of the 
model which would be expected of similar Lanchester-type combat models. 

The first assumption is that of homogeneous forces, which was 
made as a matter of convenience. The defensive forces will be modeled 
as a Tube-Launched, Optical-Guided, Wire Controlled missile (TOW) company 
made up of three TOW platoons located in three separate and fixed defen- 
sive positions. Each TOW platoon will be comprised of three TOW sections, 
and have the capability of withdrawing to an alternate position provided 
as input by the user. The aggressor force ashore will consist of the 
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consolidated surviving landing force ashore which has been redistributed 
into three offensive units. 

The second assumption is that the aggressor force units will fol- 
low three user-defined routes as they advance toward the three defensive 
force positions. Routes can be supplied by the user, or defaulted to 
preassigned routes dictated by the program. However, only three routes 
can be utilized. 

A line of sight (LOS) model written by Professor James Hartman, 
Naval Postgraduate School [Ref. 11], is used, adding great flexibility 
to the modeling of the terrain in the basic scenario. This has a direct 
impact on the probability of detection during any one time period (t,t+dt) 
which is shown as: 



P(Unit i does not detect Unit 
j in a time period t+dt) 

P(Unit i dees not detect 

Unit j in a time per- = x 

iod t+dt) 

P(Unit i does not detect Unit 
j in a time period t,t+dt) 



The first two assumptions provide a basis for applying the third 
assumption which is that attrition of opposing forces will be defined by 
direct-fire Lanchester differential equations similar to the ship-to- 
shore phase, only modeled in more detail. Figure 3-7 provides the scheme 
for the sequence and general flow of events in the model. 

2. LVA Movement Conceptualization 
a. General 

In the original model, aggressor force and defensive force 
locations were modeled in two different ways. Defensive force locations 
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Figure 3-7. Generalized Flowchart for Land Combat 
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were left as user inputs, whereas aggressor force locations had been pre- 
determined by the model builder, and could not be altered by the user. 

This allowed a flexibility of modeling defensive positions, tut flex- 
ibility was limited because of the method of determining movement routes 
for the aggressors. Glenn Mills provided a user option to the model which 
permitted the user to model a variety of aggressor force movement scenarios. 
This option allows for the choice of attack routes and vehicle speed. In 
addition, the option is highly useful to the unfamiliar user of the model 
since unit locations and attack routes can be initially set to the model's 
default values. Different user selected parameters can be input as the 
user acquires a better understanding of the model's algorithm. 

b. Model 

Three predetermined routes are provided for aggressor force 
movements. Each route is subdivided into 40-meter length intervals, 
since a nonfiring aggressor unit is assumed to move one such interval 
during a time-step of 10 seconds (i.e., average speed of 9 mph). A firing 
aggressor unit is delayed a specified number of time-steps before moving 
to the next interval by the state variable NOD. Each interval in each 
route is represented by its center point coordinates, and by its direction. 
If an aggressor unit enters an interval along its associated route, then 
it is considered to be positioned in the center of the interval, generating 
a possible location error of t 40 meters, since this is the distance 
between two consecutive intervals. 

c. User-Defined Routes 

The user is required to input the original location of each 
aggressor unit, and the locations of each of ten nodes he desires the 
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aggressor unit to move through as it advances on the defensive unit's 
position. This information, along with vehicle speed, is used to calcu- 
late route intervals that move the attacking unit through each of the 
designated nodes. A complete route would look like that depicted in 
Figure 3-8. The method used to complete the routes is as follows. 

The straight-line ground distance between the first two adja- 
cent nodes (DIST) is calculated as shown in Figure 3-8. The angle between 
the desired direction of movement and a straight west-to-east movement (a) 
is then calculated. Utilizing these quantities and the distance desired 
to be moved during each time-step (DST), the distance to be moved in the 
X and Y direction (XLN and YLN) is now computed as shown in Figure 3-9. 
These distances are added to the coordinates of the previous interval 
endpoint, point C in Figure 3-9, to determine the coordinates of the next 
interval endpoint, point D. This same distance is again added to compute 
the coordinates of the next endpoint. Point E. This process is continued 
until the distance from the last endpoint computed to the next node is 
less than DST. This general process is repeated for each pair of nodes 
until the entire route is completed, or the unit’s force level is reduced 
to zero or the battle terminated, whichever comes first. 

To insure that all intervals are of equal length, the computa- 
tion of the first interval between any two nodes must be considered separ- 
ately by taking into account the distance left over from the last 
computation between the previous two nodes. To accomplish this, the 
first interval takes the remaining distance (e) and adds it to an 
interval length (DST - e) for the first interval between any two nodes. 
This insures that each interval along the route is of length DST, which 
is the required length. 
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Figure 3-8. User Determined Routes 
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Figure 3-9. Route Computation 
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3. LOS, Detection, and Fire Allocation 



a. LOS 

The existence of a line-of-sight between any two opposing 
units is determined utilizing a line-of-sight model written and program- 
med by Professor James K. Hartman, Naval Postgraduate School [Ref. 12], 
and is listed as Subroutine LOS in the land combat phase of the model. 
Professor Hartman's model utilizes a parametric terrain model proposed 
by Needles [Ref. 13], which represented terrain by modeling individual 
hill masses. Each hill is described by a bivariate normal density function, 
and fitted together to form a section of terrain utilizing the following 
information illustrated in Figure 3-10: 



1) (Xc,Yc) - Coordinates of each hill's centerpoint 

2) PEAK - Peak height of each hill 

3) a - Standard deviation corresponding to 

x the X-axis 



4) 



y 



- Standard deviation corresponding to 
the Y-axis 



5) (p) - Rotation factor 



Once the terrain has been "mapped", the existence of a line- 
of-sight can be determined for each pair of opposing units. The informa- 
tion required to accomplish this is the location and elevation of each 
unit, as well as the height of the vehicle each unit uses. Professor 
Hartman's model yields the fraction of aggressing Unit A as seen by 
defending Unit B, and the fraction of defending Unit B as seen by 
aggressing Unit A. Figure 3-11 is used to illustrate this. 
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Figure 3-10. Terrain Conceptualization 
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Figure 3-11. Partial LOS Conceptualization 
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b. Acquisition 



The acquisition process was well-modeled in the original land 
combat model devised by Joseph Smoler. The model employs the concept of 
parallel acquisition, whereby the weapon system continuously searches for 
targets, even while engaging other targets. When such a weapon system 
kills its presently engaged target, it immediately can shift its fire to 
a new target, provided that such a target has been acquired either during 
the engagement of the previous target just killed, or earlier [Ref. 14]. 

A general description of the manner in which Smoler modeled target-acquisi- 
tion is provided here. However, a more detailed description is provided 
in his thesis. 

The probability that a Unit j is detected by a Unit i at 
time t+dt was modeled for four different combat situations in which the 
opposing forces might find themselves. These situations can be summarized 
as follows: 



Observer 

Not firing (t,t+dt) 
Not firing (t,t+dt) 
Firing (t,t+dt) 
Firing (t,t+dt) 



Target 

Not firing (t,t+dt) 
Firing (t,t+dt) 
Not firing (t,t+dt) 
Firing (t,t+dt) 



The formulas derived to compute the probability of detection 
for each of these situations have a number of common variables, therefore 
their definitions are provided beforehand for clarity: 

P . . (t=dt) = The probability that a typical firer in Unit 
i has acquired one or more targets of type 
j by time t+dt 
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QijU+dt) = [l - P i j (t+dt)D 



QV ij .(t+dt) 



e - x ij dt V t} 





The probability that target j is not 
visually detected by Unit i during 
(t,t+dt) provided Unit j does not 
fire during this time interval 


where: 


s j (t > ■ 


the number of survivors 
in Unit j at time t 


ar.d: 


X ij 


the nonfiring detection rate 
of one target in Unit j by 
one observer in Unit j 



QP ij (t+dt) = (1 - P k ) FR j dtS j (t) 

The probability that target j is not 
detected by a launch signature during 
(t,t+dt) provided that target j fires 
during this time interval 

where: P. = The probability that one 

observer in Unit i is look- 
ing in a direction which 
enables him to detect target j 

and: FR . = The firing rate of one 

J firer in Unit j 



The first situation occurs when neither the observer nor the 
target is firing during the interval (t,t+dt). This situation allows the 
observer to conduct search operations only, thereby maximizing the prob- 
ability of detecting a target in his sector of responsibility, and has 
the target maximizing his probability of not being detected by exposure 
to an observer by a launch signature. Thus, the probability of not 
detecting in time interval (t,t+dt) is 

Q.j j (t ,t+dt) = Q-jjCt) x QV i j ( t , dt ) 
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The second situation occurs when the target is firing during 
the search interval (t,t+dt) while the observer is conducting only search 
operations. This provides the observer with additional information to 
assist in detection of the target. The observer will detect the target 
by the target's launch signature. Thus, the probability of not detect- 
ing in time interval (t,t+dt) is 

Q-j j (t,t+dt) = Q..{t) x (QV.jU.ot) + QP i j ( t ,dt) - 
QV^t.dt) x QP ij (t,dt)) 

The third situation occurs when the observer is firing during 
the search interval (t,t+dt) while the target is maximizing the probabil- 
ity of not being detected by not firing during the interval. The 
observer has lowered detection probability by diverting a portion of 
his force to firing on a known target. A new factor is introduced 
which will alter the probability of detection, namely the event: 



A = The situation in which Unit j is within the 
field of view of Unit 1 , with at least one 
of the targets at which Unit i is firing 



This states that Unit j, which is not currently firing, happens to expose 
itself to firing Unit i when firing Unit i is looking and firing on at 
least one other target in j's principal direction. Thus, the probability 
of not detecting in time interval (t,t+dt) is 



0 



if event A occurs 



Q n -j(t,dt) 



g(n) if j is an_aggressor unit 
= < and event A occurs 

Q..(t) if j is a ^defending unit 
a and event A occurs 
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where: 



and: 



g(n) is an increasing function of n, 
where n is the number of time intervals 
elapsed since time t. 

9(0) = Q.j j(t) 

Q - - ( t: ) 4 g(n) 4 1.0 for all n 

* J 



The fourth situation occurs when both the observer and target 
are firing during the interval (t,t+dt). In this situation the observer 
has minimized his searching capability, and the target has maximized its 
probability of being detected. Thus, the probability of not detecting 
in time interval (t,t+dt) is 



Q i j(t,t+dt) =< 



Q-jj(t) x QV ij -(t,dt) 


if event A occurs 


g(n) 


if j is an aggressor 
unit and event A occurs 


Q lM (t) 


if j is a defender unit 


T J 


and event A occurs 



where: 



QV*j(t,dt) 



x u dts ; (t) 



and: XT- = x RF 

ij ij 

RF = Reduction Factor (the detection rate 
of Unit i has to be reduced since this 
unit fires during (t,t+dt) and the 
search for targets is not as effective 
as for a nonfiring unit) 

S^(t) - Sj(t) X ( K r k PTT iK ) 



PTT-,, = proportion of Unit i allocated to 
1 Unit K 



k = (Unit K is engaged by Unit i and Unit j 
is within the field of view of Unit i 
while observing Unit K) 
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If a 1 i ne-of-s ight does not exist between observer i and 

target j, then no accumulation of detection probability will take place 

during the current time interval (i.e., P..(t) will remain the same). 

^ J 

However, if a 1 ine-of-sight does not exist throughout more than three 

consecutive time intervals, then the P. . is set to zero (i.e., P..(t) = 

* J "l J 

0) and the accumulation process will start again from zero if a line-of- 
sight is acquired at a later point in the battle. The motivation for 
this decision rule is seen by the observation that even if observer i 
loses a 1 ine-of-sight with target j for a short period of time, he still 
probably has some idea of where to expect the target to reappear, 
c. Non-Firing Detection Rate 

The situations that occurred when the target was in a non- 
firing status had detection probability functions that had as a para- 
meter A.j, a non-firing detection rate. The manner in which the model 
derives this rate is quite detailed, and deserves attention. 

To begin, each firer in an observing unit is assigned a 
search section (or sector of responsibility) which is characterized by 
two parameters (see Figure 3-12). These parameters are the section 
width (ISECWD), and the primary direction of search (IPRDIR). Further- 
more, it is assumed that the search direction within a section of search 
has the following probability density function known as the LIMICON 
Function: 



f (9) = A + B cos 9 -D < 9 < D 



where: D = 



ISECWD/2 




Figure 3-12. Search Direction 




Figure 3-13. Observer-Target Scheme 

(A = Observer, 3 = Target) 
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9=0 corresponds to the observer's 
primary direction of search 

note: A and B are chosen such that 

/ f( 0 ) dQ = 1 

“D 



:o determine the probability that observer A is looking in 
a direction which enables him to detect target B, P k is the value of the 
LIMICON function integrated from an angular value up to 15° on either side 
of the primary direction of fire, specifically: 



-ANGLFT 

= / f(9)d9 = shaded area. Figure 3-13 

K ANGRT 



where : 



ANGLFT = < 



ANGLE + 15° 
D 



if ANGLE + 15° < D 
if ANGLE + 15° > D 



and: ANGLE = the absolute value of the angle between the 

the primary direction (IPRDIR) and the ob- 
server-target direction (OTANG) 

ANGRT = ANGLE - 15° 



Now, given that observer A is looking in a direction a such 
that ANGRT ~ a = ANGLFT, the conditional detection rate (X^l ANGRT 4 a 
4ANGLFT) is determined by a regression curve [Ref. 15] which is a 
function of the terrain, target horizontal velocity, and the equivalent 
range for a full height target. This detection rate of one observer 
detecting one target becomes 



Aab = ^ x ab I ANGRT = a = ANGLFT) x P k 
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d. Fire Allocation 

Three conditions are necessary for Unit j to be classified 
as a target for Unit i. First, a 1 ine-of-sight must exist between Unit 
i and Unit j. Second, the range between the units must be within the 
maximum range of Unit i's weapon system. Lastly, the probability that 
a detection of Unit j is made by an observer in Unit i in the time period 
t+dt must be greater than 0.C0. 

Once these conditions are satisfied, the manner in which fire 
is allocated to a target depends upon how many targets are to share in 
the firepower of Unit i, and what distance exists between i and the 
new target in relation to the other targets under fire. The priority of 
fire naturally will go to the closest target since it is of a greater 
threat to Unit i than the more distant targets. The amount of firepower 
available from Unit i is naturally a function of the percentage of 
surviving force available to fire in Unit i. 

4. Attrition 

Attrition of forces is assessed based upon variable coefficient 
Lanchester equations of modern warfare [Ref. 16]. This method of attri- 
tion assessment was used by David Chadwick in the ship-to-shore phase 
of the model, however, in less detail than was modeled by Joseph Smoler 
in the land combat phase of the model. The "extra" detail provided by 
Joseph Smoler is the generation of the conditional probability of a kill 
given a hit. This probability was stated by Chadwick as a user-supplied 
input parameter. 

The restriction of the model to aimed-fire weapons systems and 
homogeneous forces allows for the attrition of forces to be assessed 
using variable coefficient Lanchester equations of modern warfare. 
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The attrition for a defending Unit j is described by the following 
differential equation: 



dSj(t) 

dt 



(A-j j x PROP i j ) x S.(t) 



where: S^(t) 



The number of survivors in Unit k at time t 



Ajj = The rate at which one fire*' of Unit i kills 
Unit j targets (attrition rate of Unit j by 
one firer of Unit i) 

PROP^i = Proportion of Unit k allocated to fire 
against a Unit 1 



These basic differential equations of force-on-force attrition were 
approximated by the following Euler-Cauchy difference equations: 



S . (t+dt) = Max (0,S - (t) - ZA..(S.(t) x PROP.Jdt 

* J * J J * 

for each defending Unit i 

and: 

S . (t+dt) = Max (0,s • (t) - zA..(S,(t) x PROP-Jdt 

J * J I * J 

for each aggressor Unit j 



The manner in which the attrition-rate coefficient A. • is 

J 

derived stochastically already has been discussed in the model 
enhancement chapter, therefore, only a description of how the deter- 
ministic attrition-rate coefficient is derived will be mentioned here. 

The attrition-rate coefficient. A.., for each equation is 
computed according to the equation: 




E r T iP 
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where: T.. = the time for one firer of Unit i to kill 

one target of Unit j under the conditions 
in the present time interval 



T.j j is computed using the Bonder-Farrel 1 formula [Ref. 17]: 






t + t, + + 1 - p ( h lh) 

3 P(k|h) P(k|h) 



+ P (h I h) - P(h) 



where: t 

d 



t 



l 



p(h) 

P(h|h) 



P(h|m) 

P(k|h) 



Time to acquire a target 

Time to fire first round after a 
target is acquired 

Time to fire a round following a hit 

Time to fire a round following a miss 

Projectile's time of flight 

Probability of a hit on first round 

Probability of a hit on a round given 
that the prior round fired was a hit 

Probability of a hit on a round given 
that the prior round fired was a miss 

Probability of a hit on a round given 
that the round fired was a hit 



There are two assumptions of the Bonder-Farrell formula that are 
implied by the model. The first assumption is that fire is Markov- 
Dependent in that the probability of a hit of any round depends only 
upon the result of the previous round. The second assumption is that a 
Geometric Distribution describes the parameter P(k|h) in that accumulated 
damage is considered to be negligible. 

The expected value of T • • , E [ T . - ] , may now be expressed for 

I J * J 

each weapon system in the model. It is assumed that for the TOW 
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weapon system P(k h) = 1.0, and P(h m) = P(h h) = P(h), which results 
in the reduced formula: 



E[T ij ] = l a + *1 + tf + 



(t m + V x (1 - p ( h ^ 
P(h) 



If the firing weapon system is a tank, then it is assumed that P(k|h) = 1.0 
(due to a lack of information), and that t f = 0.0 (due to the velocity 
of the projectile). Thus, in this case the formula becomes: 

t 

E[T j;j ] - t a + t, + x (1 - P(h)) 

P(h|m) 



It should be noted that all targets were considered to be stationary 
throughout the attrition process. This is obvious in the case of the 
stationary defending forces, and was assumed to be the case for the 
aggressor forces due to the fact that the hit probability of a TOW 
against a moving target is almost the same as for a stationary target. 

5. Battle Termination 

Two criteria were used as rules governing battle termination. 
The first criterion was the annihilation (zero force level) of one of 
the two forces. The second criterion was that the distance between 
defender and aggressor forces is too small. 

The first criterion is an intuitively obvious reason for ter- 
minating the battle, and thus easy to model. However, although the 
reasons for the second criterion might be as obvious, the modeling of 
this is not simple. The manner in which Glenn Mills modeled it was to 
compute the distance between each attacking sub-unit on which casual- 
ties were being assessed (i.e., still alive), and each defending 
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sub-unit that was still alive. If any one of these distances between 
active sub-units was too close, the battle was considered to have 
reduced to close-in, hand-to-hand combat. The outcome of this type 
of combat is not currently provided for in the model, and for this reason 
the battle is simply terminated at this point. However, to insure that 
the aggressor units do not pass by the flanks of the remaining defending 
forces ard remain outside termination distance, a check is made of the 
location coordinates of each sub-unit. If any aggressing sub-unit's X 
coordinate places the unit beyond the location of the most forward 
defending sub-unit still in the battle, the battle also is terminated. 

The specification of the distance between forces for battle 
termination is left as a user-input, which provides added flexibility 
of breakpoint distance analysis. In particular, it lends itself to the 
study of optimum breakpoint distances for various weapons systems on 
the battlefield. 
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IV. FUTURE ENHANCEMENTS 



A small -unit amphibious operation combat model has been presented 
in this thesis which emphasizes the simplistic and avoids the abstract 
to provide an understandable and, more importantly, a useable combat 
model for students of con, bat modeling. However, the combat model pre- 
sented has the potential of being developed into a much more refined 
model which could be studied and utilized by more experienced combat 
modelers. Therefore, several enhancements which might be of some benefit 
to the more experienced modeler are mentioned here as possible approaches 
that could be taken in refining the present model. 

A. HETEROGENEOUS FORCES IN THE LAND COMBAT PHASE 

The current land combat phase of the model involves combat between 
homogeneous forces only— that is, each force is comprised of only one 
weapon system type. This type of force structure was intentionally 
modeled to maintain a relatively simple model to understand and work 
with. However, added flexibility could be attained by modeling multiple 
weapons system types for each of the opposing forces. This would allow 
the user to analyze the effect that different force mixes would have on 
battle outcome. 

The addition of different weapons system types within a single unit 
would require extensive restructuring of the attrition process currently 
used in this model. Although Lanchester equations still could be uti- 
lized in computing direct-fire weapon system attrition, separate Lan- 
chester equations would have to be provided for each weapon system. 
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Furthermore, with the addition of indirect-fire weapon systems (i.e., 
artillery, naval gunfire, and close air support) Lanchester equations 
for area-fire would have to be implemented fo»~ each area-fire weapon 
system type. The total attrition of any particular unit then would 
be the summation of the damage assessed by each weapons system type on 
the target being attrited. 

An enhancement of this type would result in more realism at the 
cost of longer execution time, and a more complicated attrition process. 
Since the original intent of the thesis was to provide a simple model 
to understand, it would be advisable to retain a copy of the original 
model prior to adding this enhancement. Then a simple model would 
still be available to the less experienced combat modeling students, 
while a more detailed model would provide the realism that more exper- 
ienced modelers would demand. 

B. LOGISTICAL SUPPORT 

Logistical support is one of the most overlooked factors of combat 
in the development of combat models. The influence that the resupply 
of fuel and ammunition alone have on the outcome of a battle is obvious 
and deserves attention. 

Ammunition and fuel consumption could be modeled along the same 
lines as attrition (i.e., through the use of expected values of con- 
sumption). When ammunition or fuel on hand reaches a specified critical 
level, a unit could be restricted in movement, or experience a reduced 
level of fighting effectiveness and maneuverability (based on a shortage 
of ammunition and fuel) until resupply of the critical resource could be 
obtained. 
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The amount expended of these resources would necessarily be a function 
of the number of surviving firers in the unit, the number of vehicles 
available to transport the unit, and the number of targets engager, by 
the unit at any one time interval. The expected values of these items 
then could be used in computing the expected rate of consumption of 
ammunition and fuel. Therefore, the overall process could be modeled by 
initially allocating specific levels of these resources (i.e., ammuni- 
tion and fuel) to each unit at the commencement of the battle, and 
subtracting the expected expenditure of the ammunition and fuel of a 
particular unit based upon the expected number of survivors firing on 
engaged targets, and the distance traveled by the expected number of 
surviving vehicles of the unit. 

C. GRAPHICAL BATTLE SUMMARY 

A graphical display of what is taking place on the battlefield can 
be worth a thousand words to the user of a combat model. Plotting 
unit locations and force levels on a display of the actual terrsir; 
fought upon would eliminate time-consuming interpretation of these 
results from a printed battle summary report. An enhancement ot this 
sort would serve both the experienced and inexperienced users ot the 
model. The inexperienced user would have results displayed in a much 
more understandable format, while the experienced user would be aole 
to study variant combat scenarios with much less effort and time 
expended. 



68 



V. FINAL REMARKS 



The purpose of the model that has been developed is to illustrate 
a number of underlying concepts of combat modeling which have been 
addressed in this study. Therefore, it seems appropriate to readdress 
these concepts to allow the reader to reflect upon them in light of 
what has just been presented. 

A. INTEGRATING INITIALLY INDEPENDENT COMBAT MODELS 

The model developed here was made up of two sub-models: ship-to- 

shore and land combat models. These sub-models, as discussed earlier, 
utilized similar combat modeling methodology (i.e., Lanchester equa- 
tions) in computing force level attrition. However, each sub-model 
was developed by different individuals, which created several problems 
when the two separate sub-models were integrated into a singular con- 
tinuous flow algorithm. In particular, individualized FORTRAN coding 
techniques and documentation of state variables within the program 
structure required the restructuring of major portions of FORTRAN 
code to make the overall combat model tractible and understandable. 
This serves to illustrate the need for a standardized programming 
technique to be applied to programming of combat models, and high- 
lights the need for proper planning and coordination in development 
of large scale combat models by teams of combat modelers. 

B. THE USER-ORIENTED APPROACH TO COMBAT MODELING 

This thesis illustrates the desirability of a user-friendly 
approach to combat modeling. It was a major contention of the thesis 
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that this approach to combat modeling has not boon closely addressed 
by combat modelers providing combat models for the United States mili- 
tary. Furthermore, it was mentioned that the lack of concern given 
to this approach of combat modeling might help to explain the lack of 
enthusiasm exhibited by the United States military in utilizing combat 
models for the training of field commanders and staffs. The thesis had 
as one of its purposes, the presentation of a combat model designed to 
be easily understood and utilized by intended users, combat modeling 
students. Combat models should be designed and documented with the 
user's capabilities and needs in mind, as opposed to those of the 
programmer. 

C. A COMBAT MODEL FOR STUDENT USE 

The small-unit amphibious operation combat model presented here is 
a basic Lanchester-type combat model which has been designed with a low 
level of complexity in order that it might be understood more easily, 
and studied by students of combat modeling. It has been recognized 
that combat modeling students may have little or no experience of the 
governing theory, and therefore would comprehend the theory of combat 
modeling more easily by utilizing and understanding its basic application. 
For this reason, enhancements that would increase the complexity of the 
model are discouraged, and enhancements that would make the model more 
understandable (i.e., graphical battle summary reports) are strongly 
encouraged. 
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I. INTRODUCTION 



The purpose of this manual is to familiarize the user with the model, 
and to provide administrative information describing how the potential 
user would access and run the model. 

The small-unit amphibious operation combat model is a two-phased 
combat model which conducts both ship-to-soore and land combat. The 
model uses both aimed and area-fire Lanchester-type equations for 
casualty assessment. The battle is initiated by an amphibious task 
force positioned 25 miles offshore from an opposing defensive force 
which is illustrated in Figure A-l . If an amphibious landing is success- 
ful, land combat will be conducted inland over a 10 x 10 km piece of 
terrain representing an area east of Fulda, West Germany, known as the 
Fulda Gap, which is depicted in Figure A-2. 
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Figure A-1 . LVA Approach Conceptualization 
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Figure A-2. Land Combat Terrain Model 
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II. AVAILABLE OPTIONS 



The model has been developed with a number of options available to 
the user to provide more model flexibility for the more experienced 
user. Each of these options, including user responsibilities, is 
discussed here with the input requirements for each being outlined in 
the next section. 

A. STOCHASTIC VS. DETERMINISTIC ATTRITION 

The user has the option of using stochastic or deterministic attri- 
tion computation. Both methods utilize Lanchester aimed-fire equations; 
the difference between the two is the method of calculating the attrition- 
rate coefficients used in the Lanchester equations. 

Deterministic attrition can be thought of as the expected value of 
attrition, and is implemented by using the Bonder-Farrel 1 method of 
calculating the attrition-rate coefficient. A... The stochastic method 

* J 

can be thought of as the randomization of attrition, and is implemented 
by using random deviates from a Beta Distribution in conjunction with 
the range of a target to generate individual attrition-rate coeffi- 
cients for each unit at each time-step. 

B. VARIANT ATTACK ROUTES 

The user has the option of providing variant aggressor force attack 
routes. The user can utilize the program's straight west-to-east routes, 
or can input desired altered routes for aggressor force units to follow. 

To select new routes, a user must input the number of nodes desired on 
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each of three routes, and the coordinates of .each of these nodes. The 
program then will compute routes through each node. The nodes must be 
inputed in order from west to east, and should not create an angle 
between the west-to-east axis and the route direction that exceeds 45°. 

C. ALTERNATE DEFENSIVE POSITIONS 

The user has the option of implementing alternative defensive unit 
locations. This option permits the user to add more realism to the 
model by allowing the defending units to withdraw to alternate positions 
when their primary positions become untenable (i.e., distance between 
opposing forces is too close). This breakpoint distance is determined 
and inputed by the user, and also is used as the distance for battle 
termination in the event that the battle reduces to close-in combat (i.e., 
hand-to-hand). The alternative to moving the defenders is to terminate 
the battle when the breakpoint distance is initially reached. 

D. BATTLE StMMARY PRINT-OUT 

The user has the option of limiting the printed output of the model. 
The user can receive a battle summary print-out at the completion of 
each 10-second time interval, or this information can be suppressed, 
printing out the results only after each phase of combat. 



76 



III. REQUIRED INPUT 



The small-unit amphibious operation combat model presented in this 
thesis has been provided with a blank data set (see Appendix D) which 
includes each variable of the model requiring input provided by the 
user, and space available following each variable for the user to place 
the desired variable value. However, the definition of each input 
variable may not be familiar to the first-time user of the program. 
Therefore, the following list of input variables and their definitions 
is provided as a quick reference for the user of the model. 



Input Variable 
IPRINT 



SPDMAX 

SPDMIN 

HTMAX 

HTMIN 

WIDTH 

TENGMX 

SENGMX 

SENGMN 

TARTM 

SARTM 



Sh ip-to-Shore Phase 

Definition 

User option for selecting type of battle 
summary report desired: 

0 - Each Time-Step 

1 - End of Battle 

Maximum speed of LVA in the water. 

Minimum speed of LVA in the water. 

Height of LVA above water at maximum speed 
Height of LVA above water at minimum speed 
Width of an LVA. 

Tank maximum engagement range. 

ATGM maximum engagement range. 

ATGM minimum engagement range. 

Tank aim-reload time. 

ATGM aim-reload time. 
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Input Variable Definition 



TVEL 


Tank projectile velocity. 


SVEL 


ATGM projectile velocity. 


TSIGV 


Standard deviation error in the vertical 
axis for Tank fire. 


TSIGH 


Standard deviation error in the horizontal 
axi s for Tank fi re. 


TMEANH 


Bias error in the horizontal axis for Tank 
fire. 


SSIGV 


Standard deviation error in the vertical 
axis for ATGM fire. 


DEFWTS 


Defensive force tactical allocation weights 


WVINT(i) 


Initial strength of assault wave I. 


DINIT (i ) 


Initial strength of defensive Tank (1=1) 
and ATGM (1=2) units. 


A(i) 


Aggressor force attrition coefficients. 


BO) 


Defensive force attrition coefficients. 


WB(i) 


Aimed-fire attrition-rate coefficients 
for defensive force Tank and ATGM units. 


GAINL 


Defensive force attrition level at which 
remaining defending forces withdraw and 
ground assault commences. 


GAMMA 


Aim-reload time suppression factor. 


DELTA 


Aiming error caused by the suppression 
factor of ATFFS 



The remaining portion of the input data refers to the terrain model 
developed by Professor James Hartman. It is suggested that this portion 
of the data set not be altered until the user has studied and fully under 
stands the Hartman terrain model. 
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Land Combat Phase 



Variable Definition 



ITRIT* 


Input variable denoting whether attrition 
will be stochastic or deterministic: 

0 - Stochastic 

1 - Deterministic 


DSEED** 


Double precision seed used in the Beta 
Distribution Random Deviate Generator. 


PP - QQ 
PD - QD 


Input parameters for the Beta Distribu- 
tion Random Deviate Generator: 

PP-QQ Aggressor force 
PD-QD Defensive force 


NBU 


Number of defensive units. 


NRU 


Number of aggressor units. 


RMINTK 


Minimum effective range of an LVA weapon 
system. 


RMAXTK 


Maximum effective range of an LVA 'weapon 
system. 


RMINTW 


Minimum effective range of a defensive 
TOW weapon system. 


RMAXTW 


Maximum effective range of a defensive 
TOW weapon system. 


IRTE 


User option for selecting type of aggressor 
force attack routes: 

0 - Program determined 

1 - User determined 


ISPD 


Speed of aggressor force units: 
1 - 9 MPH 

2-12 11 

3- 15 " 

4- 18 11 



i L 

Xic(i.-j), YIC(I,j) Coordinates of the j tn interval endpoint 



N(i) 


of the route for Unit i. 

Number of nodes for aggressor route i. 



Note: *There are two ITRIT variables in the data set. The first ITRIT 

refers to the aggressor forces. 



**There are two DSEED variables in the data set. 
refers to the aggressor forces. 
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The first DSEED 



Variable 



Definition 



XLOC ( i , j ) , YLOC ( i ,-J ) 


Coordinates of node i for aggressor route i 


X ( i ) , V ( 3 ) 


Location of defensive Unit i. 


FL(3 ) 


Force level of a defensive Unit i. 


IPRDIR(i) 


Principal direction of fire of defensive 
Unit i . 


IALT 


User option for selecting alternate 
defensive positions: 

0 - Yes 

1 - No 


BREAK 


Breakpoint distance between aggressor 
units and defensive units. 


ITEM 


Input variable denoting number of time- 
steps allowed for aggressor unit moves. 


XA(i) ,YA(i) 


Coordinates of alternate position for 
defensive Unit i. 


P(i J) 


Probability of first round hit by Unit i 
in range band j. 


PHH ( 5 , J ) 


Probability of a hit following a hit by 
Unit i in range band j. 


PHM(i,j) 


Probability of a hit following a miss by 
Unit i in range band j. 


PKH ( i , j) 


Probability of a kill given a hit by 
Unit i in range band j. 
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IV. EXPECTED OUTPUT 



The small-unit amphibious operation combat model's output is 
designed to be self-explanatory. Each phase of the amphibious opera- 
tion is reported in the output of the model. The output format for 
each phase will include an initial information section to provide the 
user with feedback concerning the operation of the model as read-in 
by the model from the user-supplied input data. This serves as a check 
and a record for the user to insure that the model was run according to 
the intended design of the user. Secondly, battle summary reports are 
provided at specific points of the battle depending upon the desires 
of the user as input by the user option variable IPRINT. An example 
of the model's output is displayed in Appendix F. 
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V. ACCESSING AND EXECUTING THE MODEL 



The prospective user who wishes to study the small-unit amphibious 
operation combat model must first contact Professor James Taylor of the 
Operations Research Department and obtain the user identification 
number and password for the disk space containing the model and its 
support programs. 

A. ACCESSING THE MODEL 

Once the required information is obtained, the user should proceed 
to LOG ON to his OWN disk space entering the CMS mode of operation. 

Upon entering CMS, the following commands should be executed: 

LINK TO (USER ID*) 191 AS 192 RR 

PASSWORD 

ACCESS 192 B/A 

COPY FILE AMPHIB FORTRAN B = = A 
COPYFILE SEA DATA B = = A 
COPYFILE LAND DATA B = = A 
COPYFILE BSEA DATA B = = A 
COPYFILE BLAND DATA B = = A 
COPYFILE WAR EXEC B = = A 
RELEASE 192 (DET 



*Note: USER ID refers to the user id provided by Professor Taylor. 
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What is received on the user's disk ib a copy of the following files: 

1. The Small-Unit Amphibious Operation Combat Model (APPENDIX B). 

2. A complete data set: SEA and LAND (APPENDIX C). 

3. A blank data set: BSEA and BLAND (APPENDIX D). 

4. The model's executive program: WAP (APPENDIX E). 

B. EXECUTING THE MODEL 

To execute the model utilizing the data set provided, the user must 
first compile the FORTRAN program, AMPHIB, by entering the following 
commands: 

DEF STOR 1M 

I CMS 

FORTGI AMPHIB 

Once the program is compiled, the user enters the name of the 
executive file WAR, which then executes the program and displays the 
listing file of output from the model, (i.e., AMPHIB1 LISTING (APPENDIX 
F) ) in the BROWSE mode of XEDIT. 

C. ALTERING THE DATA SET 

The user may desire to invoke one of the available options provided, 
or alter specific elements of the existing data set to "play out" various 
combat scenarios. To alter the existing data set, the user first decides 
whether to alter the ship-to-shore phase of combat, or the land combat 
phase. Once this has been established, the user can simply XEDIT the 
appropriate data file, replacing the old input data with the new input 
data. 

To construct an entirely new data set, the user should make use of 
the blank formatted data set provided. The user simply XEDIT ' s the 
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BSEA or BLAND data files, inputting new data by typing over the spaces 
provided. The variable names are listed in both of the data sets, as 
well as in Chapter III of this user's manual. The space provided in 
the blank data sets is designed to be compatible with the READ format 
statements of the program. 

D. EXECUTING THE MODEL AFTER ALTERING DATA 

If the user has just altered specific elements of the data set pro- 
vided without altering file names, the user will once again enter the 
name of the executive file WAR, and enter the new data set file names 
where appropriate. Once this editing of the executive file has been 
accomplished, the user simply enters the executive file name WAR to 
execute the model again. 
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VI. PROGRAM STRUCTURE 



The small-unit amphibious operation combat model is a computerized 
model written in FORTRAN. It consists of a main program and 19 sub- 
routines. To assist the user in understanding the operation of the 
model, a brief description of the function of each subroutine, as well 
as the functioning of the main program, is provided. 

A. MAIN PROGRAM 

The main program serves as a director program for the model. It call 
for the initialization of data for the ship-to-shore phase of combat, and 
then commences the execution of that phase of combat. The results of 
the ship-to-shore phase of combat as provided by subroutine SEA are 
then reviewed to determine if the land combat phase of combat should 
begin, or if the battle should be terminated. If the results warrant 
a continuation of the battle, the reason for continuation is printed 
and land combat is initiated. 

B. SUBROUTINES 

There are 19 subroutines in the model. The function of each has 
been provided at the beginning of each subroutine in the coded program, 
and also is presented here for clarity. 

1. Subroutine SEA 

This subroutine is the main driver program for the ship-to-shore 
phase of the amphibious operation. Its main purpose is to initialize 
key parameters, and to direct program flow in the ship-to-shore phase 



of combat. 
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2. Subroutine RKINT 

This subroutine provides the interface between the EULER numerical 
integration routine (RKLDEG) and the subroutine ATTR which determines 
each unit's status as time progresses throughout the amphibious operation. 

3. Subroutine ATTR 

This subroutine determines the attrition rates and updates the 
status of each unit with respect to shore movement based upon the given 
state variable strengths, and implements this information into the attri- 
tion loss-rate computation. 

4. Subroutine DTGTS 

This subroutine determines the wave numbers that are to be engaged 
by the defensive Tank and ATGM units, based upon the engagement window 
criteria and LVA wave survivor force levels. 

5. Subroutine DATAIN 

This subroutine reads in all user-supplied information required 
by the ship-to-shore phase of the model. 

6. Subroutine OUTPUT 

This subroutine provides an input summary printout based upon 
the data received by subroutine DATAIN. A printout of dispersion data 
generated as a result of data supplied also is provided. 

7. Subroutine PHIT 

This subroutine computes the probability of a hit based upon 
the range, width, and height of a given target. The type of weapon 
being employed against the target then is taken into consideration for 
computing the specific probability of a hit. 
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8. Subroutine INTRP 



This subroutine is a check to insure that the range of a target 
and the dispersion data are compatible for computing the probability of 
a hit in subroutine PHIT. 

9. Subroutine RATE 

Given the range and speed of a target, along with the type of 
weapon being used to fire upon the target, and the suppression factor 
the firer is being subjected to, subroutine RATE computes the rate of 
fire used against a particular target. 

10. This is the primary subroutine of the land combat phase of the 
amphibious operation. Information required for the operation of the 
land combat phase is read-in and printed in a summary table for user 
review. The information provided by all other subroutines used in the 
land combat phase is used in this subroutine as input to the basic land 
combat algorithm. 

11. Subroutine SETUP 

This subroutine is used to read-in the terrain data and create 
parametric terrain. The terrain data will be used when computing 
line-of-sight between targets and observers, as well as providing a 
grid system for unit locations and movement. 

12. Subroutine ROUTE 

This subroutine computes the route of each aggressor unit when 
the user has selected the option of inputting aggressor routes. It 
calculates the coordinates of each interval endpoint along the route, 
making each interval length (distance moved during a ten-second time- 
step) the same. The interval length is determined by the speed the user 
has selected and inputted for the current battle. 
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13. Subroutine LAMBDA 



This subroutine used in conjunction with the 1 ine-of-sight 
routine computes the detection rate (DETRAC) of target j by the observer 
i, given the percent of target visible (PCTVIS) to the observer. 

14. Subroutine ELEV 

This subroutine determines the terrain elevation for a given 
set of X, Y coordinates. This function is used in conjunction .vi th the 
1 ine-of-sight subroutine in computing a 1 ine-of-sight between observer 
and target. 

15. Subroutine STOCH 

This subroutine determines the attrition coefficients when a 
user has selected a stochastic attrition option. The calculation is 
a function of the original stochastically determined attrition coeffi- 
cient, as well as a function of range. 

16. Subroutine ETK 

This subroutine computes the expected time for a given firer 
to kill a given target. The calculation is a function of range, time 
of flight for a round, and hit and kill probabilities for the firing 
weapon system. It is a number that is used in computation of the 
deterministic attrition coefficients. 

17. Subroutine SORT 

This subroutine is used to sort targets in ascending range 
order. This is used to determine the priority of a target for fire 
allocation. 

18. Subroutine KOVER 

This subroutine determines what portion of a particular target 
is covered by the terrain between the target and observer. 
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This number is used both in the detection of the target, and in the 
attrition computation. 

19. Subroutine LOS 

This subroutine was written by Professor James Hartman, Naval 
Postgraduate School. It computes a percent of a target visible to a 
particular observer, given the location coordinates of both. 
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VII. DEFINITION OF VARIABLES IN COMPUTER PROGRAM 



A. VARIABLES USED IN THE SHIP-TO-SHORE PHASE 



A ( i ) 


Aggressor force attrition coefficients. 


ATFFS 


Amphibious Task Force Fire Support. 


B(i ) 


Defensive force attrition coefficients. 


CDSURV(i) 


Current strength of defensive Unit i: 

1 - Tank 

2 - ATGM 


CSURVE(i) 


Current strength of assault wave i. 


DA(i) 


Attrition rate for defensive Unit i due to 
the effects of ATFFS/TLF. 


DEFWTS 


Defensive Force Tactical Allocation Weights. 

i 


DELTA 


Aiming error caused by the suppression factor 
of ATFFS. 


DINIT ( i ) 


Initial strength of defensive Unit i. 


DS1 


That portion of the defensive force ATGM 
unit assigned to engaging the closer of two 
multiple waves in the ATGM engagement window. 


DS2 


That portion of the defensive force ATGM 
unit assigned to engaging the farther of 
two multiple waves in the ATGM engagement 
wi ndow. 


DTI 


That portion of the defensive force Tank 
unit assigned to engaging the closer of two 
multiple waves in the Tank engagement window. 


DT2 


That portion of the defensive force Tank 
unit assigned to engaging the farther of 
two multiple waves in the Tank engagement 
window. 


DT1PH 


Hit probability of rounds fired by DTI against 
the assault wave in its engagement window. 
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DT1R0F 


Rate of fire utilized by DTI against the 
assault wave in its engagement window. 


GAINL 


Defender attrition level at which remaining 
defending forces withdraw and land combat 
commences . 


GALF 


Denotes whether the landing force buildup 
is sufficient for land combat: 

0 - Insufficient 

1 - Sufficient 


GAMMA 


Aim-reload time suppression factor. 


GAIK 


Denotes whether the landing force has 
initiated the land coir-bat: 

0 - Not started yet, 

1 - Started already. 


GATM 


Time at which land combat commenced. 


IL(i) 


Denotes if wave i has reached the shore: 

0 - Wave i not ashore, 

1 - Wave i ashore. 


IPRIW 


Denotes whether the user desires battle 
summary at each time-step, or just a 
final summary: 

0 - Battle summary printed after each 

time-step , 

1 - Final battle summary only. 


IWP:! 


Weapon-type code: Tank = 1, ATGM = 2. 


IWSTAT(i) 


Current status of assault wave i: 

0 - Not engaged, 

1 - Landed, 

2 - Under fire by ATGM, 

3 - Under fire by Tank, 

4 - Under fire by both ATGM and Tank. 


RD 


Distance offshore at which waves initiate 
their transition. 


RKSURV(i) 


Concatenation of CSURV and CDSURV. 


SA( i ) 


Attrition rate for wave i due to ATGM. 


SARTM 


ATGM aim-reload time. 


SENG(i) 


Wave number of the closer of two assault 
waves in the ATGM engagement window. 
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SENGMN 


ATGM minimum engagement range. 


SENGMX 


ATGM maximum engagement range. 


SRNG(i) 


Firing range to wave SENG(i). 


SSIGH 


The standard deviation error in the 
horizontal axis for ATGM fire. 


SSIGV 


The standard deviation error in the 
vertical axis for ATGM fire. 


SVEL 


ATGM projectile velocity. 


SWTS(i) 


The proportion of the total defensive 
force ATGM strength to be allowed to 
engage wave SENG(i). 


TA 


Time first assault wave initiates its 
transition. 


TA(i) 


Attrition rate for assault wave i due 
to Tank fire. 


TARTM 


Tank aim-raload time. 


TB 


Time first assault wave completes its 
transition. 


TBW 


The intsrarrival time between waves 
arriving at the beach. 


TFF 


Time first assault wave reaches the beach 


TENG(i) 


Wave number of the closer of two assault 
waves in the Tank engagement window. 


TENGMX 


Tank maximum engagement range. 


TMEANH 


The bias error in the horizontal axis 
for Tank fire. 


TMEANV 


The bias error in the vertical axis 
for Tank fire. 


TRNG(i) 


The firing range to assault wave TENG(i). 


TSIGH 


The standard deviation error in the 
horizontal axis for Tank fire. 


TSIGV 


The standard deviation error in the 
vertical axis for Tank fire. 
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TSURV 


Total number of surviving LVA ashore at 
the current time. 


TVEL 


Tank projectile velocity. 


TWTS(i) 


The proportion of the total defensive 
force Tank strength to be allowed to 
engage wave TENG(i ) . 


WB(i) 


Aimed-fire attrition-rate coefficients 
for defensive force Tank and ATGM assets. 


WID 


Width of LVA. 


WVINT(i) 


Initial strength of assault wave i. 


WVRNG 


Firing range to assault wave i. 



B. VARIABLE USED IN THE LAND COMBAT PHASE 



ALPHA(i) 


Initial attrition-rate coefficient for 
stochastic attrition option. 


ANGH(i) 


Orientation angle of the hill ellipse 
measured in degrees counter-clockwise 
from East to the major axis. 


APOA(i.j) 


t h 

The average proportion of the j tn attacker 
of Unit i allocated to fire on Unit i. 


AVD 


Average distance. 


AVSP 


Average speed of moving aggressor units. 


BASE 


Overall terrain elevation above sea level. 


BREAK 


Breakpoint distance between aggressor units 
and defensive units. 


DISMAX 


Maximum distance allowed between aggressor 
units before the leading unit is delayed. 


DIST 


The straight-line distance between two 
movement nodes input by the user. 


DST 


The distance in meters to be moved each 
time-step by aggressor units. 


ECC(i) 


The eccentricity defined as the ratio of 
major axis length to minor axis length. 



93 



FL(i) 


Force level of Unit i. 


FO(i ) 


Initial force level of Unit i. 


IALT 


Denotes whether the user desires alternate 
defensive positions or not: 

0 - Yes, 

1 - No. 


IC 


Counts number of time units a defender 
has been moving. 


IDIR 


t h t * 

Direction of j n interval in i tn route. 


II (i ) 


Interval index for Unit i. 


I ITIME 


Current time. 


I MAX 


Maximum number of time intervals allowed. 


IMOVE 


Number of time units a defender is allowed 
for moving to an alternate position. 


IPRDIR(i) 


Primary direction of fire for defensive 
Unit i. 


IRAN 


Range . 


IRTE 


Denotes whether user wants to input routes 
or not: 

0 - Program determined routes, 

1 - User determined routes. 


ISE 


A switch variable set to 1 when the defen- 
sive force ATGM unit initiates its fire. 


ISECWD(i) 


Width of search sector for defensive Unit i 


ISPD 


Input variable to denote user's desired 
speed for aggressor force movement: 

1 - 9 MPH, 

2- 12 MPH, 

3- 15 MPH, 

4- 18 MPH. 


IT 


Current time period. 


ITE 


A switch variable set to 1 when the defen- 
sive force Tank unit initiates its fire. 


ITEM 


Input variable denoting number of time- 
steps allowed for aggressor unit move. 
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ITIME 


Current time, in seconds, of battle. 


ITRIT 


Input variable denoting whether attrition 
will be stochastic or deterministic: 

0 - Stochastic, 

1 - Deterministic. 


IUSTAT(i) 


Current status of Unit i: 

0 - Alive, not firing, 

1 - Alive and firing, 

2 - Killed, 

3 - Moving. 


LVAFR 


Firing rate for LVA weapon system. 


LATOB 


Indicator variable for one- or two-way 
LOS calls: 

0 - Do not compute LOS from Unit A 

to Unit B, 

1 - Compute LOS from Unit A to Unit B. 


LISTH(i) 


List of hill numbers for each grid square. 


LOA(i.J) 


f* h 

The number of the j u ' target of Unit i. 


LOST ( i ,j) 


Denotes whether 1 ine-of-sight exists 
between Unit i and Unit j. 


LOT ( i ,j) 


t h 

The number of the j tri target of Unit i. 


LST 


Index number for the first hill listed 
for grid square (i,j) in LISTH(i). 


MVTDIR(i) 


Movement direction of Unit i. 


N(i ) 


Number of nodes inputted by user for 
route i . 


NA(i ) 


Number of aggressors of Unit i. 


NBU 


Number of defensive force units. 


NOVELS 


Number of forest ellipses in terrain. 


NF(i) 


Number of time units a Unit i is allowed 
to fire at the same location. 


NHL ( i ,j ) 


Number of hills in each grid square (i,j). 


NHILLS 


Number of different hills to be modeled. 


NHTOT 


Total number of hills modeled on battlefield 
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NLOSC ( i ,j) 


Number of continuous time-steps that LOS 
does not exist between Unit i and Unit j. 


NOD 


Number of time intervals Unit i delayed 
in movement. 


KOI(i) 


Number of intervals in the i tn route. 


NRU 


Number of aggressor force units. 


NT ( i ) 


Number of targets of Unit i. 


CFL(i) 


Force level of Unit i during previous 
time- step. 


P (1 J) 


Probability of first round hit by Unit i 
in range band j. 


FHH(i.j) 


Probability of a hit following a hit by 
Unit i in range band j. 


FHM(i.j) 


Probability of a hit following a miss by 
Unit i in range band j. 


P KH ( i ,j ) 


Probability of a kill given a hit by 
Unit i in range band j. 


PM 


The proportion of time a moving unit is 
searching for targets. 


FOA(i ,j) 


t h 

The proportion of the j n attacker of 
Unit i allocated to fire on Unit i. 


FOL(i) 


Percent of Unit i lost during the current 
time-step . 


FTT(i) 


Proportion of surviving firepower allocated 
to the i tn target if there are j targets 
available to be engaged. 


RANGE 


Current minimum distance between aggressors 
and defenders. 


RKATTR 


Vector containing the current attrition 
loss rates to be applied within the 
Euler integration routine. 


RF 


Detection rate reduction factor for a 
firing unit (in comparison to a non- 
firing unit). 


RMINTK 


Minimum effective range for an LVA 
mounted weapon system. 
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RMINTW 


Minimum effective range for a TOW weapon system. 


RMXTK 


Maximum effective range for an LVA mounted 
weapon system. 


RMXTW 


Maximum effective range for a TOW weapon system. 


ROF 


Rate of fire. 


ROT(i.j) 


Range of the j^* 1 target of Unit i. 


SIZETK 


Size of LVA weapon system. 


SIZFTW 


Size of TOW weapon system. 


SPRD(i) 


Measure of hill size which is defined to 
be the distance in meters measured along 
the major axis from hill center to the 
contour line which is 50 meters down from 
the peak. 


SUMBO 


Total cefensive force level. 


SUMRO 


Total aggressor force level. 


SUP FAC 


Suppression factor. 


TA(k) 


Tine tc acquire a target for k tn weapon 
system type (k = 1 , 2) . 


TFl(k) 


Time of flight to 1000m for k^ weapon 
system type (k = 1 , 2) . 


TF2(k) 


Time of flight to 2000m for k tfl weapon 
system type (k = 1 , 2) . 


TF3(!<) 


Time of flight to 3000m for k 1 "* 1 weapon 
system type (k = 1, 2). 


TH(k) 


Time tc fire a round following a hit 
for weapon system type (k = 1, 2). 


T I ( k ) 


Time to fire first round after target has 
been acquired for weapon system type (k = 1,2). 


TM(k) 


Time to fire round following a miss for 
weapon system type (k = 1, 2). 


TMACI ,TMACJ 


Elevation of Unit i and Unit j in LOS model. 


TOWFR 


Firing rate for TOW weapon system. 
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TPOL(i) 


Total percentage of lost since battle began 
for Unit i . 


VISFR(i,j) 


The fraction of Unit i as seen by Unit j. 


VISFRA 


Fraction of Unit A as seen by Unit B. 


VISFRB 


Fraction of Unit 8 as seen by Unit A. 


X(i),Y(i) 


Coordinates of Unit i. 


X A ( i ) , Y A ( i ) 


Coordinates of alternate position for 
defensive Unit i. 


XC(i),YC(i ) 


Coordinates of center of hi 'f 1 i. 


X I C ( i , j ) 
Y I C ( i , j ) 


t h 

Coordinates of the j 1 interval endpoint 
of the route for Unit i. 


XL,YL 


Distance added to previous interval endpoint 
for vehicle to move DST during a time-step. 


XLOC(i,j) 
YLOC(i J) 


Coordinates of the node inputted by 

the user for the route of Unit i. 



98 



APPENDIX B 



COMPUTER PROGRAM 
for the 

SMALL-UNIT AMPHIBIOUS OPERATION COMBAT MODEL 



The small-unit amphibious operation combat model is a computerized 
model written in FORTRAN. It consists of a main program and 19 sub- 
routines. It was designed to serve as a reference to itself in order 
that the reader would not be forced to refer to various manuals outside 
of the program each time an explanation of the functioning of a parti- 
cular aspect of the program was desired. A listing of the computer 
program follows. 
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nnn ono nn nnno 



THIS PRCGPA.M IS A SMALL-UMT AMPHI8ICUS OPERATION COMBAT MODEL 
UTILIZING LANCHEST ER-T YPE EQUATIONS 70 COMPUTE ATTRITION. 

IT CONSISTS GF T*C BASIC PHASES. THE FIRST BEING THE ShIP-TO-SHOR 
COMBAT PHASE . ANC THE SECOND BEING THE LANO COMBAT PHASE. 

***** SHIP-TO-SHCRE PHASE CCMMON BLOCK VARIABLES ***** 



COMMON /AMPH/IL(5 ) , WB ( 2 ) ,A(2 1.6(2). ITE.ISE. RD» WV INT ( 5 I , WID » 
*T8W, DINIT(2) ,GAINL,IV>STAT(5) 

COMMON / ENC-R/ S PCM AX , SPDM IN , HTMAX ,HT M I N, TTS, TA A , T8 , TFF 
COMMON /DISPER/TS IGV (6 ,2 I ,TSI CH(6 ,2 ) » TME AMH ( 6 . 2 I , 

*SS IGV( 7, 2|,SSIGH( 7,2) 

COMMON /CcR/TENGMX. SENGMX .SEN GMN, TAR TM ,S ARTM , T VE L , 

*SVEL,DEFWTS( 2) 

COMMON /SUPEFT/GAMMA, DELTA 
COMMON /IOUT/TSUFV, ICR INT 

***** LANC COMBAT PHASE COMMON BLGCK VARIABLES ***** 

CCMMON /GRP1 / IPPDIR <6 I , I <ECWC(6) ,MVTDIR<6) ,X<6) ,Y<6 ) ,S?Q(6 I 
COMMON /GRP2 / TA(2 ) ,T1(2 ) ,TH( 2 ).TM( 2 ) ,TF 1( 21 ,TF2{ 2) , TF2(2 I . 
*P(2,6),PFH(2,6),FHM(2,6.),FKH<2,6),TF(2) 

COMMON / C- R P3 / NBU , MRU , FL ( < ) , F C ( 6) ,N0 I ( 3) , XI C ( 3 , 200 ) , Y I C ( 3 , 200 ) , 
*I0IR(3,2C0>, AVSF, ISPC 

*, IUSTAT (A), 11(E) , LOST ( 6. 6 ) , V I SFRA ,VI SF SB , SI ZET K , 

*S IZE T W, NT ( 6) ,NF( 6 ) ,SRF,DI SMAX , 

* N L C SC ( 6 ,6 I ,V I S F R ( 6 ,6 ), RM I HT K , PMXTK, RM INT W , RMX T W ,OP, T OWCR,LVAFR. 
*PTT (3,3) , RF, PC A ( 6 , 6 ) ,AFO A (6,6) ,L0A(6 , 6 ) , NA ( 6 ) , OF L ( 6 ) , POL ( 6 ) 
COMMON /GRF4/ TFCL ( 6 I , CL DC < 6 , c I ,Q ( 6, 6 ) 

COMMON /G 0 F5 / LOT (6 ,6) ,R01<6, 6 I 

COMMON /HILLS/ XC ( 1C C ) , YC ( 1 U C ) ,PE AK ( 100) , ANGH ( 1 00 I , S PRO ( 1 00 I 
COMMON /HILLS/ SC C < 100 I , P >X ( 1 CO I , PY Y ( 1 CO ) , PXY ( 100 1 , BASE 
COMMON /HILLS/ NH ILLS . 

COMMON /COVER/ CXC ( 1 50 I , C YC < 1 5C ) , CPE 5*K ( 1 50) , CP XX ( i 50 ) ,CPYY( 15 0) 

COMMON /COVER/ CPX Y ( 150 I , HCV c L S 

COMMON /CCLNTR/ KH.Khln.KV , K N , H G R S , KE L L , K I NT 

COMMON /CR ID/ LSI (5,4) , NHL( 5, A) ,L ISTH( 150) ,KHREP( 150) .KTREP 
COMMON /GRID/ LSTC ( 5 ,A I , N C ( 5 , A ) , L 1ST C ( *00 I ,KCR EP ( 150 ) 

COMMON /GRP6 / ALPHA(6I 

COMMON /C RP7 / XA ( 6 ) , Y A ( 6 ) , I MC VE ( 6 ) 

***** MAIN ORIVER PRCGRAM ***** 



GA TM=0. 

GATK=0. 

C *** INITIALIZE DATA FOR SH I P- 10-SH CRE PHASE 
CALL DATAIN 
C 

CALL SETLP 

C *** CCNCUCT SH IP-TC-SHCRE COMEAT PHASE 
CALL SEA (GATM, C-ATk ) 

IFtGATK.NE.O.) GC TC 5 
WRITE (6,600 ) 

STOP 

5 IF(GATK-Z.C) 1C.2C.3C 
10 KR IT E (6 ,6 10 ) 

GC 1C AO 

20 WPITE(6,62C) 

GC TC AO 

3 C API TE ( 6 ,63C ) 

AO WRI TE (6 ,oAC ) CATM 

C *** CCNCUCT LAND CCMEAT PHASE 

CALL GRCLNDt GA 1M , T SLR V ,1 PRINT , TTS ) 

C 

C 

60C FCRMATl IX, 'TOTAL LANDED LANDING FCPCS STRENGTH IS INSUFFICIENT', 
* ' FOR GROUND ATTACK ' ) 

610 FCRMAT( 1 X , 'LANC COMBAT STARTS WHILE SHCRE COMBAT IS GOING ON') 
620 FCRMAT( lx, 'LANC CCMEAT STARTS AFTER CEFENDER BREAKS C CNT ACT • ) 

630 FCFMAT ( IX , 'L ANC CC^BAT STARTS A C T ER ALL WAVES LANDED') 

6 AC FCPMATt /IX,' LAND C C M B A T ATTACK TIME =',F6.1,' SECCNCS') 

STCP 
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SLBROUT INE SEA (GA7M ,GATK ) 

*** SLEROUTINE IS THE VAIN DRIVER PROGRAM FOR THE SH I P- TO SHORE 
PHASE OF The AMPHIBIOUS OPERATION. ITS MAIN PURPOSE IS TO 
INITIALIZE KEY PARAMETERS AND TO DIRECT PROGRAM FLOW FOR THE 
SUP-TC-SHCRE PHASE CF COMBAT 

CCNMCN /AMPH/I L( 5 ) ,WB{ 2) , A( 2) , B < 2 ) , I TE , I $E * RD , WV INT ( 5 ) , WID , 
*TEVtDINn(2) ,G AIM, I US TAT (5) 

COMMON / ENGR / $P CM AX , SPD M IN , H TMAX *HT M I N, TT S , T A A , T 9 , T F F 

CALL OUTFIT 
IRC= 500 
I T EW = 1 2C 
RC=1.0*I PC 
TBW=1.0* ITBW 
TINT =0.0 



+** COMPUTATION OF FIRST WAVE TIME PARAMETERS 
TA-TIME FIRST WAVE INITIATES TRANSITQN 
TB-TIME FIRST WAVE COMPLETES TRANSITION 
TFF-T IME FIRST WAVE REACHES THE BEACH 

TAA=( 5C0C.-RD J /SPCMAX 
TE-TAA+TTS 

TFF=TB+ ( FC-( 0. 5* ( S PCMAX- S FOM I N )*TTS ) -150. )/SPDMIN 
D E L= 10 • 

W F IT p (6*600) RC,TEU 

600 FCFMATi/.IXt 1 ITERATION I N IT I ATED. .. R D= « , F10. 3, 1 X , * T8W = 
*• , F10.3 ) 

CALL RKINTfDEL ,TINT f N*GAT M,GATK ) 

RETURN 
EN C 



SUBROUTINE RKINT(H,TI*N,GATM f GATK) 

*** SL6R0UT INE P.KINT PROVIDES THE INTERFACE BETWEEN 
THF EULER NUMERICAL INTEGRATION ROUT I \E( RKLDEQ ) 

ANC THE SUBROUTINE ATTR WHICH DETERMINES EACH 
UNIT’S STATUS AS TIME PROGRESSES THROUGHOUT THE 
AMFHIBICLS OPEPATICN 

COMMON /AMPH/I L ( 5 ) tWB<2) » A(2 ) , 6(2 )f ITS* ISE»RD» WV INT( 5 ) tWIDt 
*TBWtDINI7( 2) tGAIN.Ltl wSTAT (5) 

COMMON /ICUT/TSURVt IPRINT 

DIMENSION CSURV (5 ) .C CSURV (2) , T A < 5 ) , S A ( 5 ) , DA ( 2 J 
DIMENSION RKSURV( 7) *RKATTR< 7 ) ,TA T TR (200 ,12) ,TI ME (200) 

***** VARIABLE DEFINITIONS ***** 

I v AX - MAXIMUM ALL CWAdEL NUMBER OF TIME INTERVALS 

nm - a switch variable whose element i is $e t to i when 

RAVE I ARRIVES AT THE BEACH 

ise - a switch variables set tc i when the def.atgm 

UNIT INITIATES ITS FIRE 
IT - CURRENT T IME PE^ ICO 

ITE - A SWITCH VARIABLES <ET TC 1 WHEN THE DEF.TANK 
UNIT INITIATES ITS FIRE 
T - CURRENT TIME 

TSLRV - TOTAL NUMBER OF SIRVIVING LVA AT THE CURRENT TIME 

***** STATE VARIABLE DEFINITIONS ***** 

C CSURV ( I ) - CURRENT STRENGTH CF DEFENSIVE FORCE I 
I = 1 TANK 

I = 2 A 1GM 

CSURV { I ) - CURRENT STRENGTH OF ASSAULT WAVE I 
DIMT(I) - INITIAL STRENGTH 0^ DEFENSIVE FORCE I 
RKSURV(I) - CONCATENATION CF CSURV ANC CDSURV 
WVINT(I) - INITIAL STRENGTH OF WAVE I 
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10 

15 

20 

25 

20 

120 



40 

45 



46 

50 



110 

120 

65 



I> = 235*2 
G A l F=0 • 

I . w £X=1 9 9 
I T E = 0 
1 5 E = 0 
T S L R V = 0 • 

71^E(II-C* 

T = T I 

DC 10 1=1 ,5 

CSURV(I ) = W V I M T ( I ) 

TSURV = TSUPV+CSLFV(I ) 

I L ( I ) = C 
JWSTM(I) =C 
CONTINUE 
DC 15 1=1,2 

CDSURV ( I ) =C IN IT ( I ) 

CCMINUE 
CC 20 J = 1 , 12 
T A T 7 P ( i »'J ^ - 0 • 

17 = 1 

CC 25 1=1,5 
PKSUPV (I ) “ C S 0 R V { I) 

R K SLR V ( 6 ) =C S LR VI 1 ) 

PKSURV (7 )=CSURVf 2 ) 

CC 30 1=1,7 
RKAT7P ( I )=C. 

M*0 

CALL A T7F (7, CSLRV ,CCSURV ,7A,S A ,CA,GA LF , GATK, GATM , IX ) 

***** S T AT E VARIABLE DEFINITIONS ***** 

D A ( 1 ) - ATTRITION RATE FOR DEFENSIVE UNIT I DUE 
7 HE EFFECTS CF ATFFS/TLF 
S A ( J I - ATTRITION PATE FOP WAVE I DUE TO ATGM 
TA(I) - ATTRITION RATE FCP WAVE I DUE TO TANKS 



TO 



PXATTR(I) IS A /ECTCF CONTAINING THE CURRENT ATTRITION 
LOSS RATES TC Be APPLIED W I T H I N THE EULER 
INTEGRATION ROUTINE TO THE STATE VAPIABLES. 
1=1,5 LV A WAVES 1-5 

1 = 6 D T 

1 = 7 OS 



I> = IX * 7 

IF ( I L ( 1 ) .EC. 99 ) GC TC 120 
00 4 C 1=1 ,5 

PK5URV (I ) = C S LR V ( I ) 

RKATTF (I) = (TA < I )+SA( I) )*(-1.0) 

CC 45 1=1 2 

’KSUPV i 1+5 )=CCSURV( I) 

R M T7R ( I 4 5 ) = -l .O^CA( I) 

S =PKLDE C n fRKSLRV ,RKATTR , T,H ,NT1 
LS INK = 1 

DC 50 1=1 ?5 

IF ( L5INK.EC.C) GC TO *6 

CSLRV(I) = CSUPVU) + RKATTR ( I ) 
GC TC 50 

CSURV (1 I = PKSLR V ( I ) 

CONTINUE 
DC 55 1=1,2 

CDSUPV x I ) =RKSUR V ( 1 + 5) 

CONTINUE 

I F ( S- 1 • ) 110 , 12C, 12C 
WRITE (6 ,600 ) 

STOP 

CONTINUE 
I T = IT + 1 
TSLRV=0. 

CC 65 L = 1 , 5 

TSLRV=7 HPV+CSLPV (U 

IFtTSURV.LE. C. ) T SLR V = 0* 

T I M E ( JT ) = T 

IFdPRIM.EO.l > GC TC 999 
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C *** PRINT RESLLT CF SH I P- T C- SHORE MOVEMENT AFTER EACH TIME STEP 
WPITE(6*61C5 T 
W F ITE ( 6 *620 
CC 70 I = 1 * A 

PL0ST»1.-CSLRV( I ) /WVINTf I ) 

WRI TE (6 ,6 30 ) I »CSU?V< I) f IWSTAT (I ) f OL0ST 
7C CONTINUE 

FLCST = i.-CSURV (5 )/V*V INT ( 5 ) 

V* RITE (6 ,640) C SUFV (5 ) , IWSTAT ( 5 5 * P LOS T , TSURV 
PUCS T =1.-CC$URV( 1 I / 0 I N I T ( !> 

W FITE (6 *6 50 CCSURV ( 1 J * PL CST 
FLCST= l.-CCSUR VI 2 )/CiNI T l 2) 

TASURV=CCSUR V { 1)*CCSURV( 2 ) 

V% F I T F ( 6 *660 CCSLRV( 2 > * P l CST * T AGU RV 
595 CONTINUE 



*** DETERMINE R: THE FIRING RANGE TO THE LAST INCOMING ASSAULT WAVE. 



R = PNG( T - A « *T BW ) 

DETERMINE IF ALL WAVES LAMJcC A NO LA N j COMBAT STARTED 
NCTE : THE MODEL IS TERMINATED IF: 

1. THE FIRING R A ^ GE TO THE LAST ASSAULT WAVE IS LESS 
THAN 75 METER i. 

2 . THE DEFENSIVE BREAKPOINT HAS BEEN REACHED 

3. THE MAX I MU F NLMBEF Or ITERATIONS HAS BEEN EXCEEDED 



I F (R • LT .75.) GC TC 2C0 

I F ( IT .GT. IMAX ) GO TC 2 CO 

IF ( I Li 11 .EC. 99) GO TC 2CC 
GC TO i 2C 

2C0 N = I T 

*** PRINT R ESUL~ OF THE SH I P- TO- SHORE COMBAT PHASE 
WFITE(6,6iO T 
W F I TE ( 6 , 620 
CC 90 1=1, A 

flo$t=i.-cslrv(ij/wvint( n 
WRITE (6,6 30 I , C SLR V< I ) , I W S T AT ( I I ,FLCST 
90 CONTINUE 

FLCS T = 1 . -CSUP.V (5 )/WVIN7( 5 ) 

WP ITE ( 6 , 6 ACi C S ' j F V ( 5 ) *IWSTAT( 5 5 * ? LOS T , TS JR V 
PLCST=l.-CC3URV(i ) /C IN IT ( I) 

WRITE (6 *650 CCSLFV(l) * 3 l C ST 
F L CS T = 1 • -CDS UR V ( 2 ) /O 1 N I T ( 2) 

T ASUPV=CCSuRV ( 1 ) + C C S L R V ( 2 ) 

WFITE( 6, 6cC) ClSlRV(2J,PLCST*TASURV 
W P I T E ( 6 , 670 T SLR V 
I F (GATK. G6.1 . ) GC TC 9999 

IF C TSOP V. LT.9. ) GC TO =999 
GATK=3 . 

GAT M=T 

9999 CONTINUE 

600 FORMAT (IX , ' ERRCR..S.NE.1.CR.2' ) 

61 C FORMAT ( / / / IX , ' T I M E - ’ **=6.1 *lX f 1 SECONDS' // ) 

620 FORMAT ( IX* 'WAVE' , 2 X , 4 c CR C E L E VEL 1 ,2X , ' ST A TUS 1 , 2X , ' LCST-PCT' * 
*2X * 'TOTAL SURV IV INC ' ) 

630 FORMAT ( 3X , I 1 , 2 X , F 1 C .4 , EX , I 1 ,5X ,F 8 . 3 ) 

6AC FCFMAT ( 2X. '5 1 , 3>*F 10.4 ,5X ,1 1 , f X, = 0. 3 , 7X* F5. 2) 

65 0 FORMAT (IX f ' 7 AN K ' ,2X , c 10«-* ,UX , F6. 3) 

66 C FCFMA7 ( 1 X * ? A T G M 1 ,2X*F1C.A ,1 1 X , r 0. 3 , 7 X , F5 . 2 I 
670 FCFMAT ( IXt »F INAL LVA SURVIVORS AS HCR E= »*Fi0.3) 

RETURN 

END 



FUrCTICN RKLDEOh *7, F,X, H *NT ) 
OIMENSICN 7(1) *F ( I ) * C ( 25 ) 

NT=NT+l 

GC TO (1 *2,3 *4 ) , NT 
1 hl*H 

A A = H1/ A . C 
CC 11 J =1 * N 
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11 Q(si)=0. 

X=X+AA 
GC TO 5 

2 XO + AA 
GC TO 5 

3 X = X+A A 
GC TO 5 

4 CC 6 L = 1 ? N 

6 Y ( l ) = Y ( t ) 4 AA* F (L ) 
N 7 = C 

X=X+AA 
F K IDEC =2 . 

GC TO 6 

5 CC 7 1 = 1 ,N 

7 V ( I ) =Y ( I ) + AA*F<I) 
R K10EQ= 1 • C 

8 RETURN 
ENC 



SUBROUTINE A TTR ( T , CSUR V, 0 EUR V , TA , SA t CA »G A LF »GA T !< t GAT v » I X ) 

*** GIVEN THE CURRENT TIME ANC STATE VARIABLE STRENGTHS n 

SUBROUTINE ATTR CETERM INE 5 THE ATTRITION RATES AND UPDATES 
THE STATLS OP EACH UNIT WITH FESPECT TC SHORE MOVEMENT 
ANC IMPLEMENTS THIS INFCR NAT I CN INTO THE ATTRITION LCSS RATE 
CCMPUTAT 1CN. 

***** STATE VARIAELE DEFINITIONS ***** 

C A ( I ) - CURRENT ATTRITION LOSS R A TE FOR OEF • FORCE I DUE TC 

ATFF S ( AMFH I e IOUS TASK FORCE FIRE SUPPORT )/TLF EFFECTS 
IKI) - i-HEN ECLAL TC 99 INDICATES THE DEFENSIVE BREAKPOINT 
HAS BEEN PEACHED 

S A ( I ) - CLRRENT ATTRITION LCSS RATE FCF WAVE I DUE TC A TGM FIRE 
TAU) - CLRRENT ATTRITION LCSS RATE FCR WAVE I DUE "C TANK FIRE 



CCMMON /AMPH/IK5) ,W8( 21 t A<2) ,6(2), I TE , I SE , RD , W VI NT i 5 > ,W'ID, 
*T EV.DIN IT ( 2) • G A INLt IWSTAT (5 ) 

COMMON /CEF/T6NGNX ,SENGMX ,S ENGMN, TAR TM ,S ARTM,TVEL , 
*SVEL,DEFWTS< 2) 

INTEGER T ENG ( 2 ) ,SENG<2 ) 

DIMENSION T RNG (2) ,TWT$ (2) ,SRNC(2) ,DSLRV(2) ,SWTS(2) , 

*CSLRV (5 ) ,TA< 5 ) ,SA (5 ) ,CA<2 ), ASX (20 ) 



LSINK = 1 



CC 10 1=1 ,5 
TA < I I = C • 

S A ( I )=C. 

1C CONTINUE 

***** VARI/8LE CEFINITIQNS ****** 

DTI - THAT P CR 71 CN CF THE DT UNIT ASSIGNED TO ENGAGING THE CLOSER 
OF TWO MULTIPLE WAVES IN THE TANK ENGAGEMENT WINDOW 

DT2 - THAT PCRTICN CF THE DT LNIT ASSIGNED TO ENGAGING THE FARTHER 
OF TWO MLL7IPLE WAVES IN THE TANK ENGAGEMENT WINDOW 

0S1 - THAT PCRTICN OF THE OS LNIT ASSIGNED TO ENGAGING THE CLOSER 
OF TWO MLLTIPLE WAVES IN THE A TGM ENGAGEMENT WINDOW 

CS 2 - THAT PCRTICN OF THE DS LNIT ASSIGNED TO ENGAGING THE FARTHER 
OF TWO MULTIPLE WAVES IN THE ATGM ENGAGEMENT WINDOW 



C$1=0. 
CS2=0. 
CU=0. 
CT2=0. 
F A C =1 • 



0 
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*** CETERMINE IF PART CF LAND ING FORCE ADVANCE TO ATTACK INLAND 
KEY TERR/ IN 

I F (GATK. EC #1 • 0 ) GC TC 15 

IF (GALF.E C.l.C . /NC. (CSLRV ( 1 )+DSUR V( 2) ) . LE. G A I N L* ( D I N I T ( 1 ) 

* + C IN IT ( 2 ))) G A 7M = T 

IF (CALF. EC. 1.0. AND. (DSURV( D + OSURV ( 2> > . LE. GAINL* (DINIT < 1 ) 
*+CINIT(2))> GA7K=i.O 

*** CETERMINE IF D EF • BREAKPOINT FAS BEEN REACHED 

15 IF((DSURV(1)+DSUFV(2)).LT.0.2*(DINIT ( 1 ) +D INIT( 2 ) ) ) GC TO 20 

*** CETERINE ATTRITION RATE Ci^ DEFENSIVE FORCES 3Y ATFFS 
BASED UPCN AREA CR AIMED FIRE STATUS 

► VRNG = FIRING R A f^GE TC AN ASSAULT WAVE 

c/«i) = sm 

CA <2) = B (2 ) 

IF( ITE.EC.C) DA<i)=A(l)*05URV<l) 

I F( ISE.EQ.O) C A (2 >=A( 2 )±DSLRV< 2) 

GC 1C AO 
20 CSLRV( 1 

DSL RV ( 2 ) -C • 

D A ( 1 ) = C • 

D A { 2 ) = 0 . 

IF (GATK. EC. 1 .) GC TO 35 
GAT = T 



*** CETERMINE !<= C E F. ER E AK PO I M HAS BEEN REACHED BEFORE SUFFICIENT 
LANDING FCRCE IS BLILT UP ON THE SHOFE FCR INLAND ATTACK 

25 DC 30 I=lt5 i 

WVRNG^RNG (G AT-TEWM I - 1 ) ) 

IF ( WVFNG. LT .75. ) IL ( I ) =1 

IF { I L ( IJ.EC.l) T IF = 1LF + C SUR V( I ) 

30 CONTINUE 
G A7=GAT + 1C. 

IF (TLF • L 1 • 9. 0 • ANC . IL<5). EC.l ) RETURN 

IF(TLF.LT .9.0. AND. IL( 5 l.NE. 1) GO TC 25 
G A 1 K =2 . 

GAIF«1. 

GAT M=G AT 

WRI7c(6»61C) GAT V 

610 FCFMAT (/ f lX t 'LAND C CMBAT INITIATED 

35 I L (1 )= 99 

WPITE( 6 , 62 C ) T 

620 FCFMAT ( IX, ‘BPEAKPOI NT REACHED AT TIME =* ,F7. 1 , * SECONDS 1 ) 

RETURN 

*** SLEROUTINE JTGTS CETERMINES THE FIRING STATUS FOR THE 
TWC DEFENSIVE UNITS. 

AO CALL DT GTS (T f T ENG , TRNG t TW TS » S ENG. $RN G. SW TS ,C SURV ) 

***** STATE VARIABLE DE c INI TIONS *•**■** 



T IME = * . F7 • 1 , 1 SECONDS 1 ) 



TENG(l) - THE WAV E NUMBER OF THE CLOSER C c TWO 
WAVES IN THE TAK ENGAGEMENT WINDOW 
THE FIRING RANGE TC WAVE TENG(i) 

THE FFCFCP T ICN CF THE TOTAL OT STRENGTH 
TO BE ALLOCATED TC ENGAGING TENGdl 
TENG ( 2 ) - THE WAVE NUMBER OF THE FARTHER OF TWO 
WAVES IN THE T A f K ENGAGEMENT WINDOW 
SIMILAR I N T ERP ROTATION AS TRNG <11 
SIMILAR INTERPRETATION AS TWTS(l) 

THE WAVE NUMBER CF THE CLOSE? 2 F TWO 
WAVES IN THE ATOM ENGAGEMENT WINDOW 
FIRING RANGE tq WAVE SENG( 1 ) 

THE PROPORTION CF THE TOTAL CS STRENGTH 
TO BE ALLOCATED TO ENGAGING SEMG(l) 

S ENG ( 2 ) - THE WAVE NUMBER CF THE FARTHER OF TWO 



TRNG( 1 ) 
TWTS(1 ) 



TRNGC2I 
T W TS ( 2 ) 
SENG( 1) 

S PNG ( 1 ) 
SWTS( 1) 
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WAVES IN THE AT GM ENGAGEMENT WINDOW 
SRNC*(2) - SIMILAR INTERPRETATION AS SRNG(l) 

SWTS (2 ) - SIMILAR INTERPRETATION AS SWTS ( 1 J 

*** CETERMIN6 THE CUMULATIVE hUMBER OF SURVIVING LVA’S 
THAT HAVE REACFEC THE SEA CH - TLF 

T L F=0 • 

CC 45 J*I ,5 

IF( IL (J). EQ.l) 7Lf=TLF+CSURV(J) 

45 CONTINUE 

*** CETERMINE IF TLF BUILT UP IS SUFFICE NT FOR LAND CCMBAT 

I F ( TLF . G E • 9. ) GALFM. 

IF ( GATK.EQ. 1. ) TLF=TLF-9. 

*** ALLOCATE THE FORCE STRENGTH OF TLF BETWEEN THE TWO 
CEFENSIVE FORCE UMTS 

C SUM=0 S U F V ( I ) + C SLR V ( 2 ) 

TlFI=(OSLRV(i J /CSUMJ^TLF 
UF2=(DSLRV( 2) /DSLM)*TLF 

*** ACC TO C/1 ANC CA2 THE ATTRITION LOS S PATE DUc 
TC THE EFFECTS GF TLF I ANC TLF2 

OA(l)=OA(l)+TLFl*we(U 
C A ( 2 ) = 0 A (2 ) + TLF2*W0( 2) 

IF(DSURVU).LE.O.G) CA(1 )=0. 

IF(OSLRV( 2) .LE.C.C) DA<2)=0.0 
ICA = 1 
N = 20 
MUL = I 
ISCPT = C 

CALL LRNC< IX ,ASX,N, KJL, ISCRT) 

*** CETERMINE IF THERE EXISTS AN INCOMING WAVE IN THE 
TANK ENGAGEMENT WINCCW (I.E. T ENG ( 1 ) .NE. 0 ) 

I F (TENG ( 1 ) • EC* C • ) GC TO ICO 
ITEM 



*** CETERMINE THE TI*E SINCE WAVE TENG(l) CROSSED THE 
5 C CC. MEIER OFFSHORE MARK -T1 

71 =7-TeW4 (TENG (I )-l ) 

C73=TWT$(1)*0$LRV( 1) 

F A C = 1 • 



*** CETERMINE THE SUPPRESSION EFFECT ^0 oc imposed 
CN THE CT UNIT 5ASEC CN TEE ATTRITION LOSS RATE 
CLRRENTLY IN EFFECT 

SUPFAC = ATFFS SUFRESSICN FACTOR 

S L FF AC = 0 A ( 1) 

CALL RATE (TRNG (1 ) ,SP C( T1 ) ,1 , SUFFACt D TIPOF ) 

*** OTIROF - RATE CF FIRE UTILIZED 0Y OTi AGAINST WAVE 
CALL PUT (TRNG (1 ) »WIC»HT( 711 , l,SUPFACfDTlPH» 

*** DTI Ph - HIT FwCBAEI L I T > OF RCUNOS FIRtC BY OTI 
AGAINST WAVE T ENG ( 1» 

*** DETERMINE THE ATTRITION LOSS RATE FOR WAVE TENG ( 1 ) 

CLE TO CT 1 FIRES 

TA (TENG ( 3 ) )*DT 1 PH* DT IP CF* CT1 
IF (LS INK.EQ.O J GC TC 55 

IF (ASX(ICA).GT.TA(TENO(l ) ) ) GC T C 50 
7A(TENC( 1 J ) = 1.0 
GC TC 55 

50 TA( TENG (1 J ) = C • 0 



T ENG ( 1 ) 
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55 ICA = ICA ♦ 1 

*** DETERMINE IF THEPE IS A SECCNC INCOMING WAVE THAT 
IS IN ThE TANK ENGAGEMENT WINDOW* IF THERE IS THE 
ATTRITION RATE COMPUTATIONS ARE SIMILAR IN FORM 
TC THOSE PREVIOUSLY PERFORM EC FOR THE CLOSER WAVE 

I F ( TENG ( 2 ) • E C • 0 ) GO TO IOC 
T2=T-m*M T ENG (21-1 ) 

0T2sTVTS( 2 ) *D SlR V { 1) 

CALL FATE (TFNG(2) »SPO(T2) * 1 , SU PFA C f DT 2R0F ) 

CALL RHIT (TRNG (2) , W ID * FT ( T 2 ) » 1 ,SU PFAC t 0T2 PH ) 

TA(T ENG (2 ) ) = CT 2 PH *QT2 R C C *C T 2 
I F ( LSINK.EC.O) GO TC 65 

IF(ASX(ICA).GT.TA{ T ENG ( 2 ) ) ) GO TC 60 
1 A (T EN C ( 2 ) ) = 1.0 
GC TO 65 

60 TA (TENG( 2) ) = 0.0 

65 I C A = ICA + 1 

*** CETERINE IF THERE EXISTS AN INCOMING WAVE IN THE ATOM 
ENGAGEMENT WINCCW, IF THEFE IS, DETERMINE THE ATTRITION 
EFFECTS AC-AINST THAT WAVE OUE TO ATGM THE ATTRITION 
RATE COMPUTATION ARE SIMILAR IN FORM TC THOSE FOR THE 
EFFECTS CLE ThE TANK FIRE. 



ICC 



I F ( SENG (l).EC.C) CC TO 20C 
I SE = 1 

S I - T - T 6W-* ( S ENG ( 1 ) - 1 ) 

DS1 = SHS( 1) *DSLRV<2) 

SUPFAC=CA ( 2 ) 

CALL RATE(SRNG(1),SPD( SI) * 2 * SUP«=A C , DS 1R0F ) 
CALL FHIT (SFNG (1 ) ,W ID , hT(Tl ), 2,$U PFAC, DS1PH) 
SA< SENGd ! ) =DS 1PH-DS1R CF*DS1 
I F ( LSINK. EC.U ) GO TC 75 

IF<ASX( ICA) .GT.SA( $ SNG < 1 ) ) ) GO TC 70 
SA(SENG( 1 ) ) = 1. C 
C-C TC 75 

SA ( SENG ( 1) ) = C. C 

ICA = ICA + 1 

IF <SENG<2 ) .E C • C ) GC TC 20C 
S2= T-ieM( SENG (2) -II 
DS2=SW7S( 2 ) S 0 5 V ( 2 ) 

CALL FATE (SFNG (2) *S PC ( 52 ) , 1 , SU PFA C * CS 2R0F ) 
CALL PH IT ( SRNG ( 2 ) , w I D » H T ( S 2 ) ,2,SUPFAC *DS2PH) 
SA(SENG(2 ) |sCS2Fh’CS2PCF*CS2 
IF ( LSINK. EQ.G) GO TC 35 

I F (ASX ( ICA).GT.SA( S ENG ( 2 ) ) ) GO TC 80 
SA(SENG(2 ) ) = 1.0 
GO TC 85 

80 SA <$ENG< 2 ) I = 0. 0 

E5 ICA = ICA + 1 



7 C 

75 



2CC RETURN 
END 



SL0ROUT I NE DTGTS <T tTENGt T FNGt TWTStSE NGtSRNGt SWTS , CSURV ) 

*** GIVEN ThE CURRENT TIM= ANC LVA WAVE SURVIVOR POPULATICNSt 
SUBROUTINE DTGTS DETERMINES THE WAVE NUMBERS THAT ARE 
TC BE ENGAGED BY OFFENSIVE TANK AND ATGM UNITS BASED 
CN THE ENGEMENT WINDOW CRITERIA 

COMMON /AMPH/ I L( 5 ) *WB( 2) * A( 2) , E ( 2 ) , I TE ,1 SE , RD , WV I NT ( 5 ) »WI D, 
*T E W ,DI N IT (2 ) » GAINLt I* ST AT (5 ) 

COMMON / C EF/ TENGM > tSENGMX iSENC-MN, TARTM ,S ARTM,T VELt 
*SVEL,QEFW7S< 2 ) 

INTEGER TENG (2 ) tSENG (2 ) 

DIMENSION TRNG (2) , SRNG<2) ,TWT5<2) ,SWT$(2) , CSURV (5) ,DENC(5 ) 
DO 10 1=1,2 
TENG ( I ) =0 
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ThTS ( I ) =0 . 

TRNG ( I ) =0 • 

SRNGCI>«0. 

SENG ( I»=0 
SWT S ( I ) =0 • 

10 CCNTINUE 
JT=0 
JS=0 
TSLM^O. 

SSLM=0. 

CC 100 1*1,5 

WVRNC- = FNG ( T - TE If * ( I~l) ) 

IFCWVPNG.LT. 75. ) ILCI ) =1 

IF (VVRNG.LT. 75. ) I W ST AT r I )= 1 

*** IF TFE FIRING RANGE TO A WAVE IS LESS THAN 75 METEPS, 

THE WAVE IS CC NS I C E R EG TC HAV* REACHED A COVERED AND 
CCNCEALSC POSITION CN THE 3EACH 

I F C (WVRNG.GT.TENCMX) ..CP. CCS UR VC I ).LT.0.05> .CR. 

*( WVPNG. LT.75. ) .CR. ( JT. GE . 2) i GC TO 5C 
JT= JT + 1 
TENGC JT >=I 

TWTS(JT) =0 EP W 1 S ( J‘* )*C SUR V ( I) 

7SLM*TSUM+TWTS< JT ) 

TRNG ( J T ) -W VR N C- 

50 IF (( WVPNG.GT.SENGMX) . CR. ( CSURVCI ).LT.0.05).CP. 

♦(WVRNG.LT.SENGMN).CR.(JS. CE.2) I GC TO 100 
J S T J S+ 1 
SENG ( JS i = I 
SPNG ( JS ) = W VR ^ G 
SWTS ( JS ) - DEFh TS( J$)*CSUPV(I) 

SSU M =S SUM+SWT S( J S ) 

IOC CCNTINUE 

*♦* CETEPMINE WAVE STATUS 

CC 20 1=1,2 
CC 25 J=1 ,5 

I F ( IWS7AT( jl.NS.l. AND. SENG (I). EC.J) IWSTAT (JI=2 
25 CONTINUE 

2 C CCNTINUE 
CC 20 1=1,2 
CO 35 J = 1 » 5 

IFCIWSTAT(J) .SC. 1) GO TC 35 

I F ( IWSTAT ( J I. EG. 2. AND. TENGC I).EO. J) IW$tat< j ) =4 
IF { 1WSTATC J i .NE.2 .AND .TE NG( I ) . EG. J > IWS T AT(J)=3 

25 CONTINUE 

30 CCNTINUE 

I F (TENG ( D.EQ.C) GC TC 500 
CO 2CC 1=1,2 

TV. IS (I ) =TWTS ( I J/TSU M 
CCNTIME 

IFCSENGC D.rO.O) RETURN 

oo oco :=i ,2 

SWlSCl I -SWTSC I ) /SSUM 
CONT INUE 
RETURN 
ENC 



200 

500 



600 



SL8P0UT INS DATAIN 

*♦* SLeROUTINE OATAIN RE*DS ALL USER SUPPLIED INFORMATION 
R ECU I R ED 6Y THE SHI P-TC-S FOR E PHASE OF TH£ MODEL 

COMMON /AMPH/IL( 5 ) ,WB< 2) , A ( 2 ) , B ( 2 ) , I TE , I SE , RD, WV INT ( 5 I , W I D, 
*TBIf,0INI7(2) ,GAINL, IWSTAT (5) 

COMMON / E NGR / S P CM AX , S PC M IN , H T MAX , HT M IN, TTS , TA A , TB , TFF 
COMMON /CISPER/TSIGV<6,2) ,TSIGH(6,2l , TME ANH C6 , 2 I , 

* S S IGV ( 7 , 2 ) ,SSIGH(7,2) 
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c 

c 

c 

c 

c 

c 

c 



c 

c 

c 



CCMMGN /CEF/7ENGMX ,SSNGMX ,SENGYN« TAR TM,S ART 
*SVEL,DEFVTS( 2) 

CCMMON /SLPEFT/GAMVA,DELTA 
CCMMCN / ICUT/TSUFV , I FR INT 
R E AD( 5, 5 C C) I PRINT 

SPCMAX,SPOMIN ,HT.MAX,HTMIN,WID 
TTS 

TENGMX, SENGMX ,SENGMN 
TAFT*,SA=TM,7'vEL , SVEL 
( ( 1 S I G V ( ItJ) ,1 = 1,6) ,J=i , 2) 



‘ T I V 5_ , 



5CC 
5 1 C 
520 
53 C 
540 
550 
56C 
570 



R6AC(5,51C) 
P E AO ( 5 ,5 10 
RcAC( 5,520 
READ<5 ,510 
REAO( 5, 530 
P E A C ( 5 , 530 
P E AC ( 5 ,520 
READ< 5, 320 
REAC(5,520) 
READ<5, 54C ) 
PEAC(5, 520 
FEAC(5, 54C ) 
RE A0( 5,320 
P CAD ( 5 , 5 6 C» 
RE A0(5 ,560 
READ< 5, 310 
REAC(5,57C) 



( nsiGvt 

( ( T $ IGF ( I,J) , 1=1*6) ,J = 1, 2) 

( (TM 6 ANF (I , J ) ,1 = 1 ,6 ) , J= 1 ,2 ) 

/ / C C T I T n * _ 1 t -I » 



< ( 5 5 1 G V ( I * J! , 1 = 1 , 7 ) ,J=1 , 2) 

( ( S3 IG H ( U J ) , 1=1, 7 ) ,J = 1, 2) 
(DfcF WT S ( I) ,1=1,2) 

( U V I N T ( I ) , 1= 1 ,5) 

(C IN IT ( I ) , I=i ,2) 

( A ( ! ) , 1 = 1 , 2 ) 

( e ( i ? ,i=i,2) 

(*eu ),i-i , 2 ) 

GA INL 

..w GAMMA , CELT A 

FCRMA7(I9X,I 1 ) 

FORMAT ( 10,F6. 2,9X,F4. 1,8>,F5.1,7X,F5.1, 
FCPMAT ( i6X,F6. 1,16X, *8.1 , 16X, F£.l ) 

F C RMA T ( 7 > , IF 10.3) 

FC FMAT ( 3 ; > ,F d* ltl4XtFS.ll 
FCPMAT (27X,5F1C.5 ) 
FCFMAT(46*,F9.5,iCX,F9.5 ) 

F C FM AT (17X»F7.2,12X,F7*1 ) 

RETURN 

ENC 



7X,F6.3 ) 



SUERCUT IN E OUTPUT 

*** sue ROUTUE OUTFIT RPCVIOES AN INPUT SUMMARY PP IMTOU"* RASED UPON 

THE OATA REC5IVEC BY S 13 R CUT I NE DATA IN. A PRINTOUT CF DISPERSION 
CATA GENERATED AS A RESULT OF CATA SUPPLIED IS ALSO PROVIDED 

COMMON /AMPH/I L( 5 ) ,W3( 2) , M 2) ,B( 2) ,1 TE ,J SE ,RQ , W V I NT 15) ,WI 0 , 

* TEW, DIN IT (2) tC-A I M , InSTAT (5) 

COMMON /C I SP ER /T S I GV(c, 2 ) ,T SIGH (6, 2 I , 7ME ANH ( 6 <• 2 ) , 

* S S I GV { 7 ,2 ) ,SS IGH( 7,2 ) 

COMMON / E NGR / $ P CMAX ,S PO M IN , HT M AX ,HT M IN, TTS , TA A , TB , T FF 
COMMON /CEF/ TE NGM X , SONGM X ,SE NC-MN, TAR TM r $ A R. TM * T V E L , 

* S V EL ,DE FU7S( 2 ) 

COMMON /2UPEFT/GAMMA ,DELT A 

♦♦♦4* INPUT SUMMARY PRINTOUT ***** 



UP ITS ( 6, 6 1C) 
UFITE(6,612) 
U R I TE ( 6,612) 
UR ITE(6,614) 
UPITEI6 ,611) 
U R ITE ( 6,615) 
UFITE( 6,616) 
UFITEI 6,617) 
URITEt 6, 6 IE) 
UFITE (6,619) 
UR ITE ( 6,620 
UPITE(6, 621) 
WRITE (6, 622) 
URITE( 6,6 22) 
U F ITE ( 6 ,624) 
U B I Tc ( 6,623) 
W R IT E ( 6,6 26 ) 
URITE(6,627) 
U R I TE ( 6,628) 
WR ITE (6,629) 
U R I TE ( 6 ,630 



( U»V I NT ( I) , I = I, 5) 
( C I NI T ( I) ,1 =1,2) 



S FOM AX , SPCM I N, HT MAX ? HT M IN, NIC 



T ENGMX , TA RT M , T V E L 
SENGMX ,SENGMN,S ARTM, SVEL 
OEFWTS( 1) ,DE f WTS (2) 

A ( 1 ) , 6 ( 1 ) 

A (2 ) ,E (2 ) 

Wfi( 1 ) ,W8( 2) 
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c 

c 

c 



W R I TE ( 6,621) G 4 IN l 
WR ITE < 6, 622) GAMMA , CELTA 



***** OISPERSICN CATA FRINTOtT ***** 



C 



I D I S P = 1 

IF (IDISF.EQ.C) RETURN 
WR I T E ( 6 ,6 33 ) 

WRITE (6 ,624 ) 

CQ 55 1=1,6 

WRIT E (6,63 5 ) TS IGV < 1 , 1 ), T S I GV( I, 2 ),TSIGH( 1,1) ,TSIGH(I,2) , 
*T MEANH ( 1,1), TM EANH < I ,2 ) 

55 COM INU E 
WRITE( 6,626) 

CO 56 1=1,7 

WRITE(6,637 ) SS IGV( 1 , 1 ) , SS IGV( I ,2 ) , SSIGH( I , 1 ) ,SSIGH(I,2) 

56 COMINUE 

WF I TE ( 6,62 6) 

WR I i e( 6 , 629) 



610 F ORM AT ( * ] • , 1 X , INITIAL SH I F-TQ-SH CRE PHASE INFORMATION >*») 

612 FORMAT (/6X,* INITIAL FORCE STRENGTH*) 

613 FORM AT (1), 'WAVE* ,ZX,*1*,5X**2 , ,5X,'3 1 , 5X , * 4 * , 5 X , *5*) 

61 A FORMAT ( 1 ) , *L VA * , f ( 2 X , F A. 1 J ) 

611 FORMAT (/ IX,* CEF. TANK ASSET S = ' , IX, F * . 1, 5X, 

* * DEF ♦ ATCM ASSETS = , ,1X,F^.1) 

615 FOFMAT( / nx, »LVA ENGR SPECS') 

616 FORMAT ( 1 > , * S PC M AX * , 2 X , ' S P CM I N * ,3X , * H T M AX * , 2X , * HT M IN • , 3 X, * W I D » ) 

617 FORMA T ( 2 > , F 5 . 2 , 3 X , F 5 . 2 , 3 X iF5*2,3X,F4«2,2X,F4*2) 

618 FORMAT ( /IX,* DEFENSIVE TAC1ICAI PARAMETERS') 

619 FORMAT ( 1CX,* RANGE* ,^X,*AIF-RELCA0',3X, 'PROJECTILE* ) 

6 2 C FORMAT ( 6 > , *M A X » , 2 X , * v I N ' , ' 7 IME' ,7X,* VELOCITY* ) 

621 FORMAT (1X,',TANK*,IX,F6.1, c X , F B . 2 , 7X, F6.2 ) 

622 F0FMAT(1>,*ATGM* ,1X,F6.1 , ]X,Fc.i,2X, F5.2,7X,F6.2 ) 

623 FORMAT (/,1X, 'DEFENSIVE x A CTI C A L ALLOCATION WEIGHTS:', 
*/,lX,*WAVE 1 = ' , F5,2,lX, *wAVE 2 = *,F5.2) 

624 FORMAT ( /lx,' DEFENSIVE FCRCE ATTRITION COEFFICIENTS 1 ) 

62 5 FORMAT (9 X, » ALP FA* a *,:OX,'EETA?A') 

626 FORMAT ( 1 >, *0T' ,6X , R7 * 5 »9X ,F7 • 5 ) 

627 FORMAT ( I X , *D S 1 , 6X , F 7. 5 , 9 X ,F 7. 5 ) 

628 FORMATt/,* AIMEC FIRE ATTRITION RATE COEFFICIENTS FOR • , 

*/,' DEFENSIVE TANK 4N0 ATGM ASSETS*) 

629 FORMAT ( /IX,' W8 ETA ( i) = ' ,F7.5, 1 x , * a BET A( 2) =» ,F7. 5) 

630 FORMAT (/IX,* BREAKPOINT ASSUMPTION: 0 • 2 A ( TOTAL DEF FORCE)*) 

631 FOFMAT ( /»1X, *OEFENCER ATTRITION LEVEL ALLOWING FCP LANC COMPAT* 
*/,lX,F5.2, 'MTCTAL DEFENDER FCRCE)*) 

632 FORMAT ( /IX,* AR T N SUP e ACTCR=* ,F5*i,2X, 'ERROR SUP F AC TOR = * , F 5 • 1 ) 
623 FORMAT ( / IX * ' DI SPERSI CN DATA*/) 

634 FOFMAT ( 3 X , ' R AN C E * , 2X , ' TS I CV • , 2X , * RAN GE ' , 2X, ' TS IGH • , 

*2X , •RANGE * ,2X, • T M c A N h * ) 

635 FORMAT ( lX,F7o 1 ,2X ,F5o 1 ,1X ,F7. 1 , IX ,F5 • 1 ,1 X ,F7. 1 , 1 X ,F5. I ) 

636 FOFMAT ( / 2X,' RANC-c* ,2X, *SS )GV f ,2X, 'RANGE* , 2X, 'SS IGH* ) 

627 FORMAT ( 1 X, F7 • 1 ,2X,F5,1 ,1X ,F7 #1 ,1X ,F5 •!) 

628 FORMAT (/, * CURRENT STATUS OF WAVE I VARIABLE DEF INI T IONS' , / / , 

*' C - NOT ENGAGING »,/, * 1 - LANOED*,/,* 2 - UNDER FIRE BY ATGM* 

2 - UNDER FIRE 37 TANK*,/* h - UNDER FIRE B v BOTH ATGM £ TANK 
639 FORM AT ('I', ****** THE SH I F-TO- SHO RE PHASE BEGINS *#»**•///) 
RETURN 
END 



SUBROUTINE PHI T(PANGE,W,h , I W P N , SU PF A C , PR HIT) 

C *** GIVEN THE RANGE, WIDTH ANC HEIGHT OF A TARGET, AS WELL AS, THE 
C TYRE OF WEAPON EEING USEC, SUBROUTINE PH I T COMPUTES THE 

C PROBABILITY OF A HIT 

C 



CCNMON /AMFH/ILC5) ,W6(2) , 4 ( 2 ) , E ( 2 ) , I TE , ISE , RD, WV INT ( 5 ) , W ID , 
*Te^,OIN IT( 2) ,GAINL , IWSTAT (5) 

CCMMCN /CISPER/TS IGV (6, 2 ) ,TS I GH( 6,2 ) , TME ANH( 6, 2) , 

*SSIGV( 7,1) ,SSIGH( 7 ,2 ) 

CC M MON /SLPE FT/C-AMMA , DELTA 



C 

C *** IWRN CODE : TANK 



ATGM =2 



no 



r»nnn nnn noon non nnnnnnnn 



4**** VARIABLE DEFINITIONS ***** 

TSIGH - THE STC CEV ERROR IN THE HORIZONTAL FOR TANK 
TSIC-V - THE STC CEV ERROR IN THE VER T I CA L FOR T ANK 
TMEANH - THE E IAS ERROR IN THE HORIZONTAL FOR TANK 
TNEANV - THE eiAS ERROR IN THE VERTICAL FOR TANK 
SS IGV/SS IGH - SIMILAR INTERPRETATIONS FOR THE ATGM 

P!=ARCGS (-1.0 
I F (RAN GE . LT. 25 . ) STOP 

I F ( I WFN.EC. 1 ) CO TO 50 
*** ATOM FIRING DATA COMPUTATIONS 
ViMEANH=0 .0 
WMEANV=C.C 

CALL INTFP(SSIGV,RANGE,WSIGV,7) 

CALL INT FFISSIGV, RANGE, WS IGH. 7 ) 

*** TANK FIRING DATA CONFUTATIONS 
50 RNEANV=C.C 

CAll INTFPii MEA Nr, RANGE* «i*EANF»6) 

CALL INTFPITSIGV.PANGE ,WSIGV,6i 
CALL INTRP(TSIGH, RANGE ,WS IGH , 6 ) 

*** CONVERSION TO MLS 
100 Z = ARSIN IF/RANGE) 

WS I GV= w S IGVM 1 ,+CELT4-*SUP FAC) 

WSIGH=WS IGF* ( 1 , + CELTA*SUP FAC ) 

T G TH= ( Z*64CO.C)/(2.C*PI) 

TC-TW=( ARS IN l A/ RANGE) )? <64 CO . 0 / ( 2.0*P I ) ) 

*** INSTITUTE NORMALITY ASSUMPTIONS TO COMPUTE HORIZONTAL 
AND VERTICAL HIT FRO E A e I L IT I E S 

C = -l. 0 ■* S CRT ( 1./2. ) 

HCR1=( < TGTW/2. I-RNEANHJ/W SIGH 

H0P2=< ( (-l.O^TG TV,)/ 2.01- *NEANF)/W SIGH 

F HITX= 1 . C 

IFIABS(HCRI) . GT. £ . ) GO TO 60 

PHITX-C.5M ERFCIC-'HORl )-ER FC ( C *H0 R 2 ) ) 

60 VE R1 = < < TGTH/2. I-WMEANV l/WSIGV 

VER2 = ( ( l-1.0*TGTH)/2.l-*M£ANV)/WSIGV 
PH ITY= 1 .C 

IF(ABS( VtPl) .GT.S.) GO TC 70 

PHI T Y = C. 5 * ( ERFC(C*VER 1 )-E R FC ( C :r VS R2 ) ) 

70 FRHIT = PHITX*PH ITY 
RETLRN 
ENC 



SUBROUTINE INTFP(X,ARG,VAL,N ) 

*** S L6RCUT I N E INTRP IS A CHECK SUBROUTINE TO INSURE THAT RANGE OF 
TARGET ANC DISPERSION DATA ARE CO w PA TABLE FOR PROBABILITY OF 
HIT COMPUTATION IN SUBROUTINE PHIT 

CINENSIGN X(N,2) 
hFITE(6,6C0) AFC 
IF< ARG. LT. X( 1 , 1 ) ) GO TO 3 0 
CO 10 1=1 f N 

IF(ARG.GT.X(I+1.1) ) GO TC 10 
D 1 F F=X ( I + 1, 1 )-X( I, 1 I 
DE IT A= ARG-X ( I. 1) 

VA l =X< I ,2) + (CELT A/C IFF)* (X( 1 + 1 ,2 )-X( 1,2) I 
RETURN 

10 CONTINUE 

I F (ARG . G T . X( N , 1 ) ) GO TC 2C 
VAL = X ( N ,2 ) 

RETURN 

20 WR ITE (6,010 
STCP 

30 ViF ITE( 6 ,620 

6CC FORMAT ( 1 > , 'ARG*****= ‘ , FI 0.3 1 

610 FORMAT! • ERROR IN INTRP A FG . GT . XI N, 2 ) ’ ) 

62 C FCFMATI* ERROR IN INTRP A FG. LT .X < 1, 1 ) • ) 
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$ TCP 
ENC 



SUBROUTINE RATE (RANGE, SPEED, IWPN,SUPRAC, ROP ) 

♦ ** GIVEN THE RANGE AN C SPEED OF A TARGET ALONG WITH THE TYPE OF 
WEAPON BEING USEC TC FIRE UPON THE TARGET AND THE SUPPRESSION 
FACTOR THE FIRER IS BEING SUBJECTED TG, SUBROUTINE RATE COMPUTES 
THE PATE OF FIRE USEC AGAINST A PARTICULAR TARGET • 

COMMON /CEF/7ENGMX,SENGMX ,S EN G MN , TAR 7M , S APTM ,T V EL , S V EL 
COMMON /SUPEFT /GAMMA f DELTA 
R CF=0* 0 

I F (RANG E .LT. 25 • ) PETLRN 
I F ( IWFN.EQ.2) GO TC 10 

IF ( PANGE*GT .TENGMX ) RETURN 

TPTM=TAFTM* (1.0+CAMN A*SUPFAC> 

CT=TPTM4RANGE/( TVEL^ SPEED) 

PCF =1 .0/ CT 
RETURN 

10 IF(RANGE.GT.SENCMX) RETURN 

I F ( RANGE . L T . S ENGMN’) RETLPN 

SRTM=$AF7M* ( 1.0 + GAMMAS SUPFA C ) 

DT=SRTM4RANGE/(S VEL^SPEED) 

RCF = 1 • 0 /C T 

RETURN 

ENC 



*** IN THE FUNCTIONS H T,SPD, 1110 PNG, THE ARGUMENT T 
IS THE TIME SINCE THE WAVE BEING ADDRESSED 
CROSSED THE 50C0 METER OFFSHCRE MARK 

FUNCTION SPD(T) 

CCMMCN /ENGR/ SPCMAX,SPDM IN , HT MAX , HT M IN, TTS , T A A , T B , T FF 
I F ( T *GT * T A A ) GC TO 50 
SFD = S FCMAX 
RETURN 

5C I F ( T. GT . T e ) GC TC ICC 

SPD=S FCVI N+ ( (T6-T J/TTS )*< SFCMAX-S PCMIN ) 

RETURN 

IOC' $FC = SPOM IN 
RETURN 
ENC 

FUNCTION HT( T) 

CCMMCN /ENGR/ SPCMAX , SPDM IN, HTMAX ,HT MIN, TTS , TAA , T3 ,TFF 
IF(T.GT.TAA) GC TC 50 
HT = H 7 F A X 
RETURN 

50 I F (T.GT . TE ) GC TC 100 

HT=H 7 M IN+ ( (TB-T J/TTS)* (HT MAX-HTMI N ) 

RETURN 

ICC HlsHTMIN 
RETLPN 
ENC 

FUNCTION PNG ( T ) 

CCMMCN / A N P H / IL ( 5 ) ,WB( 2) , A ( 2 ) , B ( 2 ) , I TE , I SE , RD , WV I NT< 5 ) ,WID, 
*TEV,DINIT(2) ,GAINL,IWSTAT (5) 

COMMCN /ENGR/ SPCMAX , SPDM IN, HTMAX ,HT MIN, TTS , TAA , TB , TFF 
I F (T.GT .TAA) GC TC 5C 
RNG = 5 COO. C - ( SPCMAX ^T) 

RETURN 

50 IF(T.GT.TB) GO TO 100 

PNG* PC- C. 5*(T-7AA)*(SPCMAX+SPD(T) ) 

RETURN 

100 FNG=PD- ( ( [7B-TAA)/2.0)*(S FDM I N+S PDMA X ) )-( (T-TB )*SPOM IN ) 

I F ( RNG • L T • 75 • ) RNG=0.0 
RETURN 
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ENC 

C 

c 

c 

c 

SIBROUT 1 N E GROUND < G ATM,T $ LRV ,IPRI NT , TTS) 

C *** THIS IS THE PRIMARY SUBROUTINE CF the LAND COMBAT PHASE OF 
C THE AMPHIBIOUS OPERATION. INFORMATION REQUIRED FOR THE OPERATION 

C CF THE L AN C CC MB AT PHASE IS PSAQ IN AND PRINTED IN A SUMMARY 

C TABLE FCF USER REVIEW, T V z INFORMATION PROVIDED BY ALL OTHER 

C SUBROUTINES USED IN T ME LAND COMBAT PHASE ARE USED IN THIS 

C SUBROUTINE AS INPUT TO THE BASIC LAN C COMBAT A LG OP. I TH YM 

REAL*8 CSEEQ 

PEAL TSLrV*TTS ,P (5 ) 

COMMON /GRPI / IPACIR (6 ) , I <ECWC(6) ,MV TO I R ( 6 ) , X < 6 ) T Y(6),$PD<6) 
COMMON /GRP2/ TA ( 2 ) , T >. ( 2 ) ,TH ( 2 ) , TM< 2 ) , TF 1 ( 2 ) t TP2 ( 2 ) , TF 3 < 2 ) , 
*P<2,6) , PHH(2 f 6 ) , PHM<2 ,6) , FKH<2,6) ,TF (2) 

COMMON / C* R P 3 / NdUvNRU* FL l t) , PC (o) ,Nu I i 3 i , X IC ( 3 ,200 i .YIC(3 ,200 ) , 
*ICIR(3,2CC),AVSP, ISPC 

*, IUSTAT( 6) ,1 I (6) , LOST '6, 6 ),VISFRA,VI$FFB,SIZETK, 

*$IZETW f NT«6) t.\F< 6>t SRFtDI SMAX, 

*NLCSC(6 «6 ) ,V 1SFR (6,6 ) , RM INTK, RMXTKt R M I NT W, RMXTW , OP , TOWFR , L VAF R , 
*PTT< 3,2 ) , P F , P 0 A ( t » 6 ) , A FOA 16 ,6 ) ,LOA(o ,6 ) , N A < 6 ) , OF L ( 6 ) , POL ( 6 ) 

COMMON /0PP4/ TPCL(6) , OLD C ( 6 , 6 ) ,Q ( 6, 6) 

COMMON / C R P5 / LOT (6 ,6) ,RCT<6,6 ) 

COMMGN /HILLS/ X C ( 1 C 0 > *YC ( 1 0 C ) • PE AK ( 100) , SX ( 10 0 ) , SY ( 100 ) , RHO ( 100 ) 
COMMON /HIlLS/ SCALE(IOO) ,TWCPHC( 100 ),TWOSCL( 100) , BASE 
COMMON /HILLS/ NHlLLS 

COMMON /COVER/ CXC( 1531 tC 15C) t CPE AK(1 50) tCPXX (150) tCPYY ( 150) 

COMMON /COVE F / C 5 XY( 150), NOVELS 

COMMON / CQLN TR / K H , KH * , KV ,KN , KGRS , KE LL ,K I NT 

COMMON /GRID/ LS T ( 10 , 1 0 ) ,MHL( 1C • 1 0) , L I STH ( 45 0) ,KHREP(100) , KTREP 
COMMON /GRID/ LSTC ( 10 , 10 ) ,NC ( 10, 1 0) , L 1ST C ( 400 ) , KCR EP ( 150 ) 

COMMON /GPPfe/ A L P H A ( £ j 

COMMON / C P P7 / ,\AU),VA<6) ,lMCVEf6) 

C*** IMTIALI2ATI ON. 

C *+* VARIABLE CEFINI T ICN$ ***** 

C BL - INITIAL DEFENSIVE P CFCE LEVEL 

C RL - INITIAL AGCRESSCP FCFCE LEVEL 

C Z L - ZERO LEVEL 

E l =0 • 0 
P L = C . C 
M P = 0 

RA1=3. 14159 
ZL=.00CC1 



c 

C*** 

c*** 

c 

c 

c 

c 

c 

c 



c*** 

c 

c 

c 



READ TERRAIN DA*A FOR LINE OF SIGHT 
CHECK FCR STOCHASTIC OR DETERMINISTIC ATTRITION 
I T R IT-ATTR I T ICN MODE 1=0 ET5RMIN I S TI C 

0 = ST CCH AS I TIC 
DSEFC-00U6LE PRECISICN SEED NU M 6 ER 

PP ANO CC ARE THE BETA DISTRIBUTION PARAMETERS FOR DBF UNITS 
PD AND CD APE THE BETA DISTRIBUTION PARAMETERS FOP ATK UNITS 

F E AC ( 9 , 5CC ) ITRIT ,CSEEC 
R E AD ( 9 , 5 C 1 ) PP,OC 
CC 5 1=1 ,2 

CALL GGBTR ( CSEED ,PP, Q C, 1 , R ) 

ALPHA (I > = P ( 1 ) 

CONTINUE 

R EAD( 9 » 5G0 ) 1TRIT,CSEEC 
READ<9,501) PC,CC 
DO 6 1=4 ,6 

CALL GGBTR ( CSE ED ,P0, Q 0, 1 , R ) 

ALPHA ( I )= R ( 1 ) 

CONTINUE 

READ IN NUMBER CF ATTACK ANO DEFENSE UMTS 
NBU - NUN 6 E R CF DEFEf'SIVE UNITS 
NR U - NUMBER CF AGGRESSOR UNITS 
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c 

o** 

c 

c 

c 



c 

c*** 

c*** 

C*** 

c 

c 

c**+ 

c 

c 

c 

c 

c 

c 

c 

c 



READ(9,5C2) NBt,NRU 



INITIAL I2E WEAPON SI ZES 

SIZETK - SIZE C»= IV A WEAPON SYSTEM 
SIZETW - SIZE OF TOW WEAPON SYS TEN 



5IZETK=2.5 
S I ZETW = 2*5 



READ IN EFFECTIVE WEAPON PANGES 

P Y IN T K A N C R M X T K ARE MAX AND MIN RANGES OF LVA MOUNTED WEAPON 
RMINTW A N C RMX7W FOR m>:x AND MIN RANGES OF TOW DEFENSIVE WEAPON 



R E AC ( 9 » 5C2 ) R Y INTK ♦ RMXTK , FMI NTW »R YXT W 



INITIALIZE PM,RF,7CWFR,LV/FR AND NOD 

PY - PFCPORTICN OF TIME A MOVING UNIT IS SEARCHING c OR TARGETS 

RF - CE7ECTICN PATE RECLCTICN FACTOR FOR A FIRING UNIT 

(IN CCYPARISCM TO t NCNFIRING UNIT) 

TOWF R - FIRING RATE DEFEND INC TCW WEAPON SYSTEM 

LVAFS - FIRING RATE ATTACKING LVA WEAPON SYSTEM 

NOD - NUMBER OF TIME INTERVALS UNI T I DELAYED IN MOVEMENT 
(TOO FAR IN FRONT OF OTHER JNITS) 



Py=. 352 
R F = . 5 
T 0 WFR = • C 3 
L V AFR = # I 
NOD = 2 

DO 10 I - I » NR U 
NOI ( I ) = 12 5 
10 CONTINUE 
K= NRU+ I 
L = NRL‘+N El 
00 15 I = ] ,L 
II( I ) = C 
15 CONTINUE 
C 

C*** READ IN FCRCE LEVELS CF EACH AGGRESSOR UNIT 
C 

I S UR V= INT (T SUPV/NRU) 

CC 20 1 = 1, NRU 

FL ( I )=FLOAT (ISUFVi 
20 CONTINUE 
C 

C *** CHECK FOR TYPE CF RCUTE DE TERM I NITICN 
C 

RE AD ( 9 , 5 04 ) IRTE,ISPC 
***** VARIABLE DEFINITIONS 



* 4z± ** 

IRTE - DENOTES WHETHER USER WANTS TO INPUT ROUTES OR NOT 

0 - FFOGRAy DETERMINED ROUTES 

1 - LSER DETERMINED ROUTES 

INPUT VARIABLE TC DENOTE USER'S DESIRED SPEED FOR 

1 - 9 MPH 

2 - } 2 MPH 

3-15 MPH 
4 - ]S MPH 

AGGRESSOR FORCE MOVEMENTS 

AVERAGE SPEED CF AGGRESSOR FORCE MOVEMENTS 
DISTANCE IN METER S TO BE MCVEO EACH TIME STEP BY 
AN AGGRESSOR UNIT 



ISPO - 



A VS C 
DST 



C 

o** 



IF< ISPD.EC.l ) 
IF (IS PD. EC. 1 ) 
IF (ISPO. EC. 2 ) 
I F ( ISPD.EC.2 ) 
IF (ISPD.EC.3 ) 
I F ( ISPD.EC.3) 
IFUSPD.EC.4 ) 
IF (ISPO. EC. 4) 



AVSP=S.O 

CST = 4 0 .232 

AVSF=12.0 

CST=53.6A3 

AVSF=15.0 

CST=B7.C53 

;vsp=is.o 

C S T = 80. 463 



R E A 0 IN INITIAL AGGRESSOR LNIT'S LOCATIONS 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 



30 

250 

255 



CO 25 I»],NRU 

PEAC(S,506) XICU,1),Y!C(I,1) 
CONTINUE 

IF(IRTE.EG.l) GO TC 250 



DO 



30 

DC 



1=1 ,NRU 
3 C J= 2 , 1 2 5 

YIC < I ,J)=YIC( I, J -II +CST*< J- 1) 

xicn ,j)=xic< it j-i)^cst*(j-i) 

IDI R( I »JI =0 
CONT INUE 
GO TO 255 
CALL ROUTE 
SL MR0 = 0 • C 

***** STATE VARIABLE CEF1NITICNS ***** 

FL<n - FCPCE LEVEL CF UNIT I 
SUMPC - TOTAL /C-GFESSGP FCRCE LEVEL 
MVTDIPCI) - MOVEMENT DIRECTION CF UN I T I 
I DIR ( I * J ) - DIRECTION CF THE JTH INTERVAL 
IUSTATU) - CURRENT STATUS UF UNIT I 

0 - UNIT ALIVE ANC NOT FIRING 

1 - UNIT ALIVE ANC FIRING 

2 - UNIT DES.T FCYEC 

3 - UNIT MOVING 

NUMBER CF TIME INTERVALS UNIT I IS ALLOWED TC FIRE 
AT SAME LOCATICN 
INTERVAL INDEX FCF UNIT I 



IN THE ITH ROUTE 



NFC I) - 
1 1 ( I ) - 



35 



CC 35 I = 1 ,NR U 
F0( I )*FL< I ) 

SUMRC = SUMPOFO(I) 

X< I ) = > ICC I , 1) 

Y( I )=Y ICC 1,1) 
MVTOIRC I) = I C I R ( I , 1 ) 
SPO( I ) = AV SP 
IUST A T ( I ) =0 
IPROIF ( I) = I Cl R < I ,1) 
IS EC W D( I ) = 1 20 
NFC I ) = 1 
11 ( 11*1 
CONTINUE 



C*** READ 
C 
C 
C 
C 

c 



IN DEFENSIVE LNIT'S LCCATICN3 

***** S T A T E VARIABLE cefiniticns ***** 
IPRCIF(I) - PRINCIPLE CIRECTIO.N OF FIRE FOR 
ISECWC(I) - WIDTH OF SEARCH SECTOF FCR UNIT 
SUM BC - TOTAL CEFENSIVE FCRCE LEVEL 



UN IT 
I 



AO 

C 

c*** 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 



$ LMBQ= Q • 0 
CC AO I = K , L 

READ ( « ,50 7) X < I) , Y C I ) , »L Cl) , IP PD I P < I) , IS ECWDC I > 

F0( I ) = F L ( I ) 

SUMBC = SUM SC + F0 ( I ) 

MVTDIR (I) =C 
SP 0 ( I ) = C • C 
IUST A T ( I ) =0 
IMOVE ( I)=C 
CONTINUE 

CHECK FCR ALT ER N A T E 0E C EN S IVE POSITIONS AND READ IN IF WANTED 
***** STATE VARIABLE DEFINITIONS ***** 

I ALT - INDICATES IF ALTERMTE CEFENSIVE POSITIONS DESIRED 

0 - YES 

1 - NO 

BREAK - CLOSEST DISTANCE ALLOWED BETWEEN OPPOSING FORCES 
EEFORE BREAKPOINT REACHED FOR PROGRAM 
ITEM - NUMBER CF TIME INTERVALS ALLOWED FOR DEFENDER'S MOVE 
TC THE ALTERNATE CEFENSIVE POSITIONS 
T A ( K ) - TIME TO ACQUIRE A TGT FCR KTH wEAPON SYSTEM TYPE (K = 1,2) 

T 1 ( K ) - TIME TO FIRE 1ST FOUND AFTER TGT ACQUIRED BY KTH WEAP SYS 

TH(K) - TIME TO FIRE A RCLND FOLLOWING A HIT FOR KTH WEAP SYS 

TMK) - TIME TO FIRE A ROUND FOLLOWING A MISS FOR KTH WEAP SYS 
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TF1(K) - TIME OF FLIGHT FCR KTH WEAP SYS PROJECTILE TO 1000 METERS 

T F 2 ( K ) - TIME OF FLIGHT FCR KTH WEAP SYS PROJECTILE TC 2000 METERS 

TF3(K) - TIME OF FLIGHT FCR KTH WEAP SYS PROJECTILE TO 3000 METERS 

REAO(9,5C8) I ALT ,EREAK ,1 TEM 

IF ( IALT . EC.l ) GO TC 260 
DO 4 5 I =K t L 

R E AC (9 , 5 06 ) XA(It 9 YA(I) 

45 CCNT I MJE 
26C CELT-10. 

TA ( 1 )= 2C • 

T 1 (i ) = 8 • 

T H < 1 ) = 8 • 

TM ( 1 )= 1 0 • 

T F 1 ( 1 ) =1 . 

TF2< 1 ) = 1# 

T F 2 ( 1 ) =1. 

T A ( 2 ) =20. 

T](2)=B. 

TH (2 ) = 8 • 

TM ( 2 ) = 1 5 • 

TF1(2)=10. 

TF2( 2 ) -1 2 • 

T F 2 ( 2 ) = 1 f . 

*** RE A C IN HIT AND KILL F R CB A E I L I T I ES 

***** STATE VARIABLE DEFINITIONS ***** 

P(ItJ) - F RO 8 1ST RC'JNC HIT 8Y UNIT I IN RANGE BAND J 
FHH(ItJ) - PRCE OF HIT FOLLOWED BY A HIT 

P FM ( I , J ) - PRCE OF HIT FOLLOWED BY A MSS 

PKH(I,J) - PRCE OF A KILL GIVE H A HIT 

PTT(IfJ) - PROPORTION SURVIVING R I RE POWER ALLOCATED TO 

ITH TARGET IF J TARGETS ARE AVAILABLE 
NLGSCUfJl - NUMBER CF CONTINUOUS TIM INTERVALS THAT A LINE OF 
SIGHT(LOS) DCES NOT EXIST BETWEEN UNIT I AND UNIT J 
0(1, J) - FPOEAEILITY UNIT J NOT DETECTED BY UNIT I AT CURRENT TIME 
VISFR(ItJ) - FRACTION CF FEIGHT OF T GT J VISIBILE TQ FIRER I 
IRAN - RANGE 

DO 55 1=1*2 
CC 5 C J = 1 * 6 

RE AC 19 ,509 ) P ( I, J ) , FHH ( I • J ) » PHM( I, J ) ,PKH ( I , J ) 

50 CONTINUE 
55 CONTINUE 

P 11(1,11*1.0 
P7T< 1, 2)=C.B 
FTT ( 2 , 2 ) =C .2 
P1H 1*3) = C. 9 
PTT12,3)=C.15 
FT T (2 ,3 ) =0.05 
DO 60 I = 1 * NR L 
DO 6 C J =K , L 

N LC SC ( I * J ) = 0 
NL C SC( J ,1) =C 
0 ( I* J) =1 .0 
Q ( v * I ) =1 .0 
VI <FR( I , J)=C.O 
V I S F R ( J » I)=C.O 
6 C CONTINUE 
IC = 1 

C*** PRINT INITIAL BATTLE INFOR NATION 
C 

WRITE(6,6C0) 

WRITE(6,6G1) 

CO 65 I = 1 * L 

WRI TE <6,6 03 ) I ,X( I) ,Y ( 1) , FL (I ) 

65 CONTINUE 

IF ( ITRIT .EC. 1 ) GC TO 265 
WR I TE ( 6 ,6 04 ) 

GO TC 270 
265 WFITE(6 ,605) 
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27C IF (IRTP. EC.O ) GQ TC 275 
WRITE(6*6C6) 

275 WPITE(6,6C7I AVSF 
WP I TE ( 6 , 6CE) SFEAK 
I F ( IAL T . EC.O ) GC 70 28C 
WRITE (6,609 ) 

GO TO 2 85 
280 WPITE(6,61C) 

CC 70 I = H , L 

WRI TE (6 ,6 11 ) I ,XA( I ) r YA< I ) 

70 CONTINUE 
2E5 I R AN=5 0 C 

WP ITE ( 6,6 12) 

CO 75 1=1,6 

WRITE (6,6 13 ) I P AN ,P (1 , ] ) , PHH(1 ,1) , P HM ( 1 , 1 ) , PK H ( 1 , I ) 

IR AN= IRAN+500 
75 CONTINUE 
IRAN=5C0 
WF IT 5 {6, 6 16 5 
OC 80 1=1,6 

WRI TE <6 ,612 ) IRAN ,P (2 , 1) ,FHH(2 ,1 ) , PHM (2 ,1) ,PKH(2 , 1) 

IRAN= I RAN + 5 CO 
8 C CONTINUE 

W P ITE ( 6,615) 

WPITE(6,617) 

WF I TE < 6,618) 

CM* UPDATE LCCATICN CF REC UNITS, 

NF - NUMBER OF TIME INTERVALS UNIT I COES NOT FIRE 
NOD - NUMBER OF TIME INTERVALS UNIT I DELAYED IN MOVEMENT 

01 SMAX = 5CCC# 0 
250 OO 50 1=1 ,-\RU 

I F( IU ST AT ( I ) . EG. 2 > GO TO 50 

IFdUSTAT(Ii.EG.O) GC TC 295 
NF( I ) =NF ( I I + 1 

IF ( N F ( I ) .LT.NCD) GC TO 90 
NF (I ) =1 

295 00 85 J = 1, NPU 

IF (J .EC. I) GO TO 85 

IF ( I US T A T ( J ) .EQ. 2) GO TO 85 
CIST = X{ I) - X( J ) 

IF (CIST .GT. OISMAX) GO TO 90 

85 CONTINUE 

II( I ) = IKII 4 1 
K7= I I (I) 

X( I ) = >IC( I , K 7 ) 

Y ( I I = V I C ( I, K7 ) 

MVTD IF ( I ) =1 Cl P ( I ,K7 ) 

IPRO IF(I)=ICIR(I,K7) 

90 C C NT I NU E 

LINE— OF-SIGHT CHECK BETWEEN UNI^S AND TARGETS SELECTION 
***** STATE VARIABLE DEFINITIONS ***** 

NT ( I ) - NUMBER OF TARGETS DETECTED 3V UNIT I 
X X 1 , YY 1 - COOR C KATES CF UNIT I LCCATICN 
X X2 , YY2 - COCPCINATES CF UNIT J LCCATICN 
TMAC I , TM AC J - ELEVATION CF UNIT I ANC UNIT J 
0,C - INCICATES NO UNITS UNDER GRCUNC 

SI2ETK, SIZETW - SIZE CF L V A VEHICLE ANC SIZE OF TOW VEHICLE 
LATCB - INDICATOR VARIABLE FCR ONE OR TWO WAY LOS CALLS 

0 - CC NOT COMPUTE LOS FROM UNIT A TO UNIT B 

1 - COMPUTE LCS F F CM UNIT A TC UNIT B 
(YIELDING VI SERB) 

VISFRA - FPACTICN CF HEIGHT CF T GT AS SEEN BY UNIT A 
LCST(ItJ) - INCICATES IF ICS EXISTS BETWEEN UNIT I AND UNIT J 
0 - NO LCS EXISTS 1 - LG S EXISTS 



CO 55 J=K,L 
NT ( J ) =0 
55 CONTINUE 

CO 105 1=1, NRU 
NT ( I ) = C 
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2F< IUSTAT <1 ).SC .2 ) GO TG 1C5 
DC 100 J=K,L 

IF( 1UST4K J).EQ. 2. OR. IUSTAT (J). EQ.3) GO TO 100 
XX1=X ( I ) 

Y Y 1 = Y ( I ) 

CALL ELEV ( XX 1 ,YY1,TMACI ) 

X X 2 = X ( J ) 

Y Y 2= Y ( J ) 

CALL ELEV<XX2 ,YY2,TMACJ > 

L A 10 e = 1 
L 8 TO A = 1 

CALL LOS (XXL ♦ YY1 t TMAC It 0 .0 , S I 2ETK , X X2 , YY 2 , TMAC J , 0 . C , 
*SIZETU, LATCB ,L0TOA » V I SFR A , VI S F F8 ) 

VI $ F R ( I tJ ) = V ISFKA 
\ISFR(JfI)=VISFR6 
IF < VISFRA.GT. 2L) GOTO 30 0 
LOST { I , J ) = 0 
LCS T { J , I) =0 

NLuSC(IfJ) =NLO SC t I , J Mi 
NLOSC (J , I) =NL 0 SC < I f J ) 

GO TG ICO 

200 l OST ( I i J ) = 1 

1.0 ST ( J t I Ml ‘ 

NLOSC ( I ,J) =0 
NLOSC ( J » I ) =0 

RANGE* SORT ( (X (I MX (J) M*2 + <Y ( I MY( J ) )**2 ) 
IF(RANGE.LT.RMINTK. OR. RANGE. GT.RMXTK) GC TO 305 
IF (0(1, J) . cQ.l .C ) GO TO 305 
I UST AT ( I) =1 
NT(I)=NT(I) + 1 
M =NT ( I ) 

LOTH , M )=J 
RCT ( I* M ) = R ANC-E 

I F ( M.6G.1 ) GO TO 30 5 
CALL S OR T ( I t M ) 

205 F (RANGE.LT. RMINTW. OR .RANGE. GT.RMXTW) GO TO 100 

IF ( 0 ( J*I). EC.l.C) GO TC 100 
I USTAT ( J) =1 
NT (J ) = N1( J Ml 
M = NT (J ) 

LCK J, M )=I 
ROT ( J » M )=R ANGE 
IF(M.EC.l) C-0 TC 1 CO 
CALL SORT ( J f M ) 



ICC CCN 7 I MJE 
I C 5 CONTINUE 



co no 



= 1 »NRU 



l F ( IUSTAT ' I ) . EC. 2 ) GO TO 110 
J F ( N T ( .[ I.NE.C) GO TO 11C 
IUS T A T ( I ) = 0 
NF ( [ ) =0 

110 CONT IN U E 

CO 115 J - K 1 1 

IF( Il>STAT(J).cC.2.0R. IISTAT(J).EQ.3) GC TO 115 
I F ( NT ( J ) • EC • 0 ) I UST AT ( J ) =0 
115 CONTINUE 

C*** UPCATE C F TFE ACCUMULATED CETECT ION PROBABILITIES. 

C 

I A A- 1 
I E E-NRU 
I C C = K 
I C C= L 
F P = TGWF P 
C F = PM 

CO 120 1*1,6 

12 C CONTINUE 
2C7 DO 125 I = I AA » I EB 

IF( IUSTAT ( I J.EC.2.0R. ILSTATU ).EQ .2) GO TO 135 
DO 150 J = I CC 1 1 CD 
P RO P * C. C 

I F( iij ST iT(J).EQ.2.0R. IUSTAT (J1.E0.3) GO TO 130 
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21 C 
215 



22C 

225 

220 
2 2 5 



125 



24 C 
245 

25C 

120 

125 



C 

C*** 

c 

c 

c 

c 



140 



145 



C L C C (I « v f 2 ~ C i I 
I F ( L C 5 7 ( 1 * J) . 
1 F (NT (I ) • G 
PC TV IS « 
CALL L A 
C V = E X ? ( 
IF (N7( J 
Q ( I , 
GC T 
CP = ( 1. c 
C ( I » J ) = 
GO TO 1 
N 5 =NT ( I ) 

CC 12 5 1 
K1 = LGT ( 
ANG i=A T 
ANGZ=A7 
I F ( ( A N G 
I F ( A 
A 
G 

ANG = 
I F ( A 



• J » 

£0. C ) 
'.0 ) 

MS-? 
‘.n * 



GC T C 350 
GC TO 212 
( J, I ) 



.0, PCTM StCE'RA" 
-( 0E7RA7*0P *0£L7*=L( J 
i. G T . C ) GG 7C 3 1G 
v ) = c ( ;,jmc v 

0 1 2 C 

-PSL5K ) * * ( r c * D E L 7 * F L ( 

C( : , „ )*<cv+cp-cv*c?) 



CCN7 
IF (P 
I 



NGa 
A =1 
F { A 

p 

I VJ 

RCr 

F(N 

C 



C(I 

GO 

I F ( I A A 
Oil 
GC 
I F ( 



0 

G 

• J) 
7 C 

z ^ 

* jT 

70 

\LC 

; u 



7(. 
ALL 
5T = 
V= : 
( I 

%] 



1 * N * 
Ml ) 
A NZ ( i 
AN2( > 
1*A.\G 
NS 2. L 
NG=2* 
0 TO 
I*PA I 
SG« G7 
0 70 
( AN G 2 
« p a :/ 

C7. AA 
= PRC = 

3 > • C J 
« ^ - 

<* - • 

4 la>: 
m - ( 

J) -c < 
0 12 C 



( <1 )-Y (I ) f X (<1 ) - 
U)-Y( I ) ,X( j)- <( 
1 ) • G2 • 0.0 GO 70 
".0.0) GO "0 320 
?a: +ANC1-ANG2 
225 

-ANG2- ANSI 

*P A ! ) ANG= 2* D AI — 

-ANSI ) 

130 .0 

) GO T C 12 5 
•* => 77( I 1, N5 ) 

GO 7C 2^*5 
G )SCTC 3<*C 

?C7V I S t DE7 
“ R A 7 * - F 

PRC?*! E T RA7*DEL7 



SC ( 



) * 



1,J >. 

. w 



0 120 

wE.3) GC 70 130 



CCN7 IMJ2 

CCMINJE 

I F ( lAA.EC.X) GC 70 225 
FR = L / ArR 
I A A = K 
! 63= L 
ICC=1 
I DO = N F U 
CP=1. C 
GC 70 2C7 

FIRE ALLOCATION. 

***** state variable de- :v: t :cn s 

AFCAdtol - AVERAGE — CPCMCN OF THE 
ALLOCATE: 70 FIRE ON UNIT 



4x xxx 

J7H aggr 



CC 1^0 1=1, L 
N‘a( I )=C 

DO 155 1*1 .L 

I F( ILS7A7C ) • EC • 2. OF. IL’STA* (:)• EC. 3) 
IF (NT ( I ) . EC. 0 1 GC "0 155 
00 145 J = 1 ,3 

A PCAC, „ ) = C • C 
CONTINUE 

IF (NT (I ). EC.l ) GO TG 270 
I F { N ~ ( I ). EC. 2) GO 70 365 
NOT =2 

M W I a LOT ( : ,1) 

MY 2 = LOT ( I ,2) 



, PSUEK ) 
) ) ) 

J) ) 



X ( I ) ) 
I ) ) 

- n r\ 



i.NG 



RAT , =>SUEK) 
*RL( J) > ) 



ESSCR Or UNIT I 



GC 70 155 
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36 C 



150 

365 



37 C 



MM3= LOT { 1,3) 

PP0E*(1.C-C( I,MM1 ) )*C(I,MM2 )*C< I, MM3) 

APOA (I , 1 ) =APG4( I f 1) +PTTM. tl )*ooCB 
PR0B=Q( I , M.mi )*< 1 .0-C( I .MM2) ) X Q ( 7 » MM3 ) 

APOA (I ,2 > =APCA( I ,2> + PT7 < 1 , 1 l*?PCE 
PROS =Q ( I , M M 1 ) *Q ( I » MM 2 ) * ( 1 • 0 - C ( I » M M3 ) ) 

A PC A (I , 3 )*APCA( 1,3) *PTT<1, 1 )*?RC'8 
PROS *(l.C-C(Mn)!M1.0-C( I r MV < ; I I *Q ( I » MM3 ) 
A PO A ( I , 1 ) - APCA ( 1 , 1) *P7T( I ,2 I >PFf_B 
APOA (I , 2 ) = APCA( 1,2) + P77 <2 ,2 )*PRC3 
PR0B=( 1 . C-Q< I ,M Ml ) ) *Q( I ,MM2 ) * i 1.0-Q( J , MM3 ) ) 
APOA (I , 1 ) = APOA( 1,1) +PTT ( 1,2) -*PRQB 
A PO A (I * 3 ) - A?CA ( 1 ,2) + P'T (2 ,2 » *?PC2 
P ROB =Q ( I ,MM 1 ) * ( 1. C- C ( I » MM 2 ) ) * ( 1 . 0-0 ( ) , MM3 ) ) 
APOA (I ,2 ) = APOA ( I, 2H-P 7Tfl, 2 I** P OB 
APOA (I ,2 ) =A?C4( 1,3I+PTT(2 i2 l*PRCe 
FPCB* (1.0-C( I , MM1) )♦( 1 . C- 0 ( I , MM 2 ) ) ^ ( i„G-Q( ) ,MM3) ) 
APOA il , 1 )=APCA( 1,1 ) +977 11 ,2 1*P»C3 
APOA (1,2) =AP CA ( I ,2) +PTT( 2 ,3 ) *-PfiC3 
APOA (1,3 ) = APOA { 1,3) +PTT(3, 3)** ROB 
00 150 J = 1 , NCT 
— KK= L OT ( I , J) 

N A ( K K ) = N A (KK )+l 
I N = N A ( KK ) 

LCA(*K, IN)=I 

PCA { I* K , I N ) = A PO A <I,J ) 

CONTINUE 
GO TO 1 55 
N C 7 = 2 

MM 1 = LOT { I , 1) 

MM2 = L OT (1,2) 

P F C a = (1 .0- C ( I ,M*1> )*Q ( I , M M 2 ) 

APCA( I, 1) = /POA ( I, 1 ) + ? 7T( 1, 1)*PR0B 
PFC8-C(I ,MM)*{ 1.0-Q( I , MM 2 ) ) 

APCA< I ,2) = APCA ( 1,2 )+PTT(l , 1 ) *PRC3 

PFCe= (1.0-U I ,MM1) )*( i. c-Q( I , MM2 ) ) 

A F C A { I , 1 ) = / PCA ( 1,1 )+PTT (1 ,2 )*PR C3 
A PC A ( I, 2) = APOA ( 1,2) + P IT (2 ,2) *PR OB 
GC TC 360 
NO 1 = 1 

MM 1 = LCT( I , 1) 

PRC£=1.0-C(I, FM1 ) 

APCAC1 ,1»=AP0A( I ,i)+PTT( 1 ,l)*PROe 
GO TC 360 



155 CONTINUE 
C 

C*** ATTFITICN COMPUTATION 

***** STATE VARIABLE DEFINITIONS ***** 

RANGE - CURRENT MINIMUM 0 IS7ANCE BETWEEN AGGRESSOR /NO DEFENDER 
POA - PFCFCRTICN CP THE J7H ATTACKER CF UNIT I ALLOCATED TO 
FIFE ON LNIT I 

T POL - TOTAL PERCENTAGE LOST SINCE START CF BATTLE PCR UNIT I 
AVD - AVERAGE CISTANCE 
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SUMR=0 .0 
SUM8=0. C 
CO 165 1*1, L 

IF ( IUSTAT ( I ) . EC. 2 .OR. ILST AT (I) .EG .3 ) GC TO 165 
M 6 = N A ( I ) 

SUM=C .C 

IF(M6. EC.O) GC TO 3 c 5 
CO 1 6 C J = 1 , M6 
M7=LCA{ I ,J) 

IF(M7.LT.K) GC TC 375 
1 1 7PE = 2 
GC TO 330 

iivpe*i 

RANGE = SQR T (( X(I ) - MM7 I ) * *2+< Y ( I ) -Y ( M7 ) )**2) 
IFUTRIT.EQ.l ) GC TO 385 

CALL STGCH ( I TY P E .RANG E , A J I ) 

GO TO 3 SO 

CALL 8 7K ( ITYP E, RANGE, T ) 
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A J I = 1 • Q/T 

SUM=$UM+AJl*Fl( M7 )*POA( I , J )*DELT 
CCNT INUE 
OF l ( I ) =F L ( I ) 

FLU )=FL(I ) - SUM 
IF (FL( I 5.GT.ZL ) GO TO 400 
FL ( I ) = C . 0 
IUS T A 7 ( I ) = 2 
IF (I.LT.K) GO TO 405 
SLM B = SUM B + F l ( I ) 

TPQL( I ) = ( FO (I )-FL( I )) /FO ( I ) 

GO TO 165 
S L F P =■ S U M P + r L ( I ) 

TPCLCI ) = ( FO ( I } -FL ( I ) ) / FC ( I ) 

165 CONTINUE 
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C*** PRINT ANO CHECK FCF BATTLE TERMINATION. 

C 

c 

III ME = I C < I NT (TTS) 

DC 175 I 3 K t L 

I F ( IUSTAT (I ). 6C.2 ) GO 10 175 
DO 170 J = 1 , N RU 

IF( IUSTATt J). EQ. 2) GO TO 170 
CHECK=X( I ) -X ( J ) 

A V C-SCPT ( ( X ( I )-X ( J) )**2+ (Y (I )-Y(J ) )**2 ) 

IF(AVD. L7.BRE4K.OP.CHECK.lt. 50.) GO TC410 
170 CONTINUE 

175 CONTINUE 
C-0 TO 415 
C 

C*** COMPLETE AGGRESSOR UNIT'S FGVE 
C 

410 CO 180 I = K ,L 

IF( IA17.EC. 1.0R.IMCVE( I). EC. ITEM) GO TO 440 
IF ( IUS 7 AT ( I ) . EC. 0 ) 1UST A T( I ) -3 
I MOVE ( I ) = I M CV 5 ( I ) + l 
I F ( I MOVE ( I ) . L T. I ~EM ) GO TO 180 
X( I ) = X A ( I) 

Y ( I ) = Y A ( I ) 

I F ( IUSTAT ( I ) . EG. 2 ) IUSTA7(I)=0 
180 CONTINUE 

415 IITIME=ITIME+IFI X(GA7M) 

I F ( IPP IN T. EC. 1) GO TO 430 
WP IT E (6 ,622 ) 1 1 T I N 5 

WRITE (6,623) 

WR ITE (6 i 624 ) 

00 185 1 =1 , NRU 
M6 = NT ( I ) 

IF (N6 .NE.O ) GO TO 420 

WRITE (6 ,625 ) I , X ( I ) , Y ( I ) , FL ( I ) , IUSTAT ( 15 » TPQL ( i j 
GO TC 1 1 5 

42C VRITE (6,6 26 ) I ,X( I ), Y ( I ), FL( I ) , IUSTAT( I ) ,TPOL( I ) , 

*(LOT(I,J ) , J= 1 , N 6 ) 

185 CCNTINUE 

W F I T E (6 ,627) 

WPITE(6 ,628) 

NNN = N PU+ 1 
DO ISC I =NNN , L 
N 6 = N T ( I ) 

IF(N6.NE.C) GO TO 42 5 

WR I TE ( 6 ,625 ) I ,X( I ) ,Y( I ) ,FL( I ) , IUSTAT ( I ) , TPCL ( If 
GO TO ISO 

42 5 WR I TE(6 ,6 30 ) I ,X( I ) , Y ( I ), FL( I ), IUSTAT ( I ) , TPOU I ) t 

* ( LOT ( I ,J ) , J= 1 , N6 ) 

190 CONTINUE 

4 3 C CONTINUE 
C 

C*** CHECK FCR EATTLE TERM IN AT I CN. 

C 

I C T=0 

C*+* CHECK IF AN AGGRESSOR FQRC E UNIT IS STILL ALIVE 
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00 195 1=1 ,NRU 

IF(FLd).EQ.C.C) GO TO 195 
IOT * 1 
155 CONTINUE 

IF ( IOT. £0.1) GO TO 5 25 

C*** NO AGGRESSOR FCFCE UNIT ALIVE SC OEFENOER WINS 
WRI TE (6,621 ) 

GO TO 550 
535 I CT = 0 

C**» CHECK IF A DEFENSIVE FORCE UNIT IS STILL ALIVE- 
CC 200 I = K , L 

IF( FL ( I ). EC.0.0 i C-0 TC 200 
IOT= 1 
200 CONTINUE 

IF(I0T.5C.l) GO TO 555 

C*** NO DEFENSIVE FORCE UNIT ALIVE SO AGGRESSOR WINS 
WRITE (6,622 ) 

GO TO 550 

C*** C I $ T ANC E EETWE5N FORCES TOO CLOSE 
55 C WRITE (6,632) 

GO TC 550 
555 IC=IC+1 
C-0 TO 250 
5 5 C CONTINUE 

WRITE(6,622) IITIME 
V P IT E ( 6 .6231 
WR ITE ( 6,625) 

CO 205 1 = 1, NPU 
N6 = N T ( I ) 

IF(N6.NE.C) GO TO 555 

WR ITE (6,625 ) I ,X (I ) ,Y ( I ) , F L (I ) , lUSTAT(I) , T POL ( I ) 

GO TC 205 

555 WRITE (6,6 26 )I,X( I ), Y( I ),FL (I ), IUSTAT( I) ,TPOL( I ) , 

* ( l OT ( I , J),J=1,N6) 

2C5 CONTINUE 

WRITE<6 ,627) 

WRITE(6,62EI 
NNN=NRU+ 1 
CC 210 I = NNN , L 
N6=NT (I ) 

IF ( N6 ,NE. C ) GC TC 560 

WRI T E ( 6,625 ) I , X ( I ) , Y ( I ) , F L ( I ) , IUSTATd ) ,TPOL( I ) 

GO TO 21 C 

560 WRITE (6,6 20 ) I,X( I ), Y< I ), FL ( I) , IUS TAT( I ) ,TPOL( I ) , 

* ( LOT ( I , J) ,J=1,N6) 

210 CONTINUE 

500 FORMAT ( 1 8X , I 2 , 12X , F1C. 2 ) 

501 FORMAT(25X,F12.2,7X,F12.2 ) 

502 FORMAT ( 1 1 X , I 2 » 16X , 12) 

503 FORMAT (1CX ,F7.1 ,10X , F7 .1 , 10X , F7.1 ,10X, F7 .1 ) 

5C5 FORMAT< If X, I 2 , 1SX ,1 2) 

506 FORMAT ( 12X.F3.1 , 12X, FS.l ) 

SC7 FORMAT ( 1 2 X ,r 3 . 1 , 1 2X , F3 .1 , 6X , F5 .1 , 7X , 1 3, 7X , 13 ) 

508 FORMAT (15X,I1,12X,F6.1,25X,I1) 

505 FORMAT ( 5>, F5. 2 , 5X ,F5.2 ,5X ,F5 .2,5X ,=5 .2 ) 

600 F0FMA7( ' 1 SIX, '»* INITIAL LA M3 COMBAT INFORMATION ** » , 

*////,' LCCA T ION*) 

601 FOEMATd UNIT' ,7X X' , EX , »Y *,4X, • FORCE LEVEL’) 

6 02 FORMA T( 1 X, 13 , 2 X , FT. 1 ,2 X , F T. 1 , 7X ,F5. 1 ) 

605 FGFMAT(/1X,’ ATTRITION IS STOCHASTIC’/) 

6C5 FORMAT! / IX,' ATTRITION IS CET E F M IN 1ST I C ' / ) 

6 G 6 FORMAT( 1>, 'ROLTcS CE T EPMdEO 8Y USSR'/) 

607 FORMAT (IX , 'ATTACK VEHICLE SPEEC IS ',F5.1,' M.P.H.',/) 

6C8 FCFMATdX, 'BREAKFCINT DISTANCE IS ',F6.1,' METERS',/) 

605 FCFMAT( IX , 'DEF EMCER WILL FOT MCVE TO ALTERNATE POSITIONS'/) 

610 FOFMAT (IX , 'DEFENCER WILL MOVE TO ALTERNATE POS I T I ON S • / IX , 

* • ALTER NATE POSITIONS ARE : '/IX , • UNIT » ,5 A , • X' , 8X , • Y ' ) 

611 FORMAT ( 1>, 13 ,3X,F7.1,2X,F T. 1 ) 

612 FOFMAT ( /5 X , ' ATK KILL PROBABILITIES' / 1X,'RANGE',5X, 'P', 

*5 X , 'PHH • , 2X, ' PHM ' , 3X , ' PKH ' I 

613 FCFMAT (2X«I5,5(2X,F5.2)) 

615 FOFMAT ( /5X ,' DEF. KILL FPC EAB I L I TI ES » /IX , • R ANGE • , 4X , ’ P ’ , 
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> » • PHH • , 2X, ' PHM ' , 3 X , ' PKH •) 

615 FORMAT!//, • CURRENT STATUE OF UNIT I VARIABLE OEF IN I T I CMS • . // . 

** C - ALIVE NOT FIRING*./. 4 1 - ALIVE ANO FIRING*,/, 

*♦ 2 - KILLED *,/,* 3 - MOVING’) 

617 FQFNAT (/ / , * VEHICLE SPEED VARIABLE DEFINITIONS’,//, 

*• 1 - S M PH * , / , ' 2-12 MPH • , / , * 3-15 MPH',/,* A - 18 MPH'I 

618 FORMAT ( • 1 * , 1 OX , ' ***** the LANC COMBAT PHASE BEGINS ****«',//) 

627 FORMAT (//IX, ' T I V E = * , 14, IX , ’SECONDS'/ /) 

62? FOFMAT ( / , ' A GG RE SCR UNIT INF C R NAT ION • ) 

62 A FORMAT ( 1>, ’UNIT' ,5X, *X* , 8>, • Y • ,5X , ' F CRCE LEVEL* ,2X, • STATUS’ , 

*2 X, ' LOST- FCT • ,2X , 'TARGETS ' > 

62 5 FORMAT ( 3>,Il,i>,F7.1,2X,Fl. 1,6X,F5. 1 ,?X, II ,6X, F5.3 t 

626 FOFMAT ( 3> , II , 3 X , F 7 . 1 , 2X, F T. 1 , 6 X , F 5. 1,SX, !1,6X,F5.3,3X,3(I1,1X)) 

627 FORMAT (/. * DEFENSIVE UNIT INFORMATION*) 

628 FORMAT ( 1>,*UNIT*,5X, 'X*. E>,*Y *,5X,' FORCE LEVEL’ ,2X,’ STATUS' , 

*2X , • LOST -FCT ' , 2X , 'TARGETS • ) 

62S FORMAT! 2 X, II , 3 X , F 7 . 1 ,2 X , F T . 1 , 6 X , F3. 1 ,9X , 1 1 , 6X , F5 . 3 ) 

620 FOFMAT ( 3>, !1 ,3X,F7.1,2X, = 7. 1, 6X,F 5. 1 , = X, I1,6X,F5.3 ,3X,3< II , 1 X ) ) 

631 FORMAT (1 > , •*&**■£ AGGRESCP FORCE IS ELIMINATED. END OF BATTLE.’) 

632 FORMAT! !>,* = *** CEFENSIVE FC C CE IS ELIMINATED. END CF BATTLE.') 

633 FORMAT ( IX , ****** DISTANCE BETWEEN FORCES IS TOO CLOSE. ’, 

*• END OF EATTLE ****** > 

RETURN 

END 

SUEROUTJNE SETLP 

*■** SUBROUTINE SETUP IS USEC TO READ IN THE TERRAIN DATA AMD 
C R E ATE PARAMETRIC TERRAIN, THIS TERRAIN DAT 4 rfllL BE USED 
\n HEN COMPUTING L I N E-CF -S I CHT EETwEEN TAPGcTS AND CB SERVERS 
AS NELL AS PROVIDING A GRID SYSTEM FCR UNIT LOCATIONS AMD 
MOVEMENT. 

COMMON /HILLS/ X C ( 100 ) ,YC ( 100 ) , PE AK ( 100 ) , AN GH ( 100 ) , SPRDt 100) 
COMMON /HILLS/ EC C ( 1 CO ) , P >X ( 1 C C ) , PYY ( 1 CO ) , PX Y ( 1 00 1 , e AS E 
COMMON /HILLS/ NF ILLS 

COMMON /COVE P/ CX C < 1 50 ) , C VC ! 1 5C ) , CPE AK (1 50 ) , CP XX ! 150 ) , CPY Y ! 1 50 ) 

COMMON / C CVE P / CPXY! 150 .NOVELS 

COMMON / C CUM T R /K H , K KV, , K V , KN , K C- R S , KE L L , K I N T 

COMMON /C-RID/ LST ! 5 ,2) ,NH L( 5 ,4 ) ,L 1ST H( 150 I ,KHR 5P ( 150 ) ,KTREP 
COMMON /C-RID/ LSTC(5»4) ,NC! 5, A) .LIST C ( 400 ) ,KCREP( 150) 

FA 1=5. 1 A 159 
READ (5, 5 CO) NHILLS 
READ (5,510) BASE 
CO 50 1 = 1 .NHILLS 

READ! 5,52 C) XC ( I ) , YC< I ) , PE AK ( I ) , A NGH (I ) ,SPRD ( I ) , ECC ( I) 

50 CONTINUE 

PEA0(5,52C)LS7 

R E A 0 ( 5 , 52C1NHL 

READ (5 ,5CC 1NHTCT 

READ! 5, 5 30! LI STH< I ) ,1*1 .NHTCT) 

65 CO 100 1*1, NHILLS 

AANG L E = ANGH ( I )* FA 1/190 . 

SANG* S IN! AANGLc ) 

CANC- = COS( AANGLE ) 

A=PEA MI) /( FEAMI )-50. ) 

A=ALCG (A) 

B=A* ECC (I )'**2 
SSPD = 5 FPD ( I )**2 

PXX ( I )=-( A*CANG+CANG+B< SANG* SANG) /SSPO 
PYY ( I )=-< A*SANC*S ANG + B *C A N G*C ANG ) /SSP D 
PXY ( I ) = (2.^ SANG *C ANG* < E— A ) ) /SS PD 
KHREF ( I )=-2 147482600 „ r 

C ALL VALUES NCV> IN METERS CN 0 — 10,000 C-RID 
ICC CONTINUE 

REAC(5,5CO) NOVELS 
I F (NC VE L S . EQ . 0 ) C- C TC 200 
DO 1 5 C 1=1, NOVELS 

REAC(5»54C)CXC(I ),CTC( I), CP EAK (I), CP XX (I )»CPYY!I)»CPXY(I) 

KC PE F( I ) =-214748 36 C C 
150 CONTINUE 

PE AO ( 5 » 53 C ) LSTC 
READ ( 5 , 53 C ) NC 
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READ(5,50C>NCTC7 

REAC< 5.53 0 (LISTCCI ), I*I.NCTQT ) 

2 CC * iPEP=-2 1474S3<:C0 
K H = 0 
KHW=0 
K V= C 
KN*0 
KGFS=0 
<E LL=0 
K IM = 0 

5CC FORMAT ( 30,14) 

510 FORMAT ( 1 £ X . r 5 • 1) 

520 FGFMAT(13X,F7. 1 ,3X *F7. 1 * 5 > * F6 . 1* 5X* F 6. 2* 5X. F8. 2 * 4X , F4 . 1 ) 

5 30 FGRMAT< If > f lCi 3) 

5 40 FORMAT (2 f 10.4,2813.7) 

PE URN 
END 

$ UcRGLJT I ME ROUTE 

*** SUBROUTINE ROUTE COMPUTES THE ROUTE OF EACH AGGRESSOR UNIT 
WHEN THE USER HAS SOLE CISC THE CPTION CF INPUTING AGGRESSOR 
ROUTES. IT CALCULATES THE COORDINATES OF EACH INTERVAL ENDPOINT 
ALONG THE ROUTE, HIKING EACH INTERVAL LENGTHOI STANCE MOVED DURING 
A 1C SECOND TIME STEH THE SAME. THE INTERVAL LENGTH IS DETERMINED 
6T THE SPEED THE USER HAS SELECTED AND INPUTED FOR THE CURRENT 
BATTLE . 

COMCN / G P P3 / N6U ,NRU, FL Ul . F0<6 ) ,N0 I (3) »X ICO *200) , Y IC ( 3* 200) , 

* I D IR ( 3 * 2CC), AVSP, I Sr D 

*, IUSTAT ( t ) , I I < t i , LCST( 6, 6 ) , VI SFRA ,VI SFRB , SI ZETK, 

*S!2ETW,NT<6) . NF(cj.£RF,OIS v AX* 

♦ NLCSCl 6 ,E I .VI SFC ( 6 ,6 ) .RMI M K .RNXTK. RMI NTw , RMXTW , CP . TGbFR . LV AFR, 
*PTT(3.3 ) .RF, PCAl t ,c) . ^ 3QA (1 * 6 ) .LOA( 6 ,6) * NA(6) *0FL(6) *POL< 6) 

DIMENSION XLCC(3 .50) tVLCC (3 f 2C)tN(3) 

I F ( ISPD • EC.4 J CST = eD.^c3 
!F(ISPD.EC.3) CST =67 ,U53 
IF(ISPD.c;.2 ) CS~ = 53. 3 43 
I F ( ISPD. EC. 1 ) 057 = 40 352 
DC 110 1=1, NRU 

REAO ( S . 50 C ) M ( I ) 

NL = N ( I )♦! 

DC 1C IN=2.M 

re ad (9, 5 to >lg:s,ylc:cs 
XLCC(I ,IM=XLCCS 
YLCC (I » I M = V LCDS 
CONTINUE 

XLCC ( I .1) =A IC< 1 • 1 ) 

YLOC ( ! ,1) =YIC ( l .1 ) 

ID I R ( ! , !) =C 
NL= N ( I ) 

NUM=2 

CO ICC *J = 1 » N L 

X l = XLQC< I . J - 1 ) -a LOC (:» J ) 

YL = YLOC ( I . J + 1 ) - < LGC (1 • J) 

D IST=SQR7l XL>**2~YL* *2) 

Y = A B S( YU 
Z = Y / XL 

A N C L - A 7 A N ( 7 } 

DEG=ANGL*57 .2953 
IF(J.EC. 1) GC T 0 50 

XLN=(CS7-BXTRA)*CQS( ANGL ) 

C I S T = (DIST+EX7RA )-0S7 
YLN= (CS T- EXTRA)* SIM ( ANGL ) 

X I C ( I *NUM ) = X I C ( I ,NUM-1 )+XLN + XLE 
I F( YL .GT.C. ) GO 1G 2C 
YLN=- YLN 

Y IC ( I ,M!M )=YIC(I iNUM-l)+YLN+YLE 
I F( YL.GT.C. ) GC TO 2C 

IC IP C I » N' U M ) = - IFIX(CEG) 

GO TC 40 

IDI R ( I,NLM)=IFI X (DEG) 
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40 


NUM=NlM+ 1 




50 


X L f=CST* CCS (ANGL ) 
YIN=0S7*'IN<ANC-L> 
IF (Yt.GT.O. ) GO TO 


to 




Y IN =- Y L N 




6G 


IFIDIST. U.CS7I GC 


70 9 C 



70 

80 



<90 



X IC ( I fNUi* ) = XIC< I ,NUM-n + XLN 
YICU »NUM)=YICU iNUR -1 ) + YLN 
I F ( YL.G7. CJ GO 1C 7C 

I 0 I R ( I.NUM ) = - 1FIX (CEG ) 

GG 7C 8C 

ICI R ( I, NUN )=IF!X <0£G) 

C I S T= C I S 7- CST 
NLM = N CM+ 1 
GO 70 6C 
EX 7 P A= C I ST 
XLE=EXTRA*CGS( ANGL ) 
YlE=EXTRA*SIN( ANGL ) 
IF(YL.G7.0.) GC 70 iGO 
Y L£ =- YL E 
CCNT I NUE 



100 

lie CCN7INUE 

500 FOFMAT ( 3 6X 9 1 2 ) 

510 FC'PMAT( 15X»F8.1fl2>fF8.1 ) 

P E7URN 
ENC 

SLEROUT INE L AMCA ( I * J , PC T V IS , C E7RAT, P K ) 

*** SUBROUTINE LAMCA IN CCNJLCTICN WITH THE LCS ROUTINE COMPUTES 

THE DETECTION RATE ( C ETR A T ) OF TARGET J BY th£ OBSERVER I GIVEN 
THE PERCENT CF TARGET VISIBLE ( PCTV I 5 ) TO T He: OBSERVER. 

C C N MON / C P Fl / IPFCIR <6 )t I SECWC <6 ) ,MV7D;R< 6) f X< 6) , Y< 6) , SPO( 6) 

T C FAC 7 = 1 • C 
2 E POl.= C • GO CO 1 
FA 1=3. 14159 

7 C=( ISECWCtI )*PAI / I EC. 0/2.0 

fiEE*(l .0/ <2.0* <S IN (0 )-C*CCS( C ) ) ) ) 

IF (ABS( BE9). LT.2EPCL) EBe=O.C 
A AA=(-B 6 t >*COS (0) 

IF (A6S( AAA). LT.ZERCL ) AAA=0.0 
C 7ANG= A T A N2( ( Y ( J )- Y( I) ),()( J)-X(Ii)) 

IF(CTANG.LT.-PAI/2.ANO.CTANG.G7.-PAI ) C'TANG = 2* PA I + CTANG 

F0= I PRC I P (I)*PA!/1S0.0 

I f ( (PD*CMANG ) .GE. C.O) GCTC 1 

I F ( P 0 . L 7 .C.O ) GGTC 9 

ANGLE=2*FAI+GTANG-FD 

GG7G 10 

9 ANGlE = 2*F AI+FG-GTANG 
10 IF (ANGLE .G7. PAI ) ANGlE=2* FA I -ANGLE 
GOTO 2 

1 ANGLE= A 6 S (FD-C7 ANG ) 

2 IF(ANGLE.GT.D) GC 3 
CLF=PD+C 

CLCR=PO-C 

ANCLFT=C7AN3+ ( 15. C*PAI /I 6 C. ) 

I F ( ANGL FT . GT .CUP) ANGL FT = CUP 
ANGLRT = 01ANG-( 15.* FA 1/180 .) 

IF (ANGLRT.LT.DLCW ) ANG LR T =OLC W 

FK= EBB* A ES (A BS (S IN ( ANGL FT J ) - A ES ( S IN ( ANGL RT ) ) )+ AAA* < ANGL FT- 
* ANGLR7 ) 

IF(PK.LT.C.O) GO TC 3 
IF (PK.G7.1.0 ) GC TC 5 
GO TO 8 

3 P K = 0 • 0 

0 E 7RA T =0 .C 
GO TO 6 
5 PK = 1 .0 

8 RANGE=SGRT (( X( J)-X< I ) )**2 U Y ( J >-Y (I ) )**2 ) 

RR=0.0C1*RANGE / PCTVI S 
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T0ANG=ATAN2( (Y(I)-Y(J)),(»(I)-X(J))) 

A0 = MVTCIH J)*PAi/IS0.0 
HORVEL=AES(SPD<J )* SI N( TC A NG- AC ) I 
HOPVEL = HGRVE L* 16 C <5 . 3 / 2 600 .0 

D E N0M= 1 .A52+TCFAC7MQ.5 < ?7 6 + 2 . 188* (RR**2 )-0 .503 8 *H0RVEL ) 
IF(DENCM.LE.ZEPCL) CENC^=iERCl 
CETRAT = 0.CC3+1 .CSo/OENCM 
CSTRAT=CETRA7*FK 
6 RETURN 
ENC 

SUERGUT I N E ELcV (X , Y , 7MAC ) 

*** SUBROUTINE ELEV CE7ERMINE E THE TERRAIN ELEVATION FOR A GIVEN 
SET GF X, Y COORDINATES. THIS FUNCTION IS USED IN CONJUNCTION 
WITH THE LCS SL6RCCTINE IN CONFUTING L IN E-OF-S IGHT BET «EcN 
OBSERVER AND TARGET. 

COHmjN /HILLS/ XC(IJC) ,YC UOO ) ,P=AK< 100) ,ANGH(100 ),SPRC<100) 

COMMON) /HILLS/ ECC( ICC), P >X( ICC) , PYY ( 100) ,PXY( 100) .BASE 
COMMON /HILLS/ NHILLS 

COMMON /GRID/ LST(r.A) ,NH l( 5 , A ) ,L 1ST H (15 0 ) ,KHREP(150 ) ,K7R£P 
COMMON /GRID/ LS7C (5*4) ,NC( 5, A) ,LIST C(ACO) ,<CREP( 150 
DATA 3SIZE/100C./ 

** FUNCTION TC COMPUTE TERRAIN ELEVATION FCR GIVEN X, Y COORDINATES. 

Z m/X = BAS E 
I X = 1+ I F ! X(X/GSIZE) 

iy=i+if:x(y/gsizei 

I F (NHL (IX,IY).EC.u) GO TO 150 
L S = LST ( I X , I Y ) 

LENC=LS+NHL( IX,IY)-1 
DC 100 L - L S» LEND 
I =L i STH ( L I 
Q X =X-X C ( 1 ) 

CY = Y-YC< I ) 

OXSQ=OX*CX 
OYSC = Q Y * C Y 
CXY=CX*CY 

F ACTCR= P X X ( I l*CXSC + PYY( I I =*QY SC + PX Y( I )*CXY 
IF(FACTC?.LT.-3. ) GO TO ICO 
HT = PEAK ( I )»c xp (FACTOR ) 

I f (FT.LE.ZMAX) GO TC ICO 
ZMA X=HT 
ICC CONTINUE 
150 T v AC=Z MAX 
RETURN 
ENC 

SUEROUT INS STCCh ( I , RANGE , A ) 

*** SUBROUTINE STCCH DETERMINES THE ATTRITION COEFFICIENTS WHEN 

A USER HAS SELECTED A STOCHASTIC ATTRITION OPTION. THE CALCULATION 
IS A FUNCTION CF THE ORIGINAL STCCH A ST IC A LLY DETERMINED ATTRITION 
COEFFICIENT AS WEcL AS A FUNCTICN OF RANGE. 

COMMON / GR P6 / AL PHA ( 6 ) 

COMMON / G RP3 / NEC , N RU , FL ( 6 ) , F C ( 6 ) ,N0 I(3),XIC(3,200),YIC(3,2 00) , 

* 1 CIR( 3, 2CC) , AV SP, ISPC 

*, IUS7AT ( 6) , I I( 6) ,LCS7( 6,6 ) ,VI SFRA ,VI SFRB, SIZETK , 

* SI 2£Tw, NT (6) , N F ( 6 ) ,SRF ,qt $MAX , 

*NLCSC( 6 , 6 ) .VI SFR( 6 ,6 ) ,RMI FTK . RMXTK, RM NTW ,RMXTW ,0P ,TCWFR. LVAFR, 
*PT7(2,3 ) ,RF, FCA (6,6) ,A?0 A (6, 6 ) ,L0A( 6 , 6 ) , MA ( 6 ) , OF L ( 6 ) ,POL( 6) 

IF (I. EC. 2) GC TC 10 
A = ALPHA ( I )*( ( l.C-RANGE/RM>TW)**2) 

GC TC 20 

1C A = A LPH A ( I ) <•( (1.0-RANGE/RMXTK)**2) 

20 RETURN 
END 

SUERC'JT INE ETK ( I , RANGE ,T ) 

*** SUEROUTINE ETK CCMPUTES TFE EXFECTED yi M = FOR A GIVEN FIRER TO 
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nnnn o 0000 



KILL A GIVEN TARGET. THE CALCULATION IS A FUNCTION OF RANGE, 
TIME OF FLIGhT FCP A RCUNC ANC FIT AND KILL PROBABILITIES FOR 
THE FIRING WEAPON SYSTEM. IT IS A N'JMeER THAT IS USEC IN THE 
COMPUTATION CF THE DETERMINISTIC ATTRITION COEFFICIENTS. 

CCMMCN /GR P2 / TA ( 2 1 , T 1 ( 2 ) ,TH ( 2 ) ,TM( 2 ) » Tr 1 ( 2 ) , TF 2 ( 2 ) ,TF3 ( 2 ) , 
*P(2,6),PFF(2,6),FFM(2,6),FKF(2,6),TF(2) 

IF (I. SO. 2) GO TO 5 
TF ( I ) = TF 1 ( I ) 

GOTO 6 

5 IF(RANGE.GT. iCCC.C) GO TC T 

TF< I J*TF1 ( I ) — ( T F 1 ( I )*( 100C.C-RANG E) /I COO. 0) 

GO TO 6 

7 IF (RANGE *GT. 2CC0.G1 GO TO 8 

TF( IJ=TF2 ( I )-( (TF2( D-TFK I ) )*( 20 CC.C- RANGE ) /I 000.0) 

GO TO 6 

8 TF(I) = TF2(I)-( (TF3 ( I )-TF 2(1) )* (3000. C-RANGE ) /I 000 .0 ) 

6 J = (PANGE**250.0 ) / 500 • C 
IF(J.GT.« ) J =6 

T=T A ( I ) +T 1 ( I )-TH( I ) + ( ( TH( I )+TF(I ) ) /PKH ( I , J) ) + ( ( TM ( I ) +TF ( I ) ) / 

* FHM ( I , J ) l=M(l.C-PFMI,J) l/PKM ItJ ) + PhH(I ,J)-P( I, Jl ) 

RETURN 

ENC 

SUERCUTINE SOR T ( I ,M) 

*** SUERCUTINE SORT IS USED TO SORT TARGETS IN ASCENDING RANGE 
OFCER. THIS IS USED TO DETERMINE THE PRIORITY OF A TARGET 
FGP FIRE ALLOCATION. 

CCMMCN /G FP5 / L OT ( 6 , 6 ) , R 0 T ( 6 , 6 ) 

DO 1C J= I , M 

IF(RGT( I , M ) .GE.RCKIfJ )) GCTO 10 
R = P 0T( I , J) 

NN - LOT ( I ,J> 

R QT ( I , J ) = P OT ( I , M ) 

LOT ( I , J) =LOT (I ,M) 

ROT ( I , M ) =R 
LOT (I, M) =NN 
1C CONTINUE 
RETURN 
END 

SUeRCUT INE KOVER( ZC.TMACT , S I 2 E T , Z T , S ,H TS , Z S , VI SF R T ) 

*** SUE P CUT I NE KCVER CET ERMINES W F AT PORTICN OF A PARTICULAL TARGET 
IS COVERED BY THE TERRAIN BETWEEN THE TARGET ANO OBSERVER. 

THIS NUMEER IS USEC IN THE DETECTION AND ATTRITION COMPUTATION. 

I F (S.EQ. 0. ) GC TC 10 

IF(HTS.GE.ZS) C-G TO 20 
HEXT=ZO+ (HTS-ZO) /S 
EV IS T= A MAXI ( FEXT ,TMACT) 

I F ( EVI ST.GE.ZT ) GO TO 20 

IF( EVIST.LE.ZT-SIZET ) RETURN 
VI S= ( ZT- cV I ST )/SI ZET 

IF(VIS.LT • VISFPT) VISrRT=VIS 
RETURN 

10 IF(HTS.LT.ZO) FETUFN 
20 VI SFRT = C .0 
RETURN 
ENC 

SUBROUTINE LOS ( XA , YA » TMAC A,TMICA,SI Z EA , XB , YB , TMAC B ,T M I CB , S I ZE B , 

♦ LA7CB, LB7CA, V ISFRA,V ISFPB ) 

*** THIS SLBFCIT !NE V^AS WRITTEN BY FRCFESSCR JAMES HARTMAN, NAW 

POSTGRACtATS SCFCCl. IT CGMPLTES A PERCENT OF A TARGET VISIBLE 
TO A PARTICULAR CESERVER, GIVEN THE CCCRD INATES OF BOTH 

COMMON /HILLS/ XC ( 10 C ) , YC ( 10 C ) , P E AK ( 1C C) , ANGH ( 1 0 0 ) , SPRO ( 1 00 ) 
CCMMCN /HILLS/ ECC ( 100 ), P >X ( ICC ), PYY ( 100 ), PXY ( 100 ), BASE 
COMMCN /HILLS/ NHILLS 

COMMON / CCVE R / CX C ( 1 50 ) , C YC ( 1 5C ) , CPE AK ( 1 5C ) , CP XX ( 15C ) ,CPYY( 15 0) 
COMMCN /COVER/ C F X Y ( 15 0 ) , NC V E L S 
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c*** 



c*** 



IOC 



lie 

12C 
125 
13 C 



1AC 
1 5 C 



16 C 



170 



180 
19 C 



20C 

C GRID 
C 



CGKMCN /COINTR/KHt KHW , KV, KN , KGRS , KELL, K I NT 

COMCN /GRID/ LST ( 5, A) »NHl( 5*4) ,LISTH{ 1501 , KHR E P ( 1 5C ) ,K7R6P 

CO MON /GRID/ LSTC ( 5 ,A) , NC( 5, A ) ,LIST C (AGO ),KCP EP( 150 ) 

CIMENS ION IG X ( ICC) » I GY(IOC) » I ELI 1 GO ) »CS1 (130) ,CS2 <100 ) 

DATA GS I / E/l ooc . / 

SUBROUTINE TO CCNFLTE FRAC1ICN VISIBLE FOR OBSERVER TARGET PAIRS 

V I SFR A = 1« 

V I SFP3 = 1 . 

X EA=XB- X A 

V £ A=Y'B- Y A 

IF !<X3A. EO.C- . ) . AM . ( YBA. EC.0 . )) RETURN 
IF ( SIZEA+TMICA.LE.O.) GO TO 510 
IF( S I2EB+TJ/ ice. LE.0.) GO TC 510 

if(tmica. i r.o. ) visfra=i.o+tmica/$izea 

IF(7M;C3.L7.C. ) VISFR3 s l.Q+TNIC3/SIZE3 
Z A - T , V A C. A -t 1MICA + 5 1 Z £ A 
Z 6 = I HA C B + iNICB + SIZEE 
K.TRE?=KTPEft 1 
Z EA-ZB-Z A 
xei£C--?eA*XEA 
Y 8 A Y 8 A* YE A 

X YEA-XEA*YEa 
T l»0 XBA = 2 X E A 
T KG YBA *2 . * Y f A 

CONFUTE GRID SCUAFtS CROSSED BY A TO e LINE 
NC-RSQ=0 

IFIxeAi 310, 95,100 
X3 A = G. 1 
ISGX=-L 

XI NC=G3 I ZE/XfiA 
GO TO 12 C 
I S G X= 1 

XINC=-GSI ZE/XEA 
IF(YBA) i*),125,:3C 
YBA-C.l 
I SG Y = -l 

YINC =GS JZE/Yf A 
GO ’0 15C 
I S C- Y = 1 

YlNC=-0 SI Zc/veA 

I X = 1 + IF IX (XB/GSIZK) 

I Y = l-*- i F I X >' Yb /G SI ZE ) 

XNEXT=C-SIZS« <FLGA T < IX) +0. E*< ISGX-l. ) ) 

YNcX'=GSIZE* ( FLCA* ( IY ) +0 . £*< ISGY-1. ) ) 

X SYE»'= ( XE- VEX 1) / XS A 
YSTEP= ( Y6-YMEXT ) /vEA 
NC- F SO=NG r SC* i 
IGXINGR EC ) -• x 
IGYINGRSC )= I Y 

IF((XSTSF,GT.l.) .ANO. ( YST EP.C-T.l. )) GO TO 200 
IF (XSTEP- ASTER ) 17C, 18C, 1 EC 

IX= IX * ISGX 
XS7E? = V . STEP+XlNC 
GO TC lEC 
I X = IXfl SGX 
XSirP=.<SlEP + .«INC 
IY= JY4I 3GY 
Y 3 7 F F =Y STEP* Y INC 
GO 1C IcC 
KGRS=KGRS+NOPSC 

SCUARE LIST NCW COMPLETE IN IGX, IG Y WITH NGRSO ENTRIES 



FINC WHICH COVER ELLIPSES TOUCH THE A TC B LINE, 
CHECK SLEYA7ICNS A1 SI AND S2 FCR EACH SUCH ELLIPSE 
NEL S=0 
C H1MAX = C • 

IF(NCVELE.ta.C) GC 10 270 
DO 2 1 C K=1,NGFSC 
I X * IGX ( K ) 

I Y= IGY (K ) 

N= NC ( I X , IY ) 

IF(N.EC.O) GC TC 26C 
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LS=LSTC( IX,IYI 
LENO=l$*r>-l 
CQ 25 C L - L S » L ENC 
KELl=KELL+ 1 
IC =L I STC (L) 

IF (KCPEP ( IC) . EC.KTPEP ) SC TO 250 
KCREP(IC) =*TR£C 
RX=XA-CXC < !C> 

RY=>A-CYC ( IC) 

PPXX=CPXX( IC) 

PPY Y= CPYY ( 1C ) 

PPX> = CPXY( IC J 

AA=PPXX*X6 ASG+PFYY* YS ASC+PPXY*XY8A 
BE= PPXX*TWOX EA*RX+?PYY* TWCY SA*RY+P°XY*(RX*Y8A+RY*XBA ) 
CC=PP>X*RX*RX+FPYY‘RY*RY+ PPXY^RX* RY-1 *0 
APG=EE*BB-‘.0*AA*CC 
IF(APG.LE.C. ) GC T C, 2 50 
SC = SCR T (ARG» 

S1=-(BEM$Q J/IZcC^AA) 

S2 = ( SQ- £8) / ( 2.0*AA ) 

I F ( Sl.GE. 1. ) GO TO 25C 
- IF (S2.LE.0* ) GO TO 25C 
IF(S1. IS.O. ) GO TC 510 
IF (S2.GE.1 • ) GO TO 51C 
C*** CHECK LOS AT SI AAC S2 
K I N T = KI NT + 1 
CFK=CPEAK (IC ) 

XS= > A+S2* >6A 
YS=TA+S2*‘t8A 
CALL ELEVtXS tYStHTS > 

HI S=H7S+C PK 
Z$ = 2A+S2-*Z8A 

I F ( l ATQ6 • 8 C . C ) GC TC 210 

CALL KCVER ( Z A , T M A C 8 » S I Zed ,?S »$2*HTS»ZS»VISFRB) 
IF(VISFRE.LE.O.) CQ TC 510 
I F ( l B TO A • EC * C ) GO TC 220 
S=l.-$2 

CALL XOVEM ZBtTMA (At S IZEa tZA ,S,HT$t ZSt VI SFRA) 
IFIVISFRA.LE.O.J GO TC 510 
XS*XA+S1*X8A 
Y S= TA+Si* YB A 
CALL ELE\ j (XStYS*HT$i 
HT S = HTS ♦ C PK 
Z S = ZA + S1*Z8A 

I F ( LA TO 8# EG. C) GO TC 230 

CALL KCVEP (ZAtTMACBt S IZ53.ZB »S1«.HTS»ZS»VISFR8) 
IF(VISFR8*LE«Q.) GC tq 51Q 
IF( LETOA, EG.C ) GC TO 240 
S-l* O-Sl 

CALL KOVER (Z3 t TMACA, SIZ5A fZA t S .HTSt ZS, VI SFRA)- 
IF ( V ISFRA.LE.O.) GO TG 51G 
NE L S = NE LS + 1 
IEL(NEIS) =IC 
CS1 (NELS )»Si 
CS2(NEL$) - S2 

IF (CPK.0T,»CH7MAX) CHT MAX = CPK 
CCNTINUE 
CONTINUE 
C*** ALL ELLIPSES CHECKED 
C*** START CN THE r* I L LS 
270 CO 6C0 K*lfNGRSO 
IX= I C- > ( K ) 

I Y= I G > ( K ) 

I F < NHL ( IXyIY)*EC*C) GC TO ECO 
LS=LST (IX, IY) 

LEND=lS*NHL(IX,IY) -1 
DO 5CC L = L S , L E N C 
I = L I ST H ( L ) 

I F ( KHR E P ( I ) .EC.KTREF) GO TO 50 0 
K HP E P( I ) =K TRE P 

C*** PFOCESSINC FOR FILL I STARTS HERE 
KH=KH+ 1 



?1C 



22 C 



23C 



24C 



25 C 

26 C 
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c 



28C 

2CC 

C*** 

C*** 



2 5 C 



28C 

28C 



C*** 

ACC 



COMPUTE W =TOP CF THIS HI LI ALONG G-T LINE 

TR>=XA-XC( I ) 

T R Y= YA - YC ( I ) 

T P XX=PXX ( I ) 

T P V Y=P Y Y ( I ) 

T F>Y = PXY (I ) 

FQ=TWOxeA*TFXX*TRXMWCYe.3*TPYY *TRY+TPXY*(TRX*YSA h 
GQ=TPXX* X6ASQ+7PYY* VBASC+TPXY* XyBA 
IF (GC .EC.O. ) GG TO 500 
W=-FC/ (2.*GC) 

I F ( A8S ( W ) • G7 • 5 • ) GO TO 500 
FSG-FQ*FC 

EG = 7PX X* 7P X* TR X+7PY Y*TPY*TRY+T FXY*TRX*TRY 

PC VER=E Q-FS 0/ (4.*GC ) 

IF ( PCW ER .IT. -3* ) GO TO 500 
HHV» = PE AK (I )*6XF(P0W fR ) 

KH W=Km w + I 

IF (HHW.L6.eASE) GO 10 5CC 
ZW«ZA+W* zeA 

IF<( W.LT.O. ).QR. (W.GT.l .)) GO TC300 
IFtHFW .GE. ZW ) GO* 10 5 1C 
C V H 7 V=0 • 

IF(NELS.EC.O) GO TC 30C 
DC 280 MsltNELS 

IF( (CSl(M).GE.W ).GR.(CS2(M) .LE.W)) GO TO 280 
IC=IEL(MJ 

I F( CVHTK.LT. CPE A K( IC ) ) CVH T W=C P E AK ( I C ) 

C CM IN LE 

IFKFHW + CVHTW) .GE.ZJr) GC TC 510 

I F(hFW + CHTMAX.LT .AM INK ZA-SIZE A, Z9-SIZEB) ) GO TO 
IF WE GET TC HERE THEN NEEC TC FIND LOWEST SIGHT LINE C\ 
NEWTON ITERATION A TC 6 GIVING VISFR8 
IF(LATOB.cC.O) GC T C 400 
K V - K V+ 1 
V = W 

HHV^HH W 
NC T = C 
FV = FC*V 
7 WGG V= 2 • *G C* V 

FCNV = Z A+HHV* ( TWGGV* \ + FV- 1. ) 

K N = KN+ 1 

FAC7CR=(7WCGV*TviCGVh2.*(GC+7W0GV*FG)+FSQ) 

DFCNV=hHV*V*FACTCR 

IF ( AES ( DFCNV ) .LT.1.E-1C) GC TC 250 

V= V-FC N V /DFCNV 

F V = FC* V 

TWGGV=2.*GC*V 

P0V6R = EQ + FV+GQ*V*V 

IF ( PCW ER .LT. -3.) GO TO 400 

HHV=PEAK(I)*FXR(PCW ER ) 

OHHV=HHV*(FO+TWOGV) 

ELV=ZA+CHFV*V 

IF (ABS(HHV-ELV) .L7.1.) GO TO 350 
N07 = NC T + 1 

I F (NCT.LT.10 ) GO TO 33C 

IF ( ( V.LT.C. ) .CR. ( V.C-T.l. ) ) GO TC400 

C V HTV= 0 • 

IF (NELS. EC.O ) GO TO 390 
DC 260 M=I >NELS 

IF( (CSl(M).GE.V ) .OR. (CS2(M) .LE.V) )G0 TO 360 
IC=IEL(N > 

IF( CVH7V.LT.CPEAK( I C) ) C VHT V=CP£ AK ( I C ) 

CONTINUE 
HTV=HHV+CVF7 V 
Z V = ZA+ V* ZB A 

CALL KCVER(Z/»TMACB ,SIZEE,ZB, V ,HTV, ZVtVI SFRB) 

IF (VI SFP8.L5.C. ) GC TC 510 
NEWTCN I7ERATICN B TC A GIVING VISFRA 
IF ( AeS ( V J.CT.5 . )G0 7C 4CC 
I F ( L3T0A .EC.O) GC TC 5 C C 



TR / *XB A ) 



SCO 

ER HILL 
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43C 



45C 



«ec 

4SC 

ICC 

eoc 

5 1 C 



KV=KV+ 1 

V = * 

V M 1 = V- 1 • 

HHV=hHW 
NC T =0 

F V = FG* V 
TW0GV=2.*GC*V 

FC N V=Z B + FHV* ( (FO+TWOGV . ) 

K N = K N'+ 1 

FACTGR= ( TWOGV^TWCG V *»2* + ( GQ+TWO GV*FQ )+FSC) 

DFCNV=hHV-*VMl*FACTC F 

I F ( A8S ( C FC N V ) •LT.1.E-10I GG TC 450 

V=V-FCNV/CFCNV 

IF<A6S<V).GT.5.)G0 TO 5CC 

VM 3 = V- 1. 

F V = F C* V 

TW CG V= 2 • *GC* V 

PCWER = EG ♦ F V + GQ*V* V 

IF < PCW EF • LT . -3.) CO TC 500 

HH V= PE AK ( I ) * EX P ( FCfc IR) 

OHHV=HHV * { FC+TkCGV ) 

EL V=Z3+ChhV*VNl 

IF ( AB S ( HH V- EL V) .LT.l.) GO TQ 450 
N C 7 = NCT •+ 1 

IF (NCT.LT.1C) GC TC 43C 

IF (( V. LT.O. ).CP. (V. GT. 1. ) ) GO TC 500 

C V FTV= 0 • 

IF (NELS.6G.C) GO TC 490 
CC 430 M=1»NELS 

I c ( (CS1 (P ) • G E «V ) «QR • ( CS2 <M ) .L E • V ))G0 TO 480 
IC= I E L < M ) 

IF( CVHTV.LT •CPEAK(IC) ) C VHT V=CP E AK ( I C ) 
CONTINUE 
HT\=HH V+CVHTV 
ZV=ZA+V*Z3A 
S = - VM 

CALL KOV ER ( ZB t 7MAC A ,SIZEA»ZA t $ ,H7V*2V,VI 5FRA) 

I F (VI SFRA.LE.O . ) GC TC 510 
CONTINUE 
CONTINUE 
RETURN 
VIEFRA=0. 

V 1$FR6 = C . 

RETURN 

END 
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APPENDIX C 



COMPLETE INPUT DATA SET 
for the 

SMALL-UNI T AMPHIBIOUS OPERATION COMBAT MODEL 

The small -unit amphibious operation combat model consists of two 
phases of combat, ship-to-shore and land combat, and requires data 
input for each of these phases. The data set that follows is divided 
into two parts: the first part consists of all data used as input for 

the ship-to-shore phase of combat, and the second part consists of all 
data used as input for the land combat phase of combat. The input data 
set was designed to be self-documenting in that the input variable 
names or descriptive phrases are listed alongside the data being used 
as input to the model. The purpose of this documentation was to assist 
the user in associating the input data with their respective input 
variables. A complete input data set follows. 
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IPRIM 


= I 


C FOR EACH 


TIME STEP £ 


1 FOR ENO 


OF eATTLE 


. 


L V A 1 5 


5 POM A X = 


40.00 5 


POM I N =10. 


5 HTMAX = 1. 7 


HTMIN = 


0.6 WIDTH 


= 3.533 




TTS = 


10. 


LENGTH GF 


EACH TIME STEP IN SECONDS. 




TANK 


MAX FANGE = 


I50C. 


A 7GM MAX RANGE = 2000. 


ATGM MIN 


RANGE = 


200. 


7SIGV 


T ART M = 


15. 


SARTM = 30 


. TVEL = 350 


. SVEL = 


35C. 




- -5 C 

*■ a • 




500 . 


100 0 . 


2000 . 


5000. 


10000. 


0. 




4_ . 




5. 


20. 


25. 


25. 






TS IGH 






500. 


1000 . 


2000. 


5000. 


1C000. 


0. 




c . 




c. 

o' . 


20 . 


25. 


25. 






TM6NF 


^ 25. 




5 CC • 


ICOC. 


20 GO. 


5 GOO. 


10C00. 


0. 




I. 




5. 


10 . 


15. 


15. 






SS IGV 


= 25. 




25C. 


500. 


1000. 


2500. 


5000. 


10000. 




. 




5 . 


7.5 


14. 


15. 5 


17. 


20. 


SS IGH 


- 2 5. 




250 . 


500 . 


1003. 


2500. 


5000. 


10000. 




C. 




2 . 


7.5 


14. 


15. 5 


17. 


20. 


CE F . 


WEIGHTS 


A5SIGNEC TO 


WAVE CNE 


= 2. AND 


WAVE TWO 


= 1. 






.NO. OF 


L V A 


PER fcA VE * 25 . 


20. 


15. 


10. 


5. 


SHOP 


E DCF. 


TANK 


A SSETS 


= 


10. ATGM 


ASSETS = 


10. 





ATTAIT COEF. fcf 
A 1 I R 1 1 C w h h • F C R 
AIMED FIFE 
GAIN 
GAMMA 
NO. CF HI 
9 A 

HI l C CATAiXC = 



AGGFESSOF FORCE FIRE: ALPHA> 0.00006 ALPHA2 = 0.CCC08 

utFENSiVt FORCE FIRE: BtTAi * 0.0007 EETA2 = 0.0009 

ATTRIT CCEF • : WEETAtl) = C.0005 WBETA12) = 0.0007 



LSI (5*4) = 



11 = 0.32 

= 50.0 DELTA = 

LLS IN TERRAIN MO (EL = 
SE = C.C 

2C00. YC - 1100. P EAK = 


IOC. 0 
46 

170. A NGH= 


0. 1 


SPRO = 


999.9 


I SCO. 


2200 . 


150. 


30. 




350. 


: cco. 


1 9 GO . 


150. 


130. 




2C0. 


:acq. 


1400. 


150. 


0. 1 




200. 


2450. 


17 CO • 


130 . 


80. 




5C0. 


;7co. 


ieco . 


136. 


90. 




500. 


3 200 . 


1650. 


140. 


150. 




6CC. 


4 3 00 . 


1300. 


120. 


160. 




400. 


3 750. 


17 5C. 


150. 


0. 1 




660. 


* 150. 


16C0 . 


150. 


160. 




550. 


3 2 CO . 


2150. 


130 . 


25. 




500. 


' 6C0. 


1700. 


170. 


45. 




3 CO. 


* SCO . 


1500. 


170. 


0.1 




3CC • 


i £ C 0 . 


2600. 


170. 


90. 




350. 


3 400 . 


2850. 


150. 


120. 




2C0. 


3100. 


2700. 


150. 


150. 




250. 


2 500. 


24CC. 


150. 


0. 1 




250. 


2 6 5 0 . 


2850. 


150. 


160.0 




4C0. 


27CO. 


26C0. 


150 . 


130. 




370. 


3 ECO . 


2 2 C 0 • 


150. 


0. 1 




230. 


4 5C0 • 


26C0. 


150 . 


90. 




280. 


3 6 CO. 


2 8 CO. 


150. 


145. 




500. 


2 700 . 


23C0 . 


190. 


25. 




350. 


3000 . 


33 CC. 


170 . 


15. 




400. 


3150. 


27 5 C. 


130. 


0.1 




350. 


3750. 


32 CO • 


150. 


10. 




850. 


3 ECO • 


2 8 C 0 • 


150. 


0.1 




650. 


3 ECO. 


36C0. 


150. 


160. 




320. 


4 150 . 


2950 . 


170. 


30. 




220 • 


1650. 


2 ICO. 


150. 


30. 




300. 


2250 . 


2 ICO • 


180. 


150. 




220. 


2 2C0 • 


1700. 


180 . 


80. 




150. 


2 ECO • 


1 4 C 0. 


150. 


165. 




500. 


3 600. 


2650 . 


150. 


20. 




400. 


*150. 


3 CCO . 


130. 


160. 




200. 


4*50 . 


21 50. 


120. 


140. 




350. 


4850 . 


3600. 


150 . 


45. 




40C. 


4500. 


40C0. 


150. 


60. 




450. 


4850. 


2450. 


150. 


150. 




440. 


4900 . 


3000. 


130. 


45. 




250. 


4 7C0 • 


2 4 CO. 


150. 


50. 




400. 


4900 . 


20C0 • 


130. 


0. 1 




25C. 


2 5C0 . 


11C0. 


170. 


0. 1 




950. 


3350. 


4 2 CO. 


150. 


45. 




450. 


4C00. 


2200. 


150. 


45 • 




280. 


3 9 CO. 


22 CO. 


150. 


0. 1 




300. 


= 0 


C 0 C 


0 0 


1 


7 


18 



ECC = 



27 



8.0 
2. G 

2 • 

2. 5 
2.2 
2.2 

3. 



1 . 

2 . 

2 . 

1 . 

1 . 

2 . 

1. 

3. 

1. 

1 . 

1. 

2 . 

2 . 

2 . 

2 . 



2 . 
2 . 
2 . 
1. 
1. 
2 . 
2 . 

1. 

2 . 
1. 
2 . 
1. 
1. 
2 . 
1 . 
8 . 

2.5 

2 . 

3.5 
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uiunnoo owunouiwuicooo oooouiunn Ovn 



0 


33 


39 


NHL(5,4f * C 


C 


0 


0 


A 


14 


NO. OF HILLS 7G7AL 


= 




L I S1HC I 1 = 1 


2 


3 


6 


J 2 


33 


8 


Q 


10 


8 


42 


-5 


16 


17 


18 


2 


31 


11 


46 


2 C 


L 1 


42 


45 


46 


26 


14 


25 


35 


4 4 


26 


40 







NOVELS = 



62 


0 


74 


77 


62 


93 


0 


0 


6 


11 


9 


6 


12 


C 


3 


6 


10 


9 


101 

4 


43 


1 


3 


4 


5 


11 


31 


43 


1 


6 


7 


33 


43 


10 


12 


13 


9 


30 


23 


15 


3 


14 


15 


20 


*3 


6 


23 


11 


7 


20 


22 


34 


35 


44 


45 


12 


34 


35 


36 


40 


41 


23 


15 


23 


24 


25 


15 


27 


28 


29 


24 


22 


23 


29 


35 


36 


37 


38 


39 



0 
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S3 

0 

S 

30 

7 

11 

14 

IS 

16 

1 2 

14 

26 

27 



ATTPIT 


VAR = 


0 CSEED 


as 


1 42257.0 


(FCR 


AGGRESSOR 


FORCES 1 


BETA CIST. 


INPUT 


PARAMETERS : 


FP 


= 


7. C 


QQ " 


21.0 


ATTPIT 


VAR * 


0 CSEED 


= 


122457.0 


(FOR 


DEFENSIVE 


FORCE S ) 


BETA CIST. 


INPUT 


PARAMETERS : 


FD 


= 


21.0 


QO = 


7.0 



NC. CEP UNITS = 
PM I NT K = OCCC.O 
TYPE CF ROUTE = 



03 NO. 4 TK UMTS = 03 
FMXTK = 250 C. 0 FIMNTW 
1 VEHICLE iPEEC = 4 



= 0500. C 



RMXTW = 4000.0 



XIC(1,1) = 2C00.0 


Y I C ( 1 » 1 ) = 


1 SCO.O 














1SC0.C 




2400. 0 














15C0.C 




2 ICO. 0 














NO. CF NODES FOR 


ROUTE 1 = 


0 1 














XL0C(1,1»= scco.c 


YLCC ( 1 « 1 ) - 


2 5 CO. 0 


FCR 


NODE 


1 


CF 


F CUT E 


1 


NO. CF NOOES FOR 


ROUTE 2 = 


01 














XL CC ( 2 1 1 ) = 4 SCO . 0 


YLGC (2.1) = 


2150.0 


FOR 


NODE 


1 


OF 


ROUTE 


2 


NC. CF NODES FOR 


RCUTE 2 = 


02 














XLCC ( 2 » 1 ) = 22CO.O 


Y LQC ( 3 » 1 ) = 


1 ICO. 0 


FOR 


MHOc 


* 


CF 


F CUTE 


3 


>LOC ( 2 *2 ) = ^6C0.C 


YLOC (3.2) = 


1750.0 


FOR 


MODE 


2 


OF 


F OUT E 


5 



DEF UNIT X = 

T 

LOCATION 

ALT. POS. 



36CO.C Y = 2 7CC.Q FCPCE10.0 01 R. ISC WIDTH 120 

2800.0 2 2 CO. 0 LEVEL 5.0 C'F ISC OF 120 

26 CO • 0 17C0.0 10*0 K re 180 SRCH 120 

VAR. 0 BREAK PT C5C0. C NC. TIME STEPS FCR MOVE 4 

FCR ALT. POS. 1 
* = - 2 



X A ( 1 ) 


= 4500 


.0 


YM1) = 


2600.0 


X A ( 2 ) 


= 4 5 CO 


• C 


Y A ( 2) = 


2 7 CO. 0 


X 


= 4 6C0 


.0 


Y = 


1 6 CO . 0 


0 . 6 C 


0. 


70 


C. 65 


0.85 


0.65 


C. 


90 


C. 65 


C. 90 


0.60 


0. 


85 


0.85 


G.80 


0.75 


c. 


8 C 


C . 7 5 


C.70 


0.60 


C. 


7 C 


C.65 


C.65 


0.40 


0. 


45 


0.40 


C • 5 0 


C.F5 


PHH C. 


65 


FHN C.75 


PKH C.70 


0.60 


C. 


80 


0. 75 


C.70 


C • 75 


0. 


75 


C. 70 


C .60 


0.60 


c. 


65 


C. 60 


C. 5 5 


C • 4 5 


0. 


50 


0.50 


C.35 


C.2C 


C. 


2 C 


C. 20 


C. 20 
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APPENDIX D 



BLANK INPUT DATA SET 
for the 

SMALL-UNIT AMPHIBIOUS OPERATION COMBAT MODEL 

The blank data set provided with the small-unit amphibious opera- 
tion combat model was designed to assist the more familiar user of the 
model in the development of a new input data set to be analyzed by the 
model. It is patterned after the complete input data set listed in 
Appendix C providing input variable names or descriptive phrases to 
identify the locations of required input parameters. Underlining of 
the spaces following these descriptors is intended to serve as a guide 
for inputing values for the input variables in order that they will be 
compatible with the formatted read statements of the program. The 
blank input data set follows. 
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IPRIN7 = 

LVA'S SFCMAX = 
TTS = 

TANK MAX RANGE = ' 

T ART N = _ 



0 FOR E ACH TIME STEP £ 1 FOR END CF BATTLE, 
S POM IN =„. FTMAX = HTMIN = . WIDTH = 

LENGTH CF SACH TIME STcP IN SECONDS. 

AT GM M A > RANGE = . ATGM MIN FANGE = 



SAFTM 



TV EL 



... SVEL = . 



NC. OF LVA 
SHGRE CEF. TANK 
ATTRIT COE*. F Q P 
ATTRIT CCEF. FOR 
AIMED F IFE 
GAIN 
GAMMA 
NO. CF HI 
BA 

HILL D AT A : XC 



TS IGH 


= 


-• 


-• 


— • 


— • 




THENH 


- 


- • 


- * 


. . 


— • 


— • 


SSIGV 


- 


-• 


-• 


— • 


— • 


— 1 


S SIGH 


— 


-• 


- • 


- •- 


— * 


— * 


DEF. 


WEIGHTS 


assignSc to 


WAVE ZhZ = 




AND WA7E* 



PER WAVE = 
ASSETS = 
AGGRESSOR FORCE 
DEFENSIVE p CRC E 
ATTRIT CGEF.: 



. ATGM ASSETS = . 

FIRE: 3CP KA1= _. AL PHX2 = . 

FIRE: Bti'Ai = ~ SETA2 = I.' 

WeETA(l) = _. WBt 7A { 2 ) = 



LS7<5,4) = 



L = 


DELTA = 










ILLS IN TEPP/! IN MO CEL = 










II 

OO 
• >- 
<-> 

• 

M O 
O 
LUO 

WJNJ 


11C0. P EAK = 


170. 


A NGH= 0.1 


SPRO = 


5 9 9.9 


18C0. 


2200. 


150. 


30. 




350. 


2 C CO • 


1 SCO. 


150. 


130. 




300. 


2 AGO • 


14C0. 


150. 


0.1 




300. 


2450 • 


1700. 


130. 


80. 




5CC. 


2 7 CO. 


18C0 . 


138. 


90. 




500. 


2200 . 


16 50. 


140. 


150. 




6 CO. 


1 200 • 


1 2 CO. 


130. 


160. 




400. 


2750. 


1750. 


150. 


0. 1 




660. 


4150. 


1600. 


150. 


160. 




550. 


2 2 CO • 


2150. 


130. 


25. 




500. 


4 £ CO • 


17C0. 


170. 


45. 




300. 


4800. 


1500. 


170. 


0.1 




3CC. 


2 2 CO • 


2 6 CO. 


170. 


90. 




350. 


24CO. 


2850. 


150. 


120. 




300. 


2100. 


27 CC • 


150. 


150. 




350. 


25CO. 


24C0. 


150. 


0. 1 




250. 


2650. 


2850. 


150. 


160.0 




400. 


2 7 CO • 


26CC. 


150. 


130. 




370. 


2 SCO • 


2 2C0 . 


150. 


0. 1 




23 C. 


4500 . 


26CO • 


150. 


90. 




280. 


26CO. 


28C0. 


150. 


145. 




500. 


2700. 


2 3 CO • 


190. 


25. 




35C. 


2CC0. 


2300. 


170 • 


15. 




400. 


2150. 


2750. 


120. 


0.1 




350. 


2750. 


2200. 


150. 


10. 




3 50. 


2 6 CO. 


3 8 CO. 


150. 


0.1 




650. 


2 600. 


26C0. 


150. 


160. 




22C. 


4150. 


3950. 


170. 


30. 




220. 


1650. 


21CC. 


150. 


30. 




300. 


2250 . 


2 1 CO • 


130. 


150. 




2 20. 


2 d. CO • 


17CC. 


180. 


80. 




150. 


2 E 00 • 


14C0. 


150. 


165. 




500. 


2800. 


2650. 


150 . 


20. 




4 C 0 • 


M50. 


3CC0. 


130. 


160. 




300. 


4450. 


2150. 


130. 


140. 




350. 


4850 . 


3600. 


150. 


45. 




400. 


45CO. 


4CC0. 


150. 


60. 




450. 


4850. 


2450. 


150 . 


150. 




44Q. 


4 5 CO • 


20C0. 


130 . 


45. 




250. 


4 7C0 • 


24 C 0. 


150. 


50. 




400. 


4 SCO • 


2000. 


130. 


0.1 




250. 


2 5 CO. 


1100. 


170. 


0. 1 




950. 


2250. 


4200. 


150. 


45. 




450. 


4000 . 


2200. 


150. 


45. 




280. 


2SC0 • 


22C0. 


150. 


0.1 




300. 


* 0 


0 0 0 


0 


0 1 


7 


18 



ECC = 



27 
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0 


33 


39 


Nh l ( 5 ? 4 ) = C 


C 


0 


C 


6 


14 


N0« CF HILLS TCT4L 


s 




LiMHin * i 


2 


a 


6 


32 


23 


3 


9 


10 


e 


42 


2 


16 


17 


18 


2 


31 


11 


46 


2C 


21 


42 


45 


46 


26 


14 


25 


35 


44 


26 


40 







NOVELS = 



62 


0 


74 


77 


83 


93 


0 


0 


6 


11 


9 


6 


12 


0 


3 


6 


10 


9 


101 

4 


43 


1 


3 


4 


5 


11 


31 


43 


1 


6 


7 


33 


43 


10 


12 


13 


9 


30 


23 


15 


3 


14 


15 


20 


■a 


6 


23 


11 


7 


20 


22 


34 


35 


44 


45 


12 


34 


35 


36 


40 


41 


23 


15 


23 


24 


25 


15 


27 


28 


29 


24 


22 


23 


29 


35 


36 


37 


38 


39 



0 
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53 

0 

9 

3C 

7 

11 

14 

19 

16 

22 

14 

26 

27 



£ TT R IT VAR = _ 0 SEED = 

8 ET A CIST. INFUT PAPANETEPS: FP = 

ATTRIT VAR =■ CSEED = 

BETA CIST. INFC7 PARA VETERS : FO = 

NC. CEF UNITS = _ NC. ATK JMTS = _ 

PMINTK = PNXTK = . FMfiTW 

TYPE CF ROUTE = __ VEHICLE SPEED =„ 
XI C ( 1 , 1 ) = ._ YlC(lfl) * 



(FCR AGGRESSCP FORCES) 
_• _ QQ = • 

IfZR DEFENSIVE FOR<T=ST 
00 = 

a ._ RMXTW = 



N C. CF NODE? FOa POUT E „ = 

xlcc( i*i)= nccTia 

NC. CF NO 0 E? T3R ROUTE 
XLOC ( 2 ,1 ) = YLCCl2,l 

NG. CF NODE? T5 R R 0 L T E 



XL CC ( 2 ,1) = 
>LCC (2,2)= 
CEF UNIT X = 

I 

LOCATION 

ALT. POS. 
xa m = 

X A ( 2 ) = 

X = 



_ , • _ 

varT*I 



AZV5.Z 



>LCC 73, 1 
Y L CC (3 t2 
Y 



eREAK 

Y A ( 1 ) 

Y A ( 2) 

Y 



I =■ 

) = 

) = 
) = 

PT 



FCR NODE 1 CF ROUTE 1 

FCR NODE 1 OF ROUTE 2 

FCR NODE 1 OF ROUTE 3 
FCR NODE 2 CF ROUTE 3 

FCRCE ._ DIR. WIDTH 

LEVEL ._ C c OF 

._ FIRE S.^CH 

NC. TTVE STEPS F£3 MOVE 
FCR ALT. POS. 1 
= = = 2 
= = = 3 



PFF 



F h y 



PKH I 
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APPENDIX E 



EXECUTIVE PROGR AM 
for the 

SMALL-UNIT AMPHIBIOUS OPERATION COMBAT MODEL 



The combat model presented in this thesis has been provided with an 
EXEC program which is designed to set uo and execute all of the necessary 
CMS commands for the running of the model. The EXEC program will auto- 
matically BROWSE the output listing of the model ( AMPHIB1 LISTING) 
allowing the user to review immediately the results of the battle. 

A listing of the EXEC follows. 

GLOBAL TXTLIB F0RTM0D2 M0D2EEH IMSLSP NONIMSL CMSLIB 

FILEDEF 05 DISK SEA DATA 

FILEDEF 09 DISK LAND DATA 

FILEDEF 06 DISK AMPHIB1 LISTING 

LOAD AMPHIB (START) 

BROWSE AMPHIB 1 LISTING 
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APPENDIX F 



COMPUTER OUTPUT 
for the 

SMALL-UNIT AMPHIBIOUS OPERATION COMBAT MODEL 

The computer output for the small -unit amphibious operation combat 
model was designed to be clear, concise, and identifiable to the user 
of the program. The combat model conducts two phases of combat: ship- 

to-shore and land combat. Therefore, the computer output was designee 
to report on each phase of combat. The computer output for each phase: 
of combat begins with an initial information page which lists the 
input data provided by the user of the model. The initial informatior 
page serves as a record of the battle scenario analyzed by tne model , 
as well as a check for the user to insure that the input data provided 
were read correctly by the model. In addition, a battle summary 
report is provided reporting on the status of both the aggressor and 
defender forces throughout both phases of combat. The computer output 
based upon the input data listed in Appendix C is as follows. 
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** INITIAL SHIF-TO-SHOFE FHASE INFOFMAT 10 N ** 

INITIAL FCFCE STRENGTH 
VAVE 12 3 4 

LVA 25.0 20.0 15. 0 10.0 



CEF. TANK ASSETS 



10.0 



5 

5.0 

OEF . AT GM ASSETS = 10.0 



LVA ENGR SPECS 

SFOMAX SPDMIN HTMAX HTMIN WIC 

4C.00 10. 5C 1.7C C.6C 2.53 

CEFENSIVE TACTICAL PARAMETERS 

RANGE A I M-R ELCAD PROJECTILE 
MAX MIN TIME VELOCITY 

TANK 1500.0 15. CO A5U.U0 

ATOM 2CCC.0 2CC.0 3C.CC 25C.C0 

CEFENSIVE TACTICAL ALLOCATION WEIGHTS: 

WAVE 1 = 2.0C WAVE 2 = l.OC 

CEFENSIVE FORCE ATTRITION COEFFICIENTS 
ALPHA * A BETA* A 

OT 0.0000* C. 00070 

DS O.OOCCE 0.C0090 

AIMEC FIFE ATTRITION RATE COEFFICIENTS FOR 
CEFENSIVE TANK AND AT GM ASSETS 

WEETAI 1 ) = O.OOC50 W BE TA < 2 > = C. COO TO 

BREAKPOINT ASSUMPTION: 0. 3* ( TOT AL DEF FORCE) 

CEFENDEF ATTRITION LEVEL ALLOWING FOR LAND CCM8AT 
0 • 22* ( TOTAL DEFENDER FORCE) 

ARTM SUP FACT 0R= 50.0 ERROR SUP FACTOR = 10C.O 

DISPERSION DATA 



RANGE 


TSIGV 


RANGE 


TSIGH 


R4NGE 


TMEANH 


25. 0 


0 . c 


25. C 


c.o 


25. 0 


0. 0 


500. C 


2.0 


5 CC . C 


2.0 


500. a 


1.0 


1 COC. C 


5. C 


10CC.C 


5.0 


1 COO.C 


5.0 


2CC0.0 


20. C 


20 CC . C 


20.0 


2 COO. 0 


10. 0 


5 COO • 0 


25. C 


5000.0 


25.0 


5 COO.C 


15.0 


10CCC.C 


25. C 


1C0CC.C 


25.0 


10C00.C 


15.0 


RANGE 


SSIGV 


RANGE 


SS 1GH 






25. C 


0. 0 


25. C 


C.O 






250. C 


5. 0 


25C .0 


5.0 






500. C 


7.5 


500. C 


7.5 






1CCC.C 


14. C 


10CC.C 


14.0 






2500. C 


15. 5 


2500 .C 


15.5 






5C0C.C 


17. C 


5000 .0 


17.0 






10C0C.C 


20. C 


100CC.C 


2C.0 






CURRENT 


STATUS 


OF WAVE 


I VARIABLE C E F IN IT I Of 



C - NOT ENGAGING 

1 - LANCED 

2 - UNCER FIRE BY A TGM 
2 - UNCER FIRE EV TANK 

4 - UNCER FIRE EY BGTH ATGM S TANK 
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4**** the ship-to-shcre phase begins ***** 

EREAKPOINT R E ACHED AT TIME = 502.5 SEC0N0S 



T I M€ s 5C2.5 SECONDS 



WVE 


FCFCE LEVEL 


STATUS 


LOST-PCT TOTAL SURVIVING 


1 


16. OQ 00 


1 


0.360 


2 


12. OCOC 


1 


0.35C 




11.0000 


1 


0 .26? 




10.0000 


1 


0 .0 


c 

TANK 


5. OCOO 
0.0 


0 


0.0 55.00 

1 .000 


ATC-M 


c.o 




1 .000 0.0 


final 


LVA SURVIVORS 


ASHORE 


= 55.CCC 


LANC 


COMBAT STARTS 


WHILE SHGRE COMBAT IS GOING ON 


LANC 


COMBAT ATTACK 


TIME = 


495. C SECONDS 
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** INITIAL LANO CCM6AT INFORMATION ** 



LOCATION 



UN IT 


X 


V 


FORCE LEVS 


1 


2000 • 0 


190C.0 


18.0 


2 


1900.0 


2400 .0 


1 8.0 


a 


15CC.0 


21 CO. 0 


1 6.0 


4 


2600.0 


2700.0 


1 0.0 


c 


36CO. C 


23 CO • 0 


5.0 


6 


2600.0 


17CC.0 


10.0 



ATTRITION IS STOCHASTIC 

ROUTES DETERMINED 6Y USER 

AT tack VEHICLE SPEED IS 18.0 M.P.H. 

6REAKPOINT DISTANCE IS 500.0 NETEFS 

DEFEN CEP WILL MOVE TO ALTERNATE POSITIONS 
ALTERNATE POSITIONS ARE: 

UNIT X Y 

A AfOO.C 3800.0 

i V5C0.C 2700 .0 

4 <4C0.0 180C.C 

ATK :ULL I* P 08 AS I L I T I ES 



RANGE 


P 


PFH 


FHM 


PKH 


50 ,0 


C.60 


0.70 


0.65 


0.85 


1 CCC 


0.65 


C .90 


0 .85 


C.90 


1 5CC 


0.80 


0.65 


C • 8 5 


C. 8C 


2 COO 


C.75 


C .80 


0.75 


0.70 


2 5CC 


0.6C 


0 • 7 C 


C .6 5 


0.65 


2 COO 


0.40 


G .4 5 


C • 4 0 


C. 5C 


DE F 


. KILL 


FFCBAei LIT! 


ES 


RANGE 


P 


PFH 


PFf* 


PKH 


5C0 


C • 35 


0.85 


0 .75 


0.70 


10CC 


0.60 


0 • 8 C 


C.75 


C • 7 0 


1 50C 


C.75 


C.75 


C.7C 


C.6C 


2 00 C 


o.eo 


C.65 


0 .60 


0 .55 


2 500 


0.45 


0 • 5 C 


0.5C 


C.35 


3C0C 


0.20 


C.2C 


0 .2C 


0.20 



CURRENT STATUS OF UNIT I VARIABLE DEFINITIONS 

0 - ALJ VE NOT FIRI NG 

1 - ALIVE AND FIRING 

2 - RULED 
2 - MCVING 



VEHICLE SPEED VAPIAELE DEFINITIONS 



1 - 


9 


MFH 


2 - 


12 


NFH 


3 - 


15 


MFH 


4 - 


18 


MFH 
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***** THE LAND COMSAT PHASE BEGINS ***** 
**** DEFENS I VE FORCE IS EL IM I MTEC. END OF BATTLE. 
T I ME = 745 ScCONCS 



AGGRESQR UNIT 
UN IT X 


INFORMATION 

V FORCE LEVEL 


STATUS 


LOST-PC T 


TARGETS 


1 2294.5 


1978.9 0 .0 


2 


1.000 




2 3144.3 


2246.3 4.8 


1 


0.731 


5 


3 3224.7 


H21.6 \7.<a 


0 


0. 005 
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1800.0 0 ,0 


2 


1.000 
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